-
Continuity
收敛问题
FLUENT
求解器设置
FLUENT
求解器设置主要包括：
1
、压力
-
速度耦合方程格式选择
< br>2
、对流插值
3
、
梯度插值
4
、压力插值
下面对这几种设置做详细说明。
一、
压力
-
速度耦合方程求解算法
FLUENT
中主要有四种算法：
SIMPL
E
，
SIMPLEC
，
PISO
，
FSM
(1)SIMPLE(semi-implicit method for
pressure-linked equations)
半隐式连
接压力方程方法，是
FLUENT
的默认格式。
(2)SIMPLEC(SIMPLE-consistent)
。对于简单的问题收敛非常快速，不对压力进
行修正，所以压力松弛因
子可以设置为
1
(3)Pressure-Implicit
with Splitting of Operators (PISO)
。对非定常
流动
问题或者包含比平均网格倾斜度更高的网格适用
(4)Fractional
Step
Method
(FSM)
对非定常流
的分步方法。用于
NITA
格式，与
P
ISO
具有相同的特性。
二、对流插值（动量方程）
FLUE
NT
有五种方法：一阶迎风格式、幂率格式、二阶迎风格式、
M
USL
三阶格式、
QUICK
格式
（
1
）
FLUENT
默认采用一阶格式。容易收敛，但精度较差，主要用于初值
计算。
（
2
）
Power Lar.
幂率格式，当雷诺数低于
5
时，计算精度比一阶格式要高。
（
3
）
二阶迎风格式。
二阶迎风格式相对于一阶格式来说，
使用更小的截断误差，
适用于三角形、
四面体网格或流动与网格不在同一直线上；
二阶格式收敛可能比
较慢。
（
4
）
MUSL(monotone
upstream-centered schemes for conservation laws).<
/p>
当地
3
阶离散格式。主要用于非结构网格
，在预测二次流，漩涡，力等时更精确。
（
< br>5
）
QUICK
（
Quadratic upwind interpolation
）格式。
此格式用于四边形
/
六
面体时具有三阶
精度，用于杂交网格或三角形
/
四面体时只具有二阶精度。
三、梯度插值梯度插值主要是针对扩散项。
< br>FLUENT
有三种梯度插值方案：
green-
gauss cell-based
，
Green-gauss
node-based
，
least-
quares cell based.
（
1
）格林
-
高斯基于单元体。求解方法可能会出现伪扩散
。
（
2
）格
林
-
高斯基于节点。
求解更精确，最小化伪扩散，推荐用于三角形网格
上
（
3
）
基于单
元体的最小二乘法插值。
推荐用于多面体网格，
与基于节点的格
林
-
高斯格式具有相同的精度和格式。
四、压力插值压力基分离求解器主要有五种压力插值算法。
<
/p>
（
1
）标准格式（
Standard)
。为
FLUENT
缺省格式，对大表妹边界层附近的曲线
发现压力梯度流动求解精度会降低
（但不能用于流动中压力急剧变化的地方——
此时应该使用
< br>PRESTO!
格式代替）
（
2
）
PRESTO!
< br>主要用于高旋流，压力急剧变化流（如多孔介质、风扇模型等），
或剧烈弯曲的区
域。
（
3
）
Linear
（线性格式）。当其他选项导致收敛困难或出现非
物理解时使用此
格式。
（
4
）
second
or
der
（二阶格式）。用于可压缩流动，不能用于多孔介质、阶跃、
风扇、
VOF/MIXTURE
多相流。
< br>
（
5
）
Body Force Weighted
体积力。当体积力很大时，如高雷诺数自然
对流或
高回旋流动中采用此格式。
（
1
）连续性方程不收敛是怎么回事？
在计算过程中其它指数都收敛了，就
continuity
不收
敛
是怎么回事。
这和
fluent
程序的求解方法
SIMPLE
有关。
SIMPLE
根据连续方程推导出压
力修正方法求解压力。
由于连续方程
中流场耦合项被过渡简化，
使得压力修正方
程不能准确反映流场
的变化，从而导致该方程收敛缓慢。
你可以试验
SIMPLEC
方法，应该会收敛快些。<
/p>
在计算模拟中，
continuity
总
不收敛
，除了加密网格，还有别的办法吗？
别的条件都已经收敛了，
就差它自己了，
还有收敛的标准是什么？是不是到了一
定的尺
度就能收敛了，比如
10
－
e5
具体的数量级就收敛了
continuity
是质量残差，
具体是表示本次计算结果与上次计算结果的差别，
如果别的条件
收敛了，就差它。可以点
report
，打开里面
FLUX
选项，算出进口
与出口的质量流量差，看它
是否小于
0.5%.
如果小于，可以判断它收敛
.
(2) fluent
残差曲线图中
continuity
是什么含义？
是质量守恒方程的反映，
也就是连续性的残差。
这个收敛的快并不能说明你的计
算就一定正确，
还要看动量方程的迭代计算。
表示某次迭代与上一次迭代在所有
cells
积分的差值，
continuty
表示连续性方程的残差
（
3
）要加速
continuit
y
收敛该设置那些参数？
感觉需要调整
courant number
FLUENT
中
courant
number
是在耦合求解的时候才出现的。正确的调整，可
以更好地加速收敛和
解的增强稳定性。
courant number
实际上是指时间步长和空间步长的相对关系，系统自动减
小
courant
数，这种情况一般出现在存在尖锐外形的计算域，
当局部的流速过
大或者压差过大时出错，把局部的网格加密再试一下。
< br>
在
fluent
中，用
courant number
来调节计算的稳定性与收敛性。一般来
说，随着
couran
tnumber
的从小到大的变化，收敛速度逐渐加快，但是稳定性
< br>逐渐降低。所以具体的问题，在计算的过程中，最好是把
ourant
number
从小
开始设置，
看看迭
代残差的收敛情况，
如果收敛速度较慢而且比较稳定的话，
可<
/p>
以适当的增加
courant
number
的大小，根据自己具体的问题，找出一个比较合
适的
courant number
，
让收敛速度能够足够的快，而且能够保持它的稳定性。
个人认为这应该和你采用的算法有关
SIMPLE
算法是根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，
使得压力修正方程
不能准确反映流场的
变化，从而导致该方程收敛缓慢。试着用
S
IMPLEC
算法看看。
(4)
正在学习
Fluent,
模拟圆管内的流动，速度入口，出口
outflow
运行后
xy
的速度很快就到
1
e-06
了，但是
continuity
老是降不下去，维持在
1e-00
和
1e-0
3
之间，减小松弛因子好像也没什么变化大家有什么建
议吗？
你查看了流量是否平
衡吗？在
report->flux
里面操作，
mass flow rate
，
把所有进出口都选上
，
compute
一下，看看
nut
flux
是什么水平，如果它的值
小于总进口流量的
1%
，
并且其他检测量在继续迭代之后不会发生波
动，
也可以认
为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：
1.
网格质量，
主要可能是相邻单元的尺寸大小
相差较大，
它们的尺寸之比最好控
制在
1.2
以内，不能超过
1.4.
2.
离散格式及压力速度耦合方法，
如果是结构网格，建议使用高阶
格式，如
2
阶
迎风格式等
，
如果是非结构网格，
除
p
ressure
保持
standard
格式不变外，
其他
格式改用高阶格式；压力速度耦合关系，如果
使用
SIMPLE
，
SIMPLEC<
/p>
，
PISO
等
s
egerated solver
对联系方程收敛没有提高的话，可以尝试使用
coupled
solver
。另外，对于梯
度的计算，不论使用结构或非结构网格，都可以改用
node-
based
来提高计算精度。
一些情况：
1.
监测流场某个变量来判断收敛更合理一些
.
2.
网格质量
.
ty
inlet
boundary
conditions
are
not
appropriate
for
compressible
flow
problems.
（
4
）要加速
continuity
收敛该设置那些参数？
感觉需要调整
courant
number
FLUENT
中
courant
number
是在耦合求解的时候才出现的。正确的调整，可
以更好地加速收敛和
解的增强稳定性。
courant number
实际上是指时间步长和空间步长的相对关系，系统自动减
小
courant
数，这种情况一般出现在存在尖锐外形的计算域，
当局部的流速过
大或者压差过大时出错，把局部的网格加密再试一下。
< br>
在
fluent
中，用
courant number
来调节计算的稳定性与收敛性。一般来
说，随着
couran
tnumber
的从小到大的变化，收敛速度逐渐加快，但是稳定性
< br>逐渐降低。所以具体的问题，在计算的过程中，最好是把
ourant
number
从小
开始设置，
看看迭
代残差的收敛情况，
如果收敛速度较慢而且比较稳定的话，
可<
/p>
以适当的增加
courant
number
的大小，根据自己具体的问题，找出一个比较合
适的
courant number
，
让收敛速度能够足够的快，而且能够保持它的稳定性。
个人认为这应该和你采用的算法有关
SIMPLE
算法是根据连续方程推导出压力修正方法求解压力。
由于连续方程中流场耦合项被过渡简化，
使得压力修正方程
不能准确反映流场的
变化，从而导致该方程收敛缓慢。试着用
S
IMPLEC
算法看看。
ANSYS
仿真计算代做：模态分析，瞬态动力学，谐响应分析和谱
分析、
械结构的疲劳、
损伤，
CFD
流体；
结构的强度评估和优化；
企鹅：
690294845
利用
FL
UENT
不收敛通常怎么解决？
①、
一般首先是改变初值，尝试不同的初始化，事实上好像初始化很关键，对于收敛。
②、
FLUENT
的收敛最基础的是网格的
质量，计算的时候看怎样选择
CFL
数，这个靠经验
③、首先查找网格问题，如果问题复杂比如多相流问题，与模型、边界、
初始条件都有关系。
④、有时初始条件和边界条件严重影响收
敛性，曾经作过一个计算反反复复，通过修改网格，重新
定义初始条件，包括具体的选择
的模型，
还有老师经常用的方法就是看看哪个因素不收敛，然
后寻
找和它有关的条件，改变相应参数。就收敛了
⑤、
A.
检查是否哪里设定有误：比方用
mm
的
unit
建构的
mesh
，忘了
scale
；比方给定的边界条件不
合理。
B
从算至发散前几步，
看
presure
分布，
看不出来的话，
再算几步
,
看看问题大概出在那个区域。
C
网格，
配合第二点作修正，
就重建个
更漂亮的，
或是更粗略的来处理。
D
再
找不出来的话，
换个
solver
。<
/p>
⑥、解决的办法是设几个监测点，比如出流或参数变化较大的地
方，若这些地方的参数变化很小，
就可以认为是收敛了，尽管此时残值曲线还没有降下来
。
⑦、调节松弛因子也能影响收敛，不过代价是收敛速度。
亚松弛因子对收敛的影响
所谓亚松驰
就是将本层次计算结果与上一层次结果的差值作适当缩减，
以避免由于差值过大而引起<
/p>
非线性迭代过程的发散。用通用变量来写出时，为松驰因子（
Re
laxation
Factors
）。《数值传热学
-214
》
FLUE
NT
中的亚松驰：由于
FLUENT
所
解方程组的非线性，我们有必要控制变化。一般用亚松驰方法来
实现控制，该方法在每一
部迭代中减少了变化量。亚松驰最简单的形式为：单元内变量等于原来的值
加上亚松驰因
子
a
与变化的积：
分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方
程，包括耦合解算器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。
在
FLUENT
中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是
对于一些特殊的非线性问题（如：某些湍流或者高
Rayleigh
< br>数自然对流问题），在计算开始时要慎重
减小亚松驰因子。
使用默认的亚松驰因子开始计
算是很好的习惯。
如果经过
4
到
5
步的迭代残差仍然增长，
你就需要减
小亚松驰因子。
有时候，如果发现残差开始
增加，你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现
这种情况。
最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件，
并对解的算法做几
步迭代以调节到新的参数。最典型的情况是，亚松驰因子的增加会使
残差有少量的增加，但是随着解
的进行残差的增加又消失了。
如
果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好
的数据文件。
注意：粘性和密度的亚松驰是在每一次迭代之间的。而且，如果直接
解焓方程而不是温度方程（即：
对
PDF
计算），基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值，你可以在解
< br>控制面板点击默认按钮。
对于大多数流动，不需要修改默认亚松弛因子。但是，如果出现不稳定或者发散你就需要
减小默认
的亚松弛因子了，
其中压力、
动量、
k
和
e
的亚松弛因子默认值分别为
0.2
，
0
.5
，
0.5
和
0.5
。
对于
SIMPLEC
格式一般不需要减小压力的亚松弛因子。
在密度和温度强烈耦合的问题
中，
如相当高的
Rayleigh
数的
自然或混合对流流动，应该对温度和
/
或密度（所用的亚松弛因子小于
1.0
）进行亚松弛。相反，当
温
度和动量方程没有耦合或者耦合较弱时，流动密度是常数，温度的亚松弛因子可以设为
1.0
。
对于其它的标量方程，如漩涡，组分，
PDF
变量，对于某些问
题默认的亚更松弛可能过大，尤其是对
于初始计算。你可以将松弛因子设为
0.8
以使得收敛容易。
⑧看了流量是否平衡
在
report->flux
里面操作，
mass
flow rate
，把所有进出口都选上，
compute<
/p>
一下，看看
nut flux
是什么水平
，如果它的值小于总进口流量的
1%
，并且其他检测量在继续迭
代之后不会发生波动，也可
以认为你的解是收敛的。
造成连续方程高残差不收敛的原因主要有以下几点：
1.
网格质量，主要可能是相邻单元的尺寸大小相差较大，它们的尺寸之
比最好控制在
1.2
以内，不
能超过<
/p>
1.4.
2.
离散格式及压力速度耦合
方法，如果是结构网格，建议使用高阶格式，如
2
阶迎风格式等
，如果
是非结构网格，除
pressure
保持
standard
格式不变外，其他格式改用高阶格式
；压力速度耦合关系，
如果使用
SIMPLE
< br>，
SIMPLEC
，
PISO<
/p>
等
segerated solver
对
联系方程收敛没有提高的话，可以尝试使用
coupled
s
olver
。另外，对于梯度的计算，不论使用结构或非结构网格，都可以改用
node-
based
来提
高计算精度。
在
FLUENT
的中文帮助里，对收敛有比较
详细地描述，建议去看看。
收敛性可能会受到很多因素影响。
大量的计算单元，
过于保守的亚松驰因子和复杂的流动物理性质
常常是主要的原因。有时候很难确定你是否得到了收敛阶。
没有判断收敛性的普遍准则。
残差定义对于一类问题是有用的，
但是有时候对其它类型问题会造成
误导。因此，最好的方法就是
不仅用残差来判断收敛性而且还要监视诸如阻力、热传导系数等相关的
积分量。
在
FLUENT
的
帮助文件里给出了下面几种典型的情况。
1
、
如果你对流场的初始猜测很好，
初始的连续性残差会
很小从而导致连续性方程的标度残差很大。
在这种情况下，
检查
未标度的残差并与适当的标度如入口的质量流速相比较是很有用的。
什么是标度
残差？就是选作用来标准化的残差值，一般是取第五步吧，所以，一开是残差就很小，那么，后面
的
残差和它一比，值也很难收敛到很小数。
2
、对于某些方程，如湍流量，较差的初始猜测可能会造成较高的标度因子。在
这种情况下，标度
的残差最开始会很小，随后会呈非线性增长，最后减小。因此，最好是
从残差变化的行为来判断收敛
性而不仅仅是残差的本身值来判断收敛性。你应该确认在几
步迭代（比如说
50
步）之后残差继续减小
或者仍然保持较低值，才能得出收敛的结论。
另一个判
断收敛性的流行方法就是要求未标度的残差减小到三阶量级。
为了实现这一方法，
FLUENT
提供了残差标准化，有关残差标准化的信息请参阅分离求
解器残差定义和耦合求解器残差定义两节。
在这种方法中，要求标准化的未标度残差降到
10^-3
。但是这种要求在很多情况下可能是不合适的。
1
、如果
你提供了较好的初始猜测，残差可能不会降到三阶量级。比方说，在等温流动中，如果
温
度的初始猜测非常接近最终值，那么能量残差根本就不会降到三阶量级。
2
、如果控制方程中包括的非线性
源项在计算开始时是零，但是在计算过程中缓慢增加，残差是
不会降到三阶量级的。例如
，在封闭区域内部的自然对流问题，由于初始的均一温度猜测不会产生浮
力，所以初始的
动量残差可能非常接近零。在这种情况下，初始的接近零的残差就不适合作为残差的
较好
的标度。
< br>3
、如果所感兴趣的变量在所有的地方都接近零，残差不会降到三阶量级。例如，
在完全发展
的管流中，截面上的速度为零。如果这些速度初始化为零，那么初始的和最终
的残差都接近零，因此
也就不能期待降三阶量级。
在这种情况下，最好监视诸如阻力、总热传导系数等积分量来判断解的收敛。检查非标准化未标度
的残差来确定这个残差和适当的标度相比是不是很小也是很有用的。
相反，如果初始猜测很差，初始的残差过大以至于残差下降三阶量级也不能保证
收敛。这种情况对
于初始猜测很难的
k
和
e
方程尤其常见。在这里，检查你所感兴趣的所有积分量就很
有用了。如果解是
不收敛的，你可以减少收敛公差
FLUENT
运行过程中，出现残差曲线震荡原因？如何解决残差震荡的问
题？残差震荡对计算收敛性
和计算结果有什么影响？
一
.
残差波动的主要原因：
1
、高精度格式；
2
、网格太粗；
3
、网格质量差；
4
、流场本身边界复
杂，流动复杂；
5
、模型的不恰当使用。
二
.
问：在进行稳态计算时候，开始
残差线是一直下降的，可是到后来各种残差线都显示为波形
波动，是不是不收敛阿？
答：有些复杂或流动环境恶劣情形下确实很难收敛。计算的精度（
2 <
/p>
阶），网格太疏，网格质
量太差，等都会使残差波动。经常遇到，
一开始下降，然后出现波动，可以降低松弛系数，我的问题
就能收敛，但如果网格质量不
好，是很难的。通常，计算非结构网格，如果问题比较复杂，会出现这
种情况，建议作网
格时多下些功夫。理论上说，残差的震荡是数值迭代在计算域内传递遭遇障碍物反
射
形成周期震荡导致的结果，与网格亚尺度雷诺数有关。例如，通常
压力边界是主要的反射源，换
成
OUTFLOW
边界会好些。这主要根据经验判断。所以我说网格和边界条件是主要因素。
三
. 1
、网格问题：比如流场内部存
在尖点等突变，导致网格在局部质量存在问题，影响收敛。
2
、可以调整一下
courant
number
，
courant
number
实际上是指时间步长和空间步长的相对关
系，系统自动减小
courant
数，这种情况一般
出现在存在尖锐外形的计算域，当局部的流速过大或者
压差过大时出错，把局部的网格加
密再试一下。
在
flu
ent
中，
用
courant
number
来调节计算的稳定性与收敛性。
一般来说，
随着
courant
nu
mber
的从小到大的变化，收敛速度逐渐加快，但是稳定性逐渐降低。所以具体的问题
，在计算的过程中，
最好是把
courant
number
从小开始设置，看看迭代残差的收敛情况，如果收敛速度较慢而且
比较稳定
的话，可以适当的增加
courant number
的大小，根据自己具体的问题，找出一个比较合适的
coura
nt
number
，让收敛速度能够足够的快，而且能够保持
它的稳定性。
turbulent viscosity
limited to viscosity ratio of
1.000000e+00
5 in 2 cells
1 <
/p>
这个应该是湍流模型的选取与第一层网格高度之间不满足近壁处理
关系而产生的问题，
如果你没有使用壁面函数的话，第一层网格高度尽
< br>可能地小点儿，比如取为弦长的百万分之一左右；另外，边界条件中关
于湍流量的
设置不合理也会导致这个警告。
2
（不推荐）
solve-controls-limits
Maximum Turb. Viscosity Ratio
加多两个
0
，估计
是一些单元的最大
Turb. Viscosity Ratio
超出了限定值
（）
恕我直言，
你的这个方法只是治标
不治本，
他这个问题多数是由于网格尺度太大引起的。
也可能是
边界条件上的湍流相关参数不合理导致的。
[br][br][
以下内容由
larky
在
200
7
年
06
月
2
3
日
00:00am
时添加
] [br]
调大限制值可能导致发散
3
这是一个办法，
能够解决一部分问题，
有一些
问题无论你怎么调整都没有用，
如果出现这
种情况可以通过调整
初始流场，
找到合适的初始值大部分能够解决，
其实如果只是一
开始初
现这个问题，可以不作调整，除非影响到收敛性
4
在别的论坛上看到的：
为了尽快收敛对异常的数值进行的限制，对最后收敛结果无影响
1
）如果边界条件设置合理，一般来说会在收敛后自动消除。<
/p>
2
）为了加快收敛对异常的数值进行的
限制（以引用
2
楼），是加快收敛的一种措施。
3
）但是如果你的问题中流场变化很大，有可能在最
有还会有。
4
）如果网格不好会经常
出现这种现象。
5
）如果不想看见它
总是报告而影响计算速度（写屏会降低计算速度），可以在下面把它关
闭：
solve->control->flow
limit....
具体记不住了，自己看看就知道了。
5
我也遇到这种情况
,
不过是在叠代求解的前一百多步
,
后面就没有了
.
因此我想是否是因为
前面计算的误差大引
起
?
而随着计算误差的减少
,
就消失了
.
如果是这样
,
就可以放心啦
.
6
一般是边界上或是网格质量差的地方出现了奇点
.
由
于是数值耗散
,
随着迭代次数越多
,<
/p>
影
响整个流场的范围越大
,
最终可能导致这个流场发散
.
如果是网格质量差的
地方出现
,
就只能重划网格了
如果是在边界上
,
一般是湍流相关参数设置不
合理造成的
,
改成固定湍流比可能能解决
对于
VOF
学习的几点体会
1
、对于网格的选择
尽量选择四边形或六面体网格，而不用三角形或四面体网格，以简化多相流问题。
2
、对于
VOF
< br>公式的选择
Geo-Reconstruct
：是目前最精确的界面跟踪方法，是对大多数瞬态
VOF
计算所推荐使用
的方法。
Donor-Acceptor
：在网格中包含很多扭曲的六面
体单元，可用该公式代替
Geo-
Reconstr
uct
。
Euler-
Explicit
：可解决
Donor-Acceptor
公式仅对四边形、六面体网格有效的问题，可用
于三角形或四面体网
格计算；也可以用于
Geo-
Reconstruct
不能给出满意结果的情形。
注意：
Geo-
Reconstruct
和
Donor-
Acceptor
，必须保证在区域内没有双边壁面。
3
、对于
VOF
模型
相的定义
原则上可以根据个人的喜好定义。
< br>如果其中有一相是可压缩的，
为了提高解的稳定性，
应指
定该相为主相。
注意：在
VOF
模型中，只能有一相是可压缩的。
4
、关于表面张力和壁面粘附
对于网格的选择：在表面张力有重大影响的计算区域内应使用四边形或六边形网格。
在打开
Wall
Ad
hesion
后，在
wall
边界条件
下指定接触角为每一对
相。（接触角定义：当系统达到平衡时
，在气、液、固三相交界处，气
-
液界面和固
< br>-
液界面
之间的夹角。实际反映的是液体表面张力与液<
/p>
-
固界面张力间的夹角。）
提高解的收敛性：对于涉及到表面张力的计算，建议在
Multiphas
e
Model
panel
中为
B
ody
Force
Formulation
打开
Implicit
Body
Force
。
这样做由于压力梯度和动量方程中
表面张力的部分平衡，从而提高解的收敛性
。
5
、关于运算环境的设置
对于
VOF
计算，应当选择
Specified
Operating
Density
，并且在
Operating
Density
下为最轻相设置密度。这样做排除了水力静压
的积累，提高了
round-
off
精度为动量平衡。
同样需要打开
Implicit
Body
For
ce
，部分平衡压力梯度和动量方程中体积力，提高解的
收敛性
。
Reference
Pressure
Location
（参考压强位置）应是位于流体永远是
100%
的某一相（空<
/p>
气）的区域，光滑和快速收敛是其基本条件。
6
、关于求解器的设置
压力插值方案：对于所有的
VOF
计算，应当使用<
/p>
body-force-
weighter
或者
PRESTO
!
压
力插值方案。
对于
Implicit
和
Euler
Explicit
方程，为了提高相间界面的清晰度，应当为
体积分率方程
选择采用
second-order
或者
QUICK
离散方案
速度压力耦合方案：通常瞬变流计算建议采用
PISO
方案。使用
PISO
时允许增加所有欠松
弛因子的值，而不会减弱解的稳定。对于
四
边形和三角形网格
上的计算，用
PISO
方案时，为了提高稳定性，建议为压力选欠松弛因
子为
0.7
-0.8
。
注意：
< br>当
FLUENT
进行任何模拟时，
如果前松弛因子设置为
1
时，
解出现
不稳定、
发散行为，
欠松弛因子必须减小；提高稳定性的另一个
方法是减小时间步长。（使用
稳态隐式
的
VOF
方案，为了提高稳定性，所有变量的欠松弛因子应设置在
0.2-0.5
之间。）
众所周知，
FLUENT
和
< br>CFX
是
ANSYS
中最牛的两
个流体分析软件。下面以
FLUENT
为例，
< br>说明其在多相流分析中的应用。
该例子来自于
FLUEN
T
帮助
,
但是其建模，
网格划分以及命名
集的定义方式则进行了改变。希望该例子对于大家做多相流的
分析有所帮助。
问题：一个水
-
空气混合物在管道内向上流动，在
T
型交叉
点分成两支。管道宽
25mm
，
输入部
分长
125mm
，顶部和右边都是
25
0mm
。空气和水在进口处的速度见下图，而两个出口
处的出流
权重分为为
0.38
和
0.62.
现在要求对该两相流做一个稳态分析。
使用
ANSYS
fluent
分析过程如下
（
1
）创建项目示意图（
WORKBENCH<
/p>
）
设置
geometry
单元格的属性
-
-
-
-
-
-
-
-
-
上一篇：数学专业名词
下一篇：fluent经验之谈(过来人的总结)word版本