24小时热门版块排行榜    

查看: 7429  |  回复: 28
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fangdian

木虫 (正式写手)

[求助] 一个老生常谈的话题----亚松弛因子

由于流体力学中要求解非线性的方程,在求解过程中,控制变量的变化是很必要的,这就通过松弛因子来实现的.它控制变量在每次迭代中的变化.也就是说,变量的新值为原值加上变化量乘以松弛因子。
     一般来讲,大家都是在收敛不好的时候,采用一个较小的欠松弛因子。 Fluent里面用的是欠松弛,主要防止两次迭代值相差太大引起发散。
      松弛因子的值在0~1之间,越小表示两次迭代值之间变化越小,也就越稳定,但收敛也就越慢。
   由于FLUENT所解方程组的非线性,我们有必要控制的变化。一般用亚松驰方法来实现控制,该方法在每一部迭代中减少了的变化量。亚松驰最简单的形式为:单元内变量等于原来的值加上亚松驰因子a与变化的积分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个相关的亚松驰因子。

注:在FLUENT中,所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题,但是对于一些特殊的非线性问题(如:某些湍流或者高Rayleigh数自然对流问题),在计算开始时要慎重减小亚松驰因子。
使用默认的亚松驰因子开始计算是很好的习惯。如果经过4到5步的迭代残差仍然增长,你就需要减小亚松驰因子。
    有时候,如果发现残差开始增加,你可以改变亚松驰因子重新计算。

  但是此时,我有一个问题,就是关于亚松弛因子-----liquid fraction update
  我的计算模型主要是3D的VOF和熔化凝固模型。对于亚松弛因子的设定如下:

本来都是默认设定的。liquid fraction update的默认值为0.9.
可是在计算中间,连续性残差不稳,出现突然上升的情况。


此时的Global Courant Number还是小于1的。一个解决办法是减小时间步。可是时间步减小,计算到指定时间所需要的时间步就多了。于是我将liquid fraction update的亚松弛因子由默认的0.9改为0.5接着往下计算,“躲过”了这次的不稳定。

  可是,这样直接修改松弛因子接着往下算,有没有问题?对结果有没有影响?
   有的朋友说既然改变了亚松弛因子,就要从头开始算。--------不是我不想从头开始算,主要是我前面的数据已经连续计算5天(120小时)了,重头再来的话,工程浩大。
    对于到底取多少松弛因子,有的朋友说多测试,只是问题又来了,我的一个数据计算得上十天,如果一个一个测试,那么时间。。。。。。?
   好了,我总结一下自己的问题吧。
  1 亚松弛因子是否可以在计算中途改变,改变前后结果是否有继承性?
  2 对于亚松弛因子的设定,除了按照默认设定,经验设定,还有其他办法吗?

  啰嗦了半天,谢谢大家了
回复此楼

» 收录本帖的淘帖专辑推荐

模拟相关

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jwm234

铜虫 (小有名气)

本人新手,总听说改松弛因子有可能会对计算有些帮助,我想顺带着问个问题。。。是不是不管我的松弛因子改成多少,那么最终收敛时得出的结果都是一样的,只不过最终收敛所需要的时间不一样?
9楼2012-11-06 22:07:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 29 个回答

47425804

木虫 (正式写手)


描述得真详细,不支持一下绝对过不去!
按照俺的理解,松弛因子中途改变完全没问题,我自己就经常采用这种临时调整松弛因子或调整时间步长的方法跳过一些发散点,不过后面又会再慢慢改回来。

» 本帖已获得的红花(最新10朵)

2楼2012-11-06 09:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

送鲜花一朵
引用回帖:
2楼: Originally posted by 47425804 at 2012-11-06 09:22:40
描述得真详细,不支持一下绝对过不去!
按照俺的理解,松弛因子中途改变完全没问题,我自己就经常采用这种临时调整松弛因子或调整时间步长的方法跳过一些发散点,不过后面又会再慢慢改回来。

谢谢你回帖,如果后面不手动改回来呢?影响大吗?
3楼2012-11-06 10:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sh073

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
fangdian: 金币+10, ★★★很有帮助, 谢谢,很有帮助。 2012-11-06 14:37:17
我也觉得应该没影响吧,改变松弛因子只是变了迭代变化量,过程还是那过程,有时候你要算非稳态还得先算个稳态的呢
4楼2012-11-06 10:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:47 by zNcEhTcH7ihq
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:50 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见