continuity不收敛怎么解决？

hdsanfeng

是一个瞬态的模型，采用的压力入口（动态压力入口），压力出口定义的边界条件，计算前所有松弛因子的设置都是默认，计算的时候发现continuity的收敛很不理想，比其他几个参数高出2个数量级，收敛的速度开始很快，后来就很慢，再后来就几乎平了；如果要调松弛因子，应该调那个呢？

    还有就是y轴的开始残差线很平滑，算了20步就开始跳动（每步组大迭代次数设置的是30），我的压力是y轴进口的，这个情况算正常吗？因为我的压力是瞬态变化的，就是沿着y轴方向，那么y轴残差线后来出现跳动是不是正常呢？

xrs333

你说的是残差吧。它仅供参考，看看残差的定义你就知道了。

hdsanfeng

恩，这个懂，但是在计算非定常的时候，为了保证每个时间步最大迭代内收敛（系统默认20），每步最大迭代次数设置比较大（设置成80），这样，如果continuity不收敛，必须要等到每步最大迭代次数完成了（每个时间步都迭代80次）才进入到下次迭代，很浪费时间。所以想改善该项，让他很其他参考一起收敛，可以提早进入下次迭代。
     我的情况下，松弛因子到底是哪个？denstiny不知道是不是？我这个调小的时候，好想有点用。
     谢谢xrs333！

xrs333

你监测的可能是折算残差 (scaled residual)，如果迭代初场比较好的话（比如通过差值（interpolate）或FMG初始化得到的初场），它就不会太小，可能你永远也等不到10^-3。每步迭代次数那么多也不合理，或者增大criterion，或者减小时间步长，或者直接把最大迭代次数设为10或20。不发散的话不应减小松弛因子，减小它只会使收敛更慢。

hdsanfeng

我是监视的折算残差 (scaled residual)，但是这个有影响吗？有个normalize residual跟折算残差差别很大吗？你说初始场好的话就不会太小，其实我是期望他小，那是不是说初始场对continuity的收敛影响很大呢？还有，你说通过差值（interpolate）或FMG初始化得到的初场是怎么弄的，我很想试试看，但是找不到相关介绍。
    我现在的确是减少了时间步长（增大criterion我认为不好，毕竟-03已经算大了），手册也说最大迭代次数设为10-20，看来这个数值是合理的，我就默认20了。松弛因子我决定等前面设置改好了恢复默认再试试。
    再次谢谢xrs333！！！

hdsanfeng

对了，补充下！
我把算法改成Non-Iterative Time Advancement,这样我就不用设置最大迭代步数，这个算法有点奇怪，残差的确是下降到了-03，但是残差图好像倒过来了，我本来想把迭代了15步的传上来大家看看，但是我传不了。大家分享下自己的经验。
复制点数据，看看有什么意见，期待！
Updating solution at time levels N and N-1.
done.
  iter continuity x-velocity y-velocity z-velocity          k    epsilon     time/iter
    14 7.1218e-01 5.1248e-02 9.5850e-02 2.3053e-02 3.4099e-01 7.1485e-01  0:00:03    1
    14 9.3508e-04 2.5128e-05 2.3274e-05 1.5743e-05 2.7315e-03 1.5338e-03
Flow time = 0.00027999998, time step = 14
1 more time step

Updating solution at time levels N and N-1.
done.
  iter continuity x-velocity y-velocity z-velocity          k    epsilon     time/iter
    15 7.1668e-01 4.7957e-02 8.9173e-02 2.1571e-02 3.1394e-01 6.7095e-01  0:00:03    1
    15 1.0250e-03 2.6050e-05 2.4807e-05 1.6303e-05 2.1552e-03 1.0498e-03
Flow time = 0.00029999999, time step = 15

118228

一头雾水