fluent 经典问题 Fluent 收敛判断和 solver选择

合集下载

Fluent经典问题及答疑

Fluent经典问题及答疑

Fluent经典问题及答疑Fluent经典问题及答疑1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢?(#61)2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。

(13楼)3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?(#80)4 常见离散格式的性能的对比(稳定性、精度和经济性)(#62)5 在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?(#81)6 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?(#130)7 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?(#55)8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?(#56)9 在一个物理问题的多个边界上,如何协调各边界上的不同边界条件?在边界条件的组合问题上,有什么原则?10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?(#143)11 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解?(#35)12 在GAMBIT的foreground和background中,真实体和虚实体、实操作和虚操作四个之间是什么关系?13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节?(#38)14 画网格时,网格类型和网格方法如何配合使用?各种方法有什么样的应用范围及做网格时需注意的问题?(#169)15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?(#154)16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?(#40)17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?(#170)18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?(#128)19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?(#127)20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT 是怎样使用区域的?(#41)21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?(9楼)22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?(7楼)23 在FLUENT运行过程中,经常会出现“turbulence viscous rate”超过了极限值,此时如何解决?而这里的极限值指的是什么值?修正后它对计算结果有何影响?(#28)24 在FLUENT运行计算时,为什么有时候总是出现“reversed flow”?其具体意义是什么?有没有办法避免?如果一直这样显示,它对最终的计算结果有什么样的影响?(#29)25 燃烧过程中经常遇到一个“头疼”问题是计算后温度场没什么变化?即点火问题,解决计算过程中点火的方法有哪些?什么原因引起点火困难的问题? (#183)26 什么叫问题的初始化?在FLUENT中初始化的方法对计算结果有什么样的影响?初始化中的“patch”怎么理解?(12楼)27 什么叫PDF方法?FLUENT中模拟煤粉燃烧的方法有哪些?(#197)28 在利用prePDF计算时出现不稳定性如何解决?即平衡计算失败。

fluent算法的一些说明

fluent算法的一些说明

FLUENT-manual 中解算方法的一些说明,摘录翻译了其中比较重要的细节,希望对初学FLUENT的朋友在选择设置上提供一些帮助,不致走过多的弯路。

离散1、QUICK格式仅仅应用在结构化网格上,具有比second-order upwind 更高的精度,当然,FLUENT也允许在非结构网格或者混合网格模型中使用QUICK格式,在这种情况下,非结构网格单元仍然使用second-order upwind 格式计算。

2 、MUSCL格式可以应用在任何网格和复杂的3维流计算,相比second-order upwind,third-order MUSCL 可以通过减少数值耗散而提高空间精度,并且对所有的传输方程都适用。

third-order MUSCL 目前在FLUENT 中没有流态限制,可以计算诸如冲击波类的非连续流场。

3、有界中心差分格式bounded central differencing 是LES默认的对流格式,当选择LES后,所有传输方程自动转换为bounded central differencing 。

4 、low diffusion discretization 只能用在亚音速流计算,并且只适用于implicit-time,对高Mach流,或者在explicit time公式下运行LES ,必须使用 second-order upwind 。

5、改进的HRIC格式相比QUICK 与second order 为VOF计算提供了更高的精度,相比Geo-Reconstruct格式减少更多的计算花费。

6 、explicit time stepping 的计算要求苛刻,主要用在捕捉波的瞬态行为,相比implicit time stepping 精度更高,花费更少。

但是下列情况不能使用explicit time stepping:(1)分离计算或者耦合隐式计算。

explicit time stepping只能用于耦合显式计算。

FLUENT中是否收敛的判定方法

FLUENT中是否收敛的判定方法

FLUENT中判断收敛的方法
判断计算是否收敛,没有一个通用的方法。

通过残差值判断的方法,对一些问题或许很有效,但在某些问题中往往会得出错误的结论。

因此,正确的做法是,不仅要通过残差值,也要通过监测所有相关变量的完整数据,以及检查流入与流出的物质和能量是否守恒的方法来判断计算是否收敛。

1、监测残差值。

在迭代计算过程中,当各个物理变量的残差值都达到收敛标准时,计算就会发生收敛。

Fluent 默认的收敛标准是:除了能量的残差值外,当所有变量的残差值都降到低于10-3 时,就认为计算收敛,而能量的残差值的收敛标准为低于10-6。

2、计算结果不再随着迭代的进行发生变化。

有时候,因为收敛标准设置得不合适,物理量的残差值在迭代计算的过程中始终无法满足收敛标准。

然而,通过在迭代过程中监测某些代表性的流动变量,可能其值已经不再随着迭代的进行发生变化。

此时也可以认为计算收敛。

3、整个系统的质量,动量,能量都守恒。

在Flux Reports 对话框中检查流入和流出整个系统的质量,动量,能量是否守恒。

守恒,则计算收敛。

不平衡误差少于0.1%,也可以认为计算是收敛的。

谈FLUENT如何判断收敛

谈FLUENT如何判断收敛

谈FLUENT如何判断收敛1.观察点处的值不再随计算步骤的增加而变化;2.各个参数的残差随计算步数的增加而降低,最后趋于平缓;3.要满足质量守恒(计算中不牵涉到能量)或者是质量与能量守恒(计算中牵涉到能量)。

特别要指出的是,即使前两个判据都已经满足了,也并不表示已经得到合理的收敛解了,因为,如果松弛因子设置得太紧,各参数在每步计算的变化都不是太大,也会使前两个判据得到满足。

此时就要再看第三个判据了。

还需要说明的就是,一般我们都希望在收敛的情况下,残差越小越好,但是残差曲线是全场求平均的结果,有时其大小并不一定代表计算结果的好坏,有时即使计算的残差很大,但结果也许是好的,关键是要看计算结果是否符合物理事实,即残差的大小与模拟的物理现象本身的复杂性有关,必须从实际物理现象上看计算结果。

比如说一个全机模型,在大攻角情况下,解震荡得非常厉害,而且残差的量级也总下不去,但这仍然是正确的,为什么呢,因为大攻角下实际流动情形就是这样的,不断有涡的周期性脱落,流场本身就是非定常的,所以解也是波动的,处理的时候取平均就可以。

但是,如果只满足判据三和判据一,也有可能没有收敛,或者没有计算稳定,最常见的就是溶液里面的浓度扩散模拟,一开始计算几步就可能已经满足了判据三,并且过会也会满足判据一,但是残差曲线一直在降,还没有到平缓的阶段,此时说明浓度场还没有把整个计算域计算完,仅仅只是计算了连续性方程。

在report->flux里面操作,mass flow rate,把所有进出口都选上,compute一下,看看nut flux是什么水平,如果它的值小于总进口流量的1%,并且其他检测量在继续迭代之后不会发生波动,也可以认为你的解是收敛的。

1.网格质量,主要可能是相邻单元的尺寸大小相差较大,它们的尺寸之比最好控制在1.2以内,不能超过1.42.离散格式及压力速度耦合方法,如果是结构网格,建议使用高阶格式,如2阶迎风格式等,如果是非结构网格,除pressure保持standard格式不变外,其他格式改用高阶格式;压力速度耦合关系,如果使用SIMPLE,SIMPLEC,PISO等segerated solver对联系方程收敛没有提高的话,可以尝试使用coupled solver。

三十二、Fluent收敛判断标准及方法

三十二、Fluent收敛判断标准及方法

三十二、Fluent收敛判断标准及方法1. Fluent迭代计算流程1.1 计算流程在讲述收敛之前,我们先了解一下Fluent是如何进行迭代求解的。

在建模划分网格后,通过一系列设置,单击Calculation,Fluent 开始计算,其流程如下:1.根据初始化值,假设流场物理量初值为Q02.在初值基础上,根据守恒方程,对流场进行计算,得出新的流场物理量值Q13.比较Q0与Q1的值,|Q1-Q0|称为绝对误差,|Q1-Q0|/Q1为相对误差。

4.若误差小于设定值P,则达到收敛;否则,未达到收敛,取新的迭代值Q2进行计算,重复以上步骤。

1.2 残差的理解检测的的物理量主要有速度、质量、能量、湍流参数等,如上图1 这里的误差并不是残差Residual,实际上残差的计算要复杂的多,但是为了理解方便,我们可以认为残差等同于误差。

下图是残差公式:并不是所有的工况都会收敛。

当各物理量的值基本不变时,即残差很小时,工况才可能收敛。

但是对于瞬态,各物理量的值总是变化,如何收敛?正因如此瞬态才有时间步的概念,瞬态问题在每个时间步上都认为是稳态,所以瞬态问题的残差图总是波浪线型。

2. 收敛标准2.1 残差标准对所有的工况,没有统一的判断标准。

对于大多数问题,默认的判断标准已经足够(For most problems, the default convergence criterion in ANSYS Fluent is sufficient. )建议残差达到设定值后,多算50步,确定残差之后都是减小的趋势。

2.2 监测物理量为了使结果更加精确,可根据工况辅助检测一些物理量,如速度,流量等。

当所关心的物理量基本不变时,说明达到了收敛。

某些情况下即使残差没有达到设定值,只要所检测的物理量很稳定,也可以认为收敛2.3 通量守恒除了上述两种情况,还可以通过Flux守恒来判断。

当计算完成后,通过查看Flux是否守恒来判断是否收敛。

fluent收敛判别

fluent收敛判别

这个1e-3或者1e-4的收敛标准是相对而言的。

在FLUENT中残差是以开始5步的平均值为基准进行比较的。

如果你的初值取得好,你的迭代会很快收敛,但是你的残差却依然很高;但是当你改变初场到比较不同的值时,你的残差开始会很大,但随后却可以很快降低到很低的水平,让你看起来心情很好。

其实两种情况下流场是基本相同的。

由此来看,判断是否收敛并不是严格根据残差的走向而定的。

可一般来说,压力的收敛相对比较慢一些的。

是否收敛不能简单看残差图,还有许多其他的重要标准,比如进出口流量差、压力系数波动等等尽管残差仍然维持在较高数值,但凭其他监测也可判断是否收敛。

最重要的就是是否符合物理事实或试验结论。

残差曲线是否满足只是一个表面的现象,还要看进口和出口总量差不得大于1%,而且即使这样子,收敛解也不一定准确,它和网格划分/离散化误差,以及屋里模型的准确性都有关系.所以得有试验数据做对比活着理论分析了当然最终是否正确是要看是否与实验数据相符合!但既然有残差图的话,总应该可以大概的看出是否收敛吧?是否要残差要小到敛的残差的大小不能决定是否收敛,我在用FLUENT计算时,多采用监测一个面的速度(或者是压力、紊动能等参数)基本上不随着计算时间的推移而变化,就认为基本达到收敛4 据质量守恒,收敛时进、出口的流量数值应大致相等(一般认为进出口质量差值比上入口质量的相对值小于0.5%时收敛,但是对特殊情况可能不同),但符号相反,一般出口流量是负值。

5 在进行稳态计算时候,开始残差线是一直下降的,可是到后来各种残差线都显示为波形波动,是不是不收敛阿?答:的问题就能收敛,但如果网格质量不好,是很难的。

通常,计算非结构网格,如果问题比较复杂,会出现这种情况,建议作网格时多下些功夫。

理论上说,残差的震荡是数值迭代在计算域内传递遭遇障碍物反射形成周期震荡导致的结果,与网格亚尺度雷诺数有关。

例如,通常压力边界是主要的反射源,换成OUTFLOW边界会好些。

FLUENT操作过程及参数选择

FLUENT操作过程及参数选择

振动流化床仿真操作过程及参数选择1创建流化床模型。

根据靳海波论文提供的试验机参数,创建流化床模型。

流化床直148mm,高1m,开孔率9%,孔径2mm。

在筛板上铺两层帆布保证气流均布。

因为实验机为一个圆形的流化床,所以可简化为仅二维模型。

而实际实验中流化高度远小于1m,甚至500mm,所以为提高计算时间,可将模型高度缩为500mm。

由于筛板上铺设两层帆布以达到气流均分的目的,所以认为沿整个筛板的进口风速为均匀的。

最终简化模型如下图所示:上图为流化后的流化床模型,可以看出流化床下端的网格相对上端较密,因为流化行为主要发生的流化床下端,为了加快计算时间,所以采用这种下密上疏的划分方式。

其中进口设置为velocity inlet;出口设置为outflow;左右两边分为设置为wall。

在GAMBIT中设置完毕后,输出二维模型vfb.msh。

outflow边界条件不需要给定任何入口的物理条件,但是应用也会有限制,大致为以下四点:1.只能用于不可压缩流动2.出口处流动充分发展3.不能与任何压力边界条件搭配使用(压力入口、压力出口)4.不能用于计算流量分配问题(比如有多个出口的问题)2打开FLUENT 6.3.26,导入模型vfb.msh点击GRID—CHECK,检查网格信息及模型中设置的信息,核对是否正确,尤其查看是否出现负体积和负面积,如出现马上修改。

核对完毕后,点击GRID-SCALE 弹出SCALE GRID窗口,设置单位为mm,并点击change length unit按钮。

具体设置如下:3设置求解器保持其他设置为默认,更改TIME为unsteady,因为实际流化的过程是随时间变化的。

(1)pressure based 求解方法在求解不可压流体时,如果我们联立求解从动量方程和连续性方程离散得到的代数方程组,可以直接得到各速度分量及相应的压力值,但是要占用大量的计算内存,这一方法已可以在Fluent6.3中实现,所需内存为分离算法的1.5-2倍。

FLUENT操作过程及参数选择

FLUENT操作过程及参数选择

振动流化床仿真操作过程及参数选择1创建流化床模型。

根据靳海波论文提供的试验机参数,创建流化床模型。

流化床直148mm,高1m,开孔率9%,孔径2mm。

在筛板上铺两层帆布保证气流均布。

因为实验机为一个圆形的流化床,所以可简化为仅二维模型。

而实际实验中流化高度远小于1m,甚至500mm,所以为提高计算时间,可将模型高度缩为500mm。

由于筛板上铺设两层帆布以达到气流均分的目的,所以认为沿整个筛板的进口风速为均匀的。

最终简化模型如下图所示:上图为流化后的流化床模型,可以看出流化床下端的网格相对上端较密,因为流化行为主要发生的流化床下端,为了加快计算时间,所以采用这种下密上疏的划分方式。

其中进口设置为velocity inlet;出口设置为outflow;左右两边分为设置为wall。

在GAMBIT中设置完毕后,输出二维模型vfb.msh。

outflow边界条件不需要给定任何入口的物理条件,但是应用也会有限制,大致为以下四点:1.只能用于不可压缩流动2.出口处流动充分发展3.不能与任何压力边界条件搭配使用(压力入口、压力出口)4.不能用于计算流量分配问题(比如有多个出口的问题)2打开FLUENT 6.3.26,导入模型vfb.msh点击GRID—CHECK,检查网格信息及模型中设置的信息,核对是否正确,尤其查看是否出现负体积和负面积,如出现马上修改。

核对完毕后,点击GRID-SCALE 弹出SCALE GRID窗口,设置单位为mm,并点击change length unit按钮。

具体设置如下:3设置求解器保持其他设置为默认,更改TIME为unsteady,因为实际流化的过程是随时间变化的。

(1)pressure based 求解方法在求解不可压流体时,如果我们联立求解从动量方程和连续性方程离散得到的代数方程组,可以直接得到各速度分量及相应的压力值,但是要占用大量的计算内存,这一方法已可以在Fluent6.3中实现,所需内存为分离算法的1.5-2倍。

fluent 收敛判断问题

fluent 收敛判断问题

FLUENT 收敛判断问题及解决方法1 Under-Relaxation Factors(亚松弛因子)所谓亚松弛,就是将本层次计算结果与上一层次结果的差值作适当缩减,以避免由于差值过大而引起非线性迭代过程的发散。

用通用变量来写出时,为亚松弛因子(Relaxation Factors)。

FLUENT中的亚松弛:由于FLUENT所解方程组的非线性,我们有必要控制变化。

一般用亚松弛方法来实现控制,该方法在每一部迭代中减少了变化量。

亚松弛最简单的形式为:单元内变量等于原来的值加上亚松弛因子a与变化的积。

分离解算器使用亚松弛来控制每一步迭代中的计算变量的更新。

这就意味着使用分离解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个相关的亚松弛因子。

在FLUENT中,所有变量的默认亚松弛因子都是对大多数问题的最优值。

这个值适合于很多问题,但是对于一些特殊的非线性问题(如某些湍流或者高Raleigh数自然对流问题),在计算开始时要慎重减小亚松弛因子。

使用默认的亚松弛因子开始计算是很好的习惯,如果经过4到5步的迭代残差仍然增长,你就需要减小亚松弛因子。

有时候,如果发现残差开始增加,你可以改变亚松弛因子重新计算。

在亚松弛因子过大时通常会出现这种情况。

最为安全的方法就是在对亚松弛因子做任何修改之前先保存数据文件,并对解的算法做几步迭代以调节到新的参数。

最典型的情况是,亚松弛因子的增加会使残差有少量的增加,但是随着解的进行残差的增加又消失了。

如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。

注意:粘性和密度的亚松弛是在每一次迭代之间的。

而且,如果直接解焓方程而不是温度方程(即对PDF计算),基于焓的温度的更新是要进行亚松驰的。

要查看默认的亚松弛因子的值,你可以在解控制面板点击默认按钮。

对于大多数流动,不需要修改默认亚松弛因子。

但是,如果出现不稳定或者发散你就需要减小默认的亚松弛因子了,其中压力、动量、k和e的亚松弛因子默认值分别为0.2,0.5,0.5和0.5。

fluent中判断收敛的方法

fluent中判断收敛的方法

fluent中判断收敛的方法[引用]FLUENT中判断收敛的方法判断计算是否收敛,没有一个通用的方法。

通过残差值判断的方法,对一些问题或许很有效,但在某些问题中往往会得出错误的结论。

因此,正确的做法是,不仅要通过残差值,也要通过监测所有相关变量的完整数据,以及检查流入与流出的物质和能量是否守恒的方法来判断计算是否收敛。

1、监测残差值。

在迭代计算过程中,当各个物理变量的残差值都达到收敛标准时,计算就会发生收敛。

Fluent默认的收敛标准是:除了能量的残差值外,当所有变量的残差值都降到低于10-3 时,就认为计算收敛,而能量的残差值的收敛标准为低于10-6。

2、计算结果不再随着迭代的进行发生变化。

有时候,因为收敛标准设置得不合适,物理量的残差值在迭代计算的过程中始终无法满足收敛标准。

然而,通过在迭代过程中监测某些代表性的流动变量,可能其值已经不再随着迭代的进行发生变化。

此时也可以认为计算收敛。

3、整个系统的质量,动量,能量都守恒。

在Flux Reports对话框中检查流入和流出整个系统的质量,动量,能量是否守恒。

守恒,则计算收敛。

不平衡误差少于0.1%,也可以认为计算是收敛的。

FLUENT中残差的概念残差是cell各个face的通量之和,当收敛后,理论上当单元内没有源项使各个面流入的通量也就是对物理量的输运之和应该为零。

最大残差或者RSM残差反映流场与所要模拟流场(只收敛后应该得到的流场,当然收敛后得到的流场与真实流场之间还是存在一定的差距)的残差,残差越小越好,由于存在数值精度问题,不可能得到0残差,对于单精度计算一般应该低于初始残差1e-03以下才好,当注意具体情况,看各个项的收敛情况(比方说连续项不易收敛而能量项容易)。

一般在FLUENT中可以进行进出口流量监控,当残差收敛到一定程度后,还要看进出口流量是否稳定平衡,才可确定收敛与否(翼型计算时要监控升阻力的平衡)。

残差在较高位震荡,需要检查边界条件是否合理,其次检查初始条件是否合理,比如激波的流场,初始条件的不合适会造成流场的振荡。

fluent 收敛判断问题

fluent 收敛判断问题

FLUENT 收敛判断问题及解决方法1 Under-Relaxation Factors(亚松弛因子)所谓亚松弛,就是将本层次计算结果与上一层次结果的差值作适当缩减,以避免由于差值过大而引起非线性迭代过程的发散。

用通用变量来写出时,为亚松弛因子(Relaxation Factors)。

FLUENT中的亚松弛:由于FLUENT所解方程组的非线性,我们有必要控制变化。

一般用亚松弛方法来实现控制,该方法在每一部迭代中减少了变化量。

亚松弛最简单的形式为:单元内变量等于原来的值加上亚松弛因子a与变化的积。

分离解算器使用亚松弛来控制每一步迭代中的计算变量的更新。

这就意味着使用分离解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个相关的亚松弛因子。

在FLUENT中,所有变量的默认亚松弛因子都是对大多数问题的最优值。

这个值适合于很多问题,但是对于一些特殊的非线性问题(如某些湍流或者高Raleigh数自然对流问题),在计算开始时要慎重减小亚松弛因子。

使用默认的亚松弛因子开始计算是很好的习惯,如果经过4到5步的迭代残差仍然增长,你就需要减小亚松弛因子。

有时候,如果发现残差开始增加,你可以改变亚松弛因子重新计算。

在亚松弛因子过大时通常会出现这种情况。

最为安全的方法就是在对亚松弛因子做任何修改之前先保存数据文件,并对解的算法做几步迭代以调节到新的参数。

最典型的情况是,亚松弛因子的增加会使残差有少量的增加,但是随着解的进行残差的增加又消失了。

如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。

注意:粘性和密度的亚松弛是在每一次迭代之间的。

而且,如果直接解焓方程而不是温度方程(即对PDF计算),基于焓的温度的更新是要进行亚松驰的。

要查看默认的亚松弛因子的值,你可以在解控制面板点击默认按钮。

对于大多数流动,不需要修改默认亚松弛因子。

但是,如果出现不稳定或者发散你就需要减小默认的亚松弛因子了,其中压力、动量、k和e的亚松弛因子默认值分别为0.2,0.5,0.5和0.5。

FLUENT操作过程及参数选择

FLUENT操作过程及参数选择

振动流化床仿真操作过程及参数选择1创建流化床模型。

根据靳海波论文提供的试验机参数,创建流化床模型。

流化床直148mm,高1m,开孔率9%,孔径2mm。

在筛板上铺两层帆布保证气流均布。

因为实验机为一个圆形的流化床,所以可简化为仅二维模型。

而实际实验中流化高度远小于1m,甚至500mm,所以为提高计算时间,可将模型高度缩为500mm。

由于筛板上铺设两层帆布以达到气流均分的目的,所以认为沿整个筛板的进口风速为均匀的。

最终简化模型如下图所示:上图为流化后的流化床模型,可以看出流化床下端的网格相对上端较密,因为流化行为主要发生的流化床下端,为了加快计算时间,所以采用这种下密上疏的划分方式。

其中进口设置为velocity inlet;出口设置为outflow;左右两边分为设置为wall。

在GAMBIT中设置完毕后,输出二维模型vfb.msh。

outflow边界条件不需要给定任何入口的物理条件,但是应用也会有限制,大致为以下四点:1.只能用于不可压缩流动2.出口处流动充分发展3.不能与任何压力边界条件搭配使用(压力入口、压力出口)4.不能用于计算流量分配问题(比如有多个出口的问题)2打开FLUENT 6.3.26,导入模型vfb.msh点击GRID—CHECK,检查网格信息及模型中设置的信息,核对是否正确,尤其查看是否出现负体积和负面积,如出现马上修改。

核对完毕后,点击GRID-SCALE 弹出SCALE GRID窗口,设置单位为mm,并点击change length unit按钮。

具体设置如下:3设置求解器保持其他设置为默认,更改TIME为unsteady,因为实际流化的过程是随时间变化的。

(1)pressure based 求解方法在求解不可压流体时,如果我们联立求解从动量方程和连续性方程离散得到的代数方程组,可以直接得到各速度分量及相应的压力值,但是要占用大量的计算内存,这一方法已可以在Fluent6.3中实现,所需内存为分离算法的1.5-2倍。

Fluent 中判断收敛的方法、残差的概念及不收敛通常的解决方式

Fluent 中判断收敛的方法、残差的概念及不收敛通常的解决方式

fluent中判断收敛的方法[引用]FLUENT中判断收敛的方法判断计算是否收敛,没有一个通用的方法。

通过残差值判断的方法,对一些问题或许很有效,但在某些问题中往往会得出错误的结论。

因此,正确的做法是,不仅要通过残差值,也要通过监测所有相关变量的完整数据,以及检查流入与流出的物质和能量是否守恒的方法来判断计算是否收敛。

1、监测残差值。

在迭代计算过程中,当各个物理变量的残差值都达到收敛标准时,计算就会发生收敛。

Fluent默认的收敛标准是:除了能量的残差值外,当所有变量的残差值都降到低于10-3 时,就认为计算收敛,而能量的残差值的收敛标准为低于10-6。

2、计算结果不再随着迭代的进行发生变化。

有时候,因为收敛标准设置得不合适,物理量的残差值在迭代计算的过程中始终无法满足收敛标准。

然而,通过在迭代过程中监测某些代表性的流动变量,可能其值已经不再随着迭代的进行发生变化。

此时也可以认为计算收敛。

3、整个系统的质量,动量,能量都守恒。

在Flux Reports对话框中检查流入和流出整个系统的质量,动量,能量是否守恒。

守恒,则计算收敛。

不平衡误差少于0.1%,也可以认为计算是收敛的。

FLUENT中残差的概念残差是cell各个face的通量之和,当收敛后,理论上当单元内没有源项使各个面流入的通量也就是对物理量的输运之和应该为零。

最大残差或者RSM残差反映流场与所要模拟流场(只收敛后应该得到的流场,当然收敛后得到的流场与真实流场之间还是存在一定的差距)的残差,残差越小越好,由于存在数值精度问题,不可能得到0残差,对于单精度计算一般应该低于初始残差1e-03以下才好,当注意具体情况,看各个项的收敛情况(比方说连续项不易收敛而能量项容易)。

一般在FLUENT中可以进行进出口流量监控,当残差收敛到一定程度后,还要看进出口流量是否稳定平衡,才可确定收敛与否(翼型计算时要监控升阻力的平衡)。

残差在较高位震荡,需要检查边界条件是否合理,其次检查初始条件是否合理,比如激波的流场,初始条件的不合适会造成流场的振荡。

fluent算法的一些说明

fluent算法的一些说明

fluent算法的⼀些说明FLUENT-manual 中解算⽅法的⼀些说明,摘录翻译了其中⽐较重要的细节,希望对初学FLUENT的朋友在选择设置上提供⼀些帮助,不致⾛过多的弯路。

离散1、 QUICK格式仅仅应⽤在结构化⽹格上,具有⽐second-order upwind 更⾼的精度,当然,FLUENT也允许在⾮结构⽹格或者混合⽹格模型中使⽤QUICK格式,在这种情况下,⾮结构⽹格单元仍然使⽤second-order upwind 格式计算。

2 、MUSCL格式可以应⽤在任何⽹格和复杂的3维流计算,相⽐second-order upwind,third-order MUSCL 可以通过减少数值耗散⽽提⾼空间精度,并且对所有的传输⽅程都适⽤。

third-order MUSCL ⽬前在FLUENT中没有流态限制,可以计算诸如冲击波类的⾮连续流场。

3、有界中⼼差分格式bounded central differencing 是LES默认的对流格式,当选择 LES后,所有传输⽅程⾃动转换为bounded central differencing 。

4 、low diffusion discretization 只能⽤在亚⾳速流计算,并且只适⽤于implicit-time,对⾼Mach流,或者在explicit time公式下运⾏LES ,必须使⽤ second-order upwind 。

5、改进的HRIC格式相⽐QUICK 与second order 为VOF计算提供了更⾼的精度,相⽐Geo-Reconstruct格式减少更多的计算花费。

6 、explicit time stepping 的计算要求苛刻,主要⽤在捕捉波的瞬态⾏为,相⽐implicit time stepping 精度更⾼,花费更少。

但是下列情况不能使⽤explicit time stepping:(1)分离计算或者耦合隐式计算。

explicit time stepping只能⽤于耦合显式计算。

fluent收敛

fluent收敛

计算不收敛的话可能是什么原因呢如果经过长时间不收敛有两种情况:(1) 由于网格质量不好引起,这需要改进网格;或者(2) 也有可能已经收敛,但残差曲线并没有下降到要求的标准,此时需要检验特征点上的求解变量的变化,如果很小(达到自己的要求),可以认为已经收敛。

如果开始迭代次数不多即发散:(1) 网格质量不好,需要改进网格,尤其要找到发散点的位置,重点改进该处的网格质量;或者(2) 调整导致发散的初始变量的值(通常为k,e和能量),甚至有关方程的松弛因子,使其绕过发散区间,继续运行。

如果仍然发散,你所选用的物理和化学模型稳定性太差,建议在牺牲先进性的前提下,改用其他成熟的模型。

在fluent里,三角形网格的质量一般要小于0.8,计算结果才容易收敛。

转载:利用FLUENT不收敛通常怎么解决?①、一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛。

②、FLUENT的收敛最基础的是网格的质量,计算的时候看怎样选择CFL数,这个靠经验③、首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。

④、有时初始条件和边界条件严重影响收敛性,曾经作过一个计算反反复复,通过修改网格,重新定义初始条件,包括具体的选择的模型,还有老师经常用的方法就是看看哪个因素不收敛,然后寻找和它有关的条件,改变相应参数。

就收敛了⑤、A.检查是否哪里设定有误:比方用mm的unit建构的mesh,忘了scale;比方给定的边界条件不合理。

B从算至发散前几步,看presure分布,看不出来的话,再算几步, 看看问题大概出在那个区域。

C网格,配合第二点作修正,就重建个更漂亮的,或是更粗略的来处理。

D再找不出来的话,换个solver。

⑥、解决的办法是设几个监测点,比如出流或参数变化较大的地方,若这些地方的参数变化很小,就可以认为是收敛了,尽管此时残值曲线还没有降下来。

⑦、调节松弛因子也能影响收敛,不过代价是收敛速度。

fluent经典问题之一

fluent经典问题之一
本文版权归原作者所有。版权所有者对本文转移、存储、 复制等方式提出异议,同时援引有效、适用的法律、法规条款, westlife 及个人主义社区愿意立即终止为其本文的发布、传播 而所做出的一切形式努力。
注:个人主义社区唯一指定为域名,其它域名与个人主义社区无 关。
QQ:9184482(请注明:版权信息) E-mail:wust#(请把#改成@)

个人主义社区声明
本文由westlife 从网上搜集而来,仅为在教育与科研领域 从事研究工作的人员提供参考与帮助,无意于将其用于其它目 的。
对本文在个人主义社区并通过个人主义社区在网络的转 移、存储、复制,westlife 及个人主义社区不提供基有损于原 文版权所有者的利益、形象等权益的帮助或便利。
我们经常遇到计算区是对称的问题,如同心圆环内的自然对流,圆柱绕流,我们为了节 省计算资源,许多时候都把计算区域趣味一半,但有些问题的真实情况是 两步分的流场及 物理量的分布并不对称呀,问我们如何判断该不该区一般呢?
对秤的问题一般用在流场稳态解..需满足 1.几何图形对秤..2.边界条件对秤.. 也就是物理条件对秤..3.structral 网格..所以对秤轴的 Flux 和 properties gradient 必须为 0... 3 按照算例学习了一段时间,有些简单的问题还可以分析对错,但是对于一些头脑里没 有概念的问题,是做出了很多图行了,矢量图了,但是如何比较仔细全面的分析其合理性, 觉得有些困难,望师兄指点~~~~~~~~~~~~
答: 一般来讲计算应该辅助以高精度的实验作为证明,无法或不容易用实验实现的往往是计算 一个经典的或别人算过的例子对比一下。既然已经作出了很多图,可以试试分析一下, 看看跟经典的理论一致否。图是做出来了,但是真的想说明问题,恐怕还有一段路。

fluent 经典问题 Fluent 收敛判断和 solver选择

fluent 经典问题 Fluent 收敛判断和 solver选择

Fluent 收敛判断和solver选择从上转过来的,虽然是英语,但是静下心来慢慢读一读能学到很多问题:----------------------------Hi! I have tried an external aerodynamic problem in Flunet. In it, I want to know pressure distribution over the wing assembly.I have used Coupled-Implicit-Spalart Allamaras solver with courant number 1 initially. I gave pressure-far-field BC in elliptical boundary around wing assembly which is 10 times larger.After 5000 iterations also, my solution is not converging or continuity and momentum residuals are not coming below 1e-3. They oscillating between 1e-2 and 1e-3. Viscous residual is less than 1e-3.I have changed under relaxation factors, discretization schemes also. Still, I am not able to achieve residual lesser than 1e-4.I want any one users help. I am conveying my advance thanks ..........with regards, vengi.某人的回答There's a few things that could be going on.One possible answer is that your model is converged (that's always the happiest answer, isn't it?). The residuals you are looking at are normalized based on the residuals of the first iteration. So if your initial guess is pretty accurate, then your first residuals will be small, and all of your following residuals will be small as well, but since they are normalized according to that first small value, they look large. This typically shows up in the continuity and momentum residuals, and sometimes even in the x, y, and z velocity residuals (at least in the coupled solver). One thing you should be doing with your model is monitoring other factors besides your residuals. If you're looking for the pressure distribution, then define a few points along your airfoil and monitor the pressure at these points. You should also monitor at least the lift of your airfoil. You can find these monitors undersolve->monitors. Judge convergence by when these have leveled off. While your model is solving, you will probably have to go in and clear the data in the monitors or adjust the scale of the axis to get a better idea of when they've truly leveled off. That can all be done in the windows where you defined the monitors.Another possibility is that your model isn't converged (the less happy of the answers). If that's the case, then there's lots of possible reasons. One common one is the use of the Coupled Solver in low speed flows. Since the coupled solveris a density based solver, it can get hung up in incompressible flow regimes. Typically, I only use the coupled solver for flows over Mach 0.7, but I've used the segregated solver from Mach 0.05 up to Mach 1.2 (paying CAREFUL attention to the mesh where shocks form). Another possible problem is that its an unsteady problem. If you've stalled, you could be shedding vortices at some frequency. The SA turb model does alright with small separation regions, but a large separation region (say behind a shock at some angle of attack) can cause it to fall apart. It was originally designed for 2D airfoils without any separation. They've modified it some to try and make it work in 3D, and to try and help it handle separation, but I still haven't had much luck with it. There could also be some issues with your mesh. Pay attention to your y+ values and the rules concerning them.Either way, you really should be monitoring more than the residuals to judge convergence. I've seen it a lot here, where someone will call a model converged because the residuals dropped below 1e-03, but when I've taken the model and continued with the iterations, I've seen a dramatic change in the forces. I've also seen it where someone will be 8 or 9 thousand iterations in trying to get the residuals to drop, but the forces have been steady.Hope this helps, and good luck,Jason又有一个人来提问Hi Jason, could you tell me more about using the segregated for transonic on ligthly supersonic flows? I tried it with the AGARD 445.6 wing for flutter determination and I had very good results for transonic flows. So, If we take care of under-relaxation factors and we make a good mesh, can segregated solver be used even for transonic flows?lucaSysnoise's Users_Manual接着回答You can use the segregated solver for transonic flows. It tends to diffuse the shocks compared to the coupled solver, but you can fix that by using a refined mesh. A lot of times that's what I have to do here, because I'm running into memory limits, and increasing the mesh by 10-20% still fits in the available memory, where switching to coupled solver doesn't (coupled solver uses 1.5x to 2x the memory because it stores the solution to the last iteration).As far as solver settings... it's very important to set your control limits. For temp and pressure, I calculate the delta between freestream and stagnation, and Idouble that. If I have a problem getting the model to converge, I may cut that down to about 1.5 to 1.2 times the difference between freestream and stagnation. (So my Pressure limits are Pstatic + 2*Q and Pstatic - 2*Q) I set the pressure and momentum URFs to 0.5 and 0.4, leaving energy at 1. The SIMPLEPressure-Velocity coupling tends to work well... I've gotten some recommendations on switching to PISO, so I'm actually trying that right now. I tend to run for 10 to 50 iterations with all the default discretizations and the turb eq turned off. Then I turn the turb eq back on (I don't go all the way to convergence like some people recommend... I haven't found any benefit... I usually only do about 50 iterations with these settings). After that, I set everything to 2nd order discretization and run to convergence. These settings have worked well for me on aerodynamic models with little or no separation, and they've worked (with a little playing around in the URFs) all the way to Mach 1.2.I have run bluff bodies at subsonic compressible to transonic speeds (typically a symmetrical model to get more of a "time averaged" solution... this avoids the oscillating vortices and cuts out the need for the unsteady solver... I haven't run transonic of a full bluff body, but I have run Mach 0.5ish with a full model and the unsteady solver... that was a while ago though, and I don't remember if I had made any changes to my solver settings). The coupled solver is a poor solver for this kind of model because of the large separation region aft of the body. This becomes a difficult model for the segregated solver as well, but I've had good luck running it with default discretization and the turb model turned off for about 100 to 200 iterations... then turning on the turb model, and lowering the energy URF to about 0.7 and running for another 100 or so iterations. Then upping the energy URF back to 1 for another 50 or so, and then switching to second order and running to convergence.The most important things I have found are paying attention to your mesh (or using adaption or dynamic adaption to resolve shocks... you don't need to refine them all the way for overall forces, but you might for flutter analysis) and setting your control limits.Hope this helps, and good luck,JasonThank you Jason, your explanations are complete and well detailed as usual. I run the Agard 445.6 at mach 1.141 and had no problem with convergence. I just set momentum URF to 0.4. Yes you're rigth when you say coupled solver requires more memory. In fact if I use the segregated I can make a more refined grid because it requires less memory. This is great! SIMPLE scheme seems to work well and I didn't need to switch to PISO...and I had no problem to solve the flow with transonic or supersonic flow. Thank you again for your answer. I just needed to have a confirm somebody else tried to use the segregated solver for not subsonic external flows. Luca1.我一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛~2.FLUENT的收敛最基础的是网格的质量,计算的时候看怎样选择CFL数,这个*经验3.首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。

如何解决FLUENT计算不收敛

如何解决FLUENT计算不收敛

如何解决FLUENT计算不收敛FLUENT专题2009-09-15 09:26:34 阅读261 评论1 字号:大中小①、一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛。

②、FLUENT的收敛最基础的是网格的质量,计算的时候看怎样选择CFL数,这个靠经验③、首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。

④、有时初始条件和边界条件严重影响收敛性,曾经作过一个计算反反复复,通过修改网格,重新定义初始条件,包括具体的选择的模型,还有老师经常用的方法就是看看哪个因素不收敛,然后寻找和它有关的条件,改变相应参数。

就收敛了⑤、A.检查是否哪里设定有误:比方用mm的unit建构的mesh,忘了scale;比方给定的边界条件不合理。

B从算至发散前几步,看presure分布,看不出来的话,再算几步, 看看问题大概出在那个区域。

C网格,配合第二点作修正,就重建个更漂亮的,或是更粗略的来处理。

D再找不出来的话,换个solver。

⑥、解决的办法是设几个监测点,比如出流或参数变化较大的地方,若这些地方的参数变化很小,就可以认为是收敛了,尽管此时残值曲线还没有降下来。

⑦、调节松弛因子也能影响收敛,不过代价是收敛速度。

可压流动-1在高速流动或者压力变化剧烈的流动中会遇到压缩性影响.当流体流动速度接近或者大于气体的声速时或者系统的压力变化D p /p较大时,气体随压力变化密度会对流动速度,温度有较大的影响.可压缩性流动具有独特的物理性质当马赫数远小于1.0时(M < 0.1),可压缩性影响可以忽略,在模拟中气体密度随压力的变化可以忽略.当马赫数接近1.0时(被称为跨声速流动区域),可压缩性影响变得越来越重要了.当马赫>1.0时,称为超声速流动,其中可能会包含激波或者膨胀波,它们会对流动有较大的影响.可压缩性流动最为典型的是使用总压p_0和总温T_0来描述.一、可压流动的基本方程:可压缩性流动用Fluent所解的标准连续性和动量方程来描述,你不必设定任何特殊的物理模型(除了气体定律的可压缩性形式一节中介绍的密度的可压缩性处理.FLUENT所解的能量方程很好的处理了流动速度和静温之间的耦合,不管你什么时候解可压缩性流动都必须激活能量方程.除此之外,如果你使用分离解算器,你需要激活能量方程一节中方程的粘性耗散项,该项在高马赫数流动中会变得很重要.二、在FLUENT中设定可压流动步骤如下:1. 在操作条件面板中设定操作压力@Define/Operating Conditions....(你可以认为p_op为流动中标准压力p为零的点处的绝对静压.)2. 在能量面板中激活能量方程的解@Define/Models Energy....3. 如果模拟湍流流动,请在粘性面板打开粘性热传导来激活能量方程的粘性耗散项(只用于分离解算器).这些项在高速流动中可能会十分重要@Define/Models/Viscous....如果使用耦合解算器这一步就不必要了,因为耦合解算器在能量方程中总是包含粘性耗散项.4. 在使用材料面板中设定@Define/Materials...(1). 在密度后面的下拉菜单中选择理想气体(2). 定义所有的相关属性(比热,分子量,热传导系数等).5. 设定边界条件,要确保边界条件与流动区域很好的结合。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Fluent 收敛判断和solver选择从上转过来的,虽然是英语,但是静下心来慢慢读一读能学到很多问题:----------------------------Hi! I have tried an external aerodynamic problem in Flunet. In it, I want to know pressure distribution over the wing assembly.I have used Coupled-Implicit-Spalart Allamaras solver with courant number 1 initially. I gave pressure-far-field BC in elliptical boundary around wing assembly which is 10 times larger.After 5000 iterations also, my solution is not converging or continuity and momentum residuals are not coming below 1e-3. They oscillating between 1e-2 and 1e-3. Viscous residual is less than 1e-3.I have changed under relaxation factors, discretization schemes also. Still, I am not able to achieve residual lesser than 1e-4.I want any one users help. I am conveying my advance thanks ..........with regards, vengi.某人的回答There's a few things that could be going on.One possible answer is that your model is converged (that's always the happiest answer, isn't it?). The residuals you are looking at are normalized based on the residuals of the first iteration. So if your initial guess is pretty accurate, then your first residuals will be small, and all of your following residuals will be small as well, but since they are normalized according to that first small value, they look large. This typically shows up in the continuity and momentum residuals, and sometimes even in the x, y, and z velocity residuals (at least in the coupled solver). One thing you should be doing with your model is monitoring other factors besides your residuals. If you're looking for the pressure distribution, then define a few points along your airfoil and monitor the pressure at these points. You should also monitor at least the lift of your airfoil. You can find these monitors undersolve->monitors. Judge convergence by when these have leveled off. While your model is solving, you will probably have to go in and clear the data in the monitors or adjust the scale of the axis to get a better idea of when they've truly leveled off. That can all be done in the windows where you defined the monitors.Another possibility is that your model isn't converged (the less happy of the answers). If that's the case, then there's lots of possible reasons. One common one is the use of the Coupled Solver in low speed flows. Since the coupled solveris a density based solver, it can get hung up in incompressible flow regimes. Typically, I only use the coupled solver for flows over Mach 0.7, but I've used the segregated solver from Mach 0.05 up to Mach 1.2 (paying CAREFUL attention to the mesh where shocks form). Another possible problem is that its an unsteady problem. If you've stalled, you could be shedding vortices at some frequency. The SA turb model does alright with small separation regions, but a large separation region (say behind a shock at some angle of attack) can cause it to fall apart. It was originally designed for 2D airfoils without any separation. They've modified it some to try and make it work in 3D, and to try and help it handle separation, but I still haven't had much luck with it. There could also be some issues with your mesh. Pay attention to your y+ values and the rules concerning them.Either way, you really should be monitoring more than the residuals to judge convergence. I've seen it a lot here, where someone will call a model converged because the residuals dropped below 1e-03, but when I've taken the model and continued with the iterations, I've seen a dramatic change in the forces. I've also seen it where someone will be 8 or 9 thousand iterations in trying to get the residuals to drop, but the forces have been steady.Hope this helps, and good luck,Jason又有一个人来提问Hi Jason, could you tell me more about using the segregated for transonic on ligthly supersonic flows? I tried it with the AGARD 445.6 wing for flutter determination and I had very good results for transonic flows. So, If we take care of under-relaxation factors and we make a good mesh, can segregated solver be used even for transonic flows?lucaSysnoise's Users_Manual接着回答You can use the segregated solver for transonic flows. It tends to diffuse the shocks compared to the coupled solver, but you can fix that by using a refined mesh. A lot of times that's what I have to do here, because I'm running into memory limits, and increasing the mesh by 10-20% still fits in the available memory, where switching to coupled solver doesn't (coupled solver uses 1.5x to 2x the memory because it stores the solution to the last iteration).As far as solver settings... it's very important to set your control limits. For temp and pressure, I calculate the delta between freestream and stagnation, and Idouble that. If I have a problem getting the model to converge, I may cut that down to about 1.5 to 1.2 times the difference between freestream and stagnation. (So my Pressure limits are Pstatic + 2*Q and Pstatic - 2*Q) I set the pressure and momentum URFs to 0.5 and 0.4, leaving energy at 1. The SIMPLEPressure-Velocity coupling tends to work well... I've gotten some recommendations on switching to PISO, so I'm actually trying that right now. I tend to run for 10 to 50 iterations with all the default discretizations and the turb eq turned off. Then I turn the turb eq back on (I don't go all the way to convergence like some people recommend... I haven't found any benefit... I usually only do about 50 iterations with these settings). After that, I set everything to 2nd order discretization and run to convergence. These settings have worked well for me on aerodynamic models with little or no separation, and they've worked (with a little playing around in the URFs) all the way to Mach 1.2.I have run bluff bodies at subsonic compressible to transonic speeds (typically a symmetrical model to get more of a "time averaged" solution... this avoids the oscillating vortices and cuts out the need for the unsteady solver... I haven't run transonic of a full bluff body, but I have run Mach 0.5ish with a full model and the unsteady solver... that was a while ago though, and I don't remember if I had made any changes to my solver settings). The coupled solver is a poor solver for this kind of model because of the large separation region aft of the body. This becomes a difficult model for the segregated solver as well, but I've had good luck running it with default discretization and the turb model turned off for about 100 to 200 iterations... then turning on the turb model, and lowering the energy URF to about 0.7 and running for another 100 or so iterations. Then upping the energy URF back to 1 for another 50 or so, and then switching to second order and running to convergence.The most important things I have found are paying attention to your mesh (or using adaption or dynamic adaption to resolve shocks... you don't need to refine them all the way for overall forces, but you might for flutter analysis) and setting your control limits.Hope this helps, and good luck,JasonThank you Jason, your explanations are complete and well detailed as usual. I run the Agard 445.6 at mach 1.141 and had no problem with convergence. I just set momentum URF to 0.4. Yes you're rigth when you say coupled solver requires more memory. In fact if I use the segregated I can make a more refined grid because it requires less memory. This is great! SIMPLE scheme seems to work well and I didn't need to switch to PISO...and I had no problem to solve the flow with transonic or supersonic flow. Thank you again for your answer. I just needed to have a confirm somebody else tried to use the segregated solver for not subsonic external flows. Luca1.我一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛~2.FLUENT的收敛最基础的是网格的质量,计算的时候看怎样选择CFL数,这个*经验3.首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。

相关文档
最新文档