[讨论]亚松弛因子与残差的对应关系

nust

    我现在有个问题,就是非定常计算时,不是有亚松弛因子和courant number等参数要求设置嘛,现在我的问题是可压缩气流,湍流模型问K-ε模型，想问你courant number，turbulence Kinetic Energy,turbulence Dissipation Rate,turbulence Viscosity,solid的参数设置是如何控制残差continuity  X-Velocity,Y-Velocity,Z-Velocity,Energy,K,ε，其具体的对应关系是什么？我现在就是Z-Velocity的残差比较大，0.1的数量级，其他都可以到0.00001，
我这个问题是气体重大圆柱往小圆柱里跑,直径比为340/50,开始在x方向上有初始恒定压力,圆柱的轴向沿Z方向.
我从其他的模型得到初步的对应关系
1.   courant number=5,T.K.E=0.1,C.D.R=0.1,T.V=0.1,Solid=0.1时,
  continuity  X-Velocity  Y-Velocity    Z-Velocity    Energy     K         ε
            13          6           21            52           13      26       111

2.    courant number=5,T.K.E=0.2,C.D.R=0.2,T.V=0.2,Solid=0.2时,
  continuity  X-Velocity  Y-Velocity    Z-Velocity    Energy     K         ε
            11          6           18             47           11     17       67

3.   courant number=5,T.K.E=0.4,C.D.R=0.4,T.V=0.4,Solid=0.4时,
  continuity  X-Velocity  Y-Velocity    Z-Velocity    Energy     K         ε
            10          5            16             42          10     10       32

4.   courant number=10,T.K.E=0.4,C.D.R=0.4,T.V=0.4,Solid=0.4时,
  continuity  X-Velocity  Y-Velocity    Z-Velocity    Energy     K         ε
            5           3           9              30           5      9        25

5.   courant number=10,T.K.E=0.6,C.D.R=0.6,T.V=0.6,Solid=0.6时,
  continuity  X-Velocity  Y-Velocity    Z-Velocity    Energy     K         ε
           6            3          8               25          6       6        15
continuity,X-Velocity,Y-Velocity,Z-Velocity,Energy,K,ε下面的数字表示达到收敛时所需要的步数.

从上面的数据我得到,当courant number不变时,亚松弛因子的变化对 ε的残差起调整的作用.而亚松弛因子保持不变时,courant number的改变对continuity  X-Velocity  Y-Velocity , Energy有显著改变,对其他残差也有不同程度的调整.