不收敛怎么办？

雨烟子

不收敛应该调整那些地方，还有我的计算没有到设置的迭代次数，计算就停止了，是什么原因？

moitoi

溢出,检查网格边界条件

雨烟子

不太明白，网格边界条件，看趋势，计算值好像也平稳了，右边误差曲线也是减小的，但是没有到设置的迭代次数，计算就停止了。

bjs30

这种情况我也见过!!!我在塞扬上过不去,在P4上过了!
为什么?我一直没找到原因!

moitoi

各种机器允许的数值范围不同

bjs30

能不能自己控制?在运行程序时,自己确定!?

phoenics

有些问题收敛确实很难，但总有办法解决
1、松弛因子
2、网格
3、边界条件
4、给定合理的初始条件
5、改变计算当中的最大最小值限制
最主要还是上面3中方法

雨烟子

如果不收敛，值老是在一定范围内上下波动，怎么才能调整，松弛因子取大取小有什么影响。

易拉罐

一般说，不容易收敛的话应该把松弛因子改小一点为好

ching

边界条件怎么影响的？具体点好吗？

windymm

楼上二位大虾的疑问我也非常非常想知道！！！
我在试验中，松弛因子调得过小，反而不能收敛，我一头雾水。

phoenics

收敛是一个比较难回答的问题，只能靠经验。如果可能将你的Q1文件提供出来大家一起探讨一下。

windymm

我的q1文件：
TALK=T;RUN( 1, 1)

************************************************************
   Q1 created by VDI menu, Version 3.3, Date 03/05/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN   =       1 ;LIBREF =       0
************************************************************
  Group 1. Run Title
TEXT(No title has been set for this run.     )
************************************************************
  Group 2. Transience
STEADY  =    T
************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,48,38,19)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
XSI= 1.000000E+04; YSI= 1.100000E+04; ZSI= 2.800000E+03
RSET(D,CHAM    )
************************************************************
  Group 6. Body-Fitted coordinates
************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
NAME(150) =EPKE
    * Solved variables list
SOLVE(P1  ,U1  ,V1  ,W1  )
    * Stored variables list
STORE(EPKE)
    * Additional solver options
SOLUTN(P1  ,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)

************************************************************
  Group 8. Terms & Devices
************************************************************
  Group 9. Properties
SETPRPS(1,  0)
RHO1    = 1.189000E+00
PRESS0  = 1.000000E+05
TEMP0   = 2.730000E+02
CP1     = 1.005000E+03
ENUL    = 1.544000E-05
DVO1DT  = 3.410000E-03
PRT   (EP  ) =  1.314000E+00
************************************************************
  Group 10.Inter-Phase Transfer Processes
************************************************************
  Group 11.Initialise Var/Porosity Fields
FIINIT(EPKE) =  1.000000E+00
   No PATCHes used for this Group


INIADD  =    F
************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes used for this Group

EGWF    =    T
************************************************************
  Group 14. Downstream Pressure For PARAB
************************************************************
  Group 15. Terminate Sweeps
LSWEEP  =     800
RESFAC  = 1.000000E-03
************************************************************
  Group 16. Terminate Iterations
************************************************************
  Group 17. Relaxation
RELAX(P1  ,LINRLX, 1.000000E+00)
RELAX(U1  ,FALSDT, 1.000000E+02)
RELAX(V1  ,FALSDT, 1.000000E+02)
RELAX(W1  ,FALSDT, 1.000000E+02)
RELAX(KE  ,LINRLX, 5.000000E-01)
RELAX(EP  ,LINRLX, 5.000000E-01)
KELIN   =       3
************************************************************
  Group 18. Limits
VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
************************************************************
  Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
GENK    =    T
ASAP    =    T
************************************************************
  Group 20. Preliminary Printout
ECHO    =    T
************************************************************
  Group 21. Print-out of Variables
OUTPUT(P1  ,Y,Y,Y,N,Y,Y)
OUTPUT(U1  ,Y,Y,Y,N,Y,Y)
OUTPUT(V1  ,Y,Y,Y,N,Y,Y)
OUTPUT(W1  ,Y,Y,Y,N,Y,Y)
************************************************************
  Group 22. Monitor Print-Out
IXMON   =      18 ;IYMON  =       5 ;IZMON  =       8
NPRMON  =  100000
NPRMNT  =       1
TSTSWP  =      -1
************************************************************
  Group 23.Field Print-Out & Plot Control
NPRINT  =  100000
NXPRIN  =       1
IXPRF   =      14 ;IXPRL  =      20
NYPRIN  =       1
IYPRF   =      28 ;IYPRL  =      28
NZPRIN  =       1
IZPRF   =       1 ;IZPRL  =      13
ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
************************************************************
  Group 24. Dumps For Restarts
NOWIPE  =    T

GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        1.000000E+04, 1.100000E+04, 2.800000E+03
> DOM,    MONIT,       3.852000E+03, 1.934430E+03, 1.117910E+03
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     11,-1.300000E+00
> GRID,   RSET_X_2,      2, 1.000000E+00
> GRID,   RSET_X_3,      5, 1.000000E+00
> GRID,   RSET_X_4,      2, 1.000000E+00
> GRID,   RSET_X_5,     -5, 1.100000E+00
> GRID,   RSET_X_6,     -5, 1.200000E+00
> GRID,   RSET_X_7,      6,-1.200000E+00
> GRID,   RSET_X_8,     12, 1.300000E+00
> GRID,   RSET_Y_1,      9,-1.500000E+00
> GRID,   RSET_Y_2,      2, 1.000000E+00
> GRID,   RSET_Y_3,    -15, 1.400000E+00
> GRID,   RSET_Y_4,      2, 1.000000E+00
> GRID,   RSET_Y_5,     10, 1.500000E+00
> GRID,   RSET_Z_1,      6, 1.000000E+00
> GRID,   RSET_Z_2,     -7, 1.100000E+00
> GRID,   RSET_Z_3,      3, 1.000000E+00
> GRID,   RSET_Z_4,      3, 1.000000E+00
> DOM,    RELAX,       5.000000E-01

> OBJ1,   NAME,        B1
> OBJ1,   POSITION,    3.000000E+03, 3.000000E+03, 0.000000E+00
> OBJ1,   SIZE,        3.840000E+03, 4.740000E+03, 2.800000E+03
> OBJ1,   CLIPART,     q
> OBJ1,   ROTATION24,        0
> OBJ1,   TYPE,        BLOCKAGE
> OBJ1,   MATERIAL,      198

> OBJ2,   NAME,        B2
> OBJ2,   POSITION,    4.920000E+03, 3.000000E+03, 9.000000E+02
> OBJ2,   SIZE,        1.000000E+03, 2.400000E+02, 1.500000E+03
> OBJ2,   CLIPART,     cube14
> OBJ2,   ROTATION24,        0
> OBJ2,   TYPE,        BLOCKAGE
> OBJ2,   MATERIAL,      198

> OBJ3,   NAME,        B3
> OBJ3,   POSITION,    3.920000E+03, 3.000000E+03, 0.000000E+00
> OBJ3,   SIZE,        2.000000E+03, 2.400000E+02, 9.000000E+02
> OBJ3,   CLIPART,     c-low
> OBJ3,   ROTATION24,        0
> OBJ3,   TYPE,        BLOCKAGE
> OBJ3,   MATERIAL,      198

> OBJ4,   NAME,        B4
> OBJ4,   POSITION,    3.920000E+03, 3.000000E+03, 2.400000E+03
> OBJ4,   SIZE,        2.000000E+03, 2.400000E+02, 4.000000E+02
> OBJ4,   CLIPART,     c-up
> OBJ4,   ROTATION24,        0
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      198

> OBJ5,   NAME,        B5
> OBJ5,   POSITION,    3.240000E+03, 7.500000E+03, 2.000000E+03
> OBJ5,   SIZE,        9.000000E+02, 2.400000E+02, 8.000000E+02
> OBJ5,   CLIPART,     m-up
> OBJ5,   ROTATION24,        0
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      198

> OBJ6,   NAME,        IN
> OBJ6,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ6,   SIZE,        1.000000E+04, 0.000000E+00, 2.800000E+03
> OBJ6,   CLIPART,     cube3t
> OBJ6,   ROTATION24,        1
> OBJ6,   TYPE,        INLET
> OBJ6,   PRESSURE,      0.000000E+00
> OBJ6,   VELOCITY,      0.000000E+00, 5.000000E+00, 0.000000E+00
> OBJ6,   TEMPERATURE,   0.000000E+00
> OBJ6,   TURB-INTENS,   5.000000E+00

> OBJ7,   NAME,        OUTLET
> OBJ7,   POSITION,    0.000000E+00, 1.100000E+04, 0.000000E+00
> OBJ7,   SIZE,        1.000000E+04, 0.000000E+00, 2.800000E+03
> OBJ7,   CLIPART,     cube12t
> OBJ7,   ROTATION24,        1
> OBJ7,   TYPE,        OUTLET
> OBJ7,   PRESSURE,      0.000000E+00
> OBJ7,   TEMPERATURE,  -1.026000E+04
> OBJ7,   COEFFICIENT,   1.000000E+03
> OBJ7,   TURBULENCE,   -1.026000E+04,-1.026000E+04
STOP
收敛情况：

horse

各位高手说说怎么把Q1里面的设置改动后使之生效吧，我一直没弄成功，太失败了……

moitoi

这样改一下试试
RELAX(P1  ,LINRLX, 0。5)
RELAX(U1  ,FALSDT, 0。01)
RELAX(V1  ,FALSDT, 0。01)
RELAX(W1  ,FALSDT, 0。01)
RELAX(KE  ,LINRLX, 5.000000E-01)
RELAX(EP  ,LINRLX, 5.000000E-01)
KELIN   =       3