方腔中稳态层流自然对流的热传导和对流的耦合问题讨论

vibration

方腔中稳态层流自然对流的热传导和对流的耦合问题讨论:
本问题其实为冰柜中温度场的计算雏形
此方腔尺寸很小，0.036258m×0.036258m，为二维稳态层流，
其中左右两边为plate，温度都为－70度
上下表面分别为外界环境30度，上下面都有固体隔热材料，
其中上表面温度设定为plate，设定为30度
而下表面温度用celltype，设定也为30度
计算后，收敛很好，
但是我有以下问题：
1，为什么用celltype设定的表面温度可以再计算后基本保持不变（29度左右），而用plate设定的表面温度（同为30度）就相差很远（25度左右）？其中区别在哪里？
2，此方腔尺寸为0.036258m×0.036258m，我发现方腔越大越难收敛，比如0.4m×0.4m后，计算就非常难收敛，我的网格也加密了，但是很难收敛，请问高手有什么办法吗？

vibration

TALK=T;RUN( 1, 1)

************************************************************
Q1 created by VDI menu, Version 3.5, Date 20/10/03
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This case is the three-dimensional counterpart of case 251,
but with a somewhat coarser grid in the y-z plane. The cavity
is now a cubical box, with adiabatic end walls. Only one half
of the domain is included, by reason of symmetry.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 278
************************************************************
Group 1. Run Title
TEXT(Laminar Natural Convec In A Cavity )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,30,30,1.000000E-05)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(150) =TEM1
SOLVE(P1 ,V1 ,W1 ,TEM1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)

************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+04 ;TEMP0 = 2.730000E+02
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -2.580000E-02
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1 ) = 1.000000E-06 ;FIINIT(V1 ) = 1.000000E-06
FIINIT(W1 ) = 1.000000E-06 ;FIINIT(TEM1) = 1.000000E+01
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

PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,V1 , FIXFLU , GRND3 )
COVAL (BUOYANCY,W1 , FIXFLU , GRND3 )

PATCH (DOWM ,CELL ,3,0,0,0,0,0,1,1)
COVAL (DOWM ,TEM1, 1.000000E+04, 3.000000E+01)

BUOYA = 0.000000E+00 ; BUOYB = 0.000000E+00
BUOYC =-9.810000E+00
BUOYD = 2.874000E-01
BUOYE = 2.000000E+01
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2000
SARAH = 5.000000E-03
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 30 ;LITER (V1 ) = 20
LITER (W1 ) = 20 ;LITER (TEM1) = 30
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
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
ASAP = T
PARSOL = F
CSG3 ='CNGR'
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 9 ;IZMON = 25
NPRMON = 100000
NPRMNT = 1
TSTSWP = -5
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NYPRIN = 4
NZPRIN = 6
NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000

PATCH (PROF ,PROFIL,1,1,1,20,15,15,1,1)
PLOT(PROF ,W1 , 0.000000E+00, 0.000000E+00)

PATCH (CONT ,CONTUR,1,1,1,20,1,30,1,1)
************************************************************
Group 24. Dumps For Restarts

GVIEW(P,9.945982E-01,-9.355936E-03,1.033781E-01)
GVIEW(UP,-1.034247E-01,-4.642130E-03,9.946265E-01)

> DOM, SIZE, 1.813000E-02, 3.625800E-02, 3.626000E-02
> DOM, MONIT, 9.065000E-03, 1.027310E-02, 2.994020E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_Y_1, 30, 1.000000E+00
> GRID, RSET_Z_1, 3, 1.000000E+00

> OBJ, NAME, HOT
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.813000E-02, 0.000000E+00, 3.626000E-02
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00,-7.000000E+01

> OBJ, NAME, COLD
> OBJ, POSITION, 0.000000E+00, 3.625800E-02, 0.000000E+00
> OBJ, SIZE, 1.813000E-02, 0.000000E+00, 3.626000E-02
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00,-7.000000E+01

> OBJ, NAME, LOWAL
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.813000E-02, 3.625800E-02, 0.000000E+00
> OBJ, CLIPART, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, CELLTYPE

> OBJ, NAME, FIXPRESS
> OBJ, POSITION, 0.000000E+00, 1.571180E-02, 1.708908E-02
> OBJ, SIZE, 1.813000E-02, 1.208599E-03, 1.223917E-03
> OBJ, CLIPART, cubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, PRESSURE_RELIEF
> OBJ, PRES_RELIEF, 1.000000E+03, 0.000000E+00

> OBJ, NAME, DOWNOBJ
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.813000E-02, 3.625800E-02, 3.625996E-03
> OBJ, CLIPART, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 134

> OBJ, NAME, CELLUP
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.626000E-02
> OBJ, SIZE, 1.813000E-02, 3.625800E-02, 0.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 3.000000E+01

> OBJ, NAME, UPOBJ
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.300000E-02
> OBJ, SIZE, 1.813000E-02, 3.625800E-02, 3.260000E-03
> OBJ, CLIPART, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 134
STOP