关于ANSYS的FLOTRAN CFD 稳态热分析中的网格问题

liudong

在ANSYS的FLOTRAN CFD 稳态热分析中，有一个例题4.7. Examples of a Laminar, Thermal, Steady-State FLOTRAN Analysis，大家可能都作过，但我发现当网格数加大结果就变化很大，比如例题中是２Ｄ分析，网格数是２５*２５的计算结果和４０*４０的差别很大，请高手指教，很急，先谢了！以下是此题的命令流：
/TITLE,Buoyancy driven flow in a square cavity
/PREP7
ET,1,FLUID141                    ! Set 2-D fluid element type
RECTNG,,0.03,,0.03               ! Create square area
LPLOT
LESIZE,ALL, , ,25,1,0            ! Define line divisions and mesh area
AMESH,1
EPLOT
NSEL,S,EXT                       ! Apply velocity b.c. on exterior nodes
D,ALL,,,,,,VX,VY
NSEL,S,NODE,,1                   ! Apply thermal b.c. on left edge
NSEL,A,NODE,,52
NSEL,A,NODE,,77,100
D,ALL,TEMP,320
NSEL,S,NODE,,2                   ! Apply thermal b.c. on right edge
NSEL,A,NODE,,27,51
D,ALL,TEMP,280
NSEL,ALL
ESEL,ALL
FLDATA1,SOLU,TEMP,1              ! Set solution option to thermal
FLDATA2,ITER,EXEC,200            ! Set FLOTRAN execution controls
FLDATA2,ITER,OVER,50
FLDATA5,OUTP,SUMF,50
!*
FLDATA13,VARY,DENS,1             ! Set fluid properties
FLDATA7,PROT,DENS,AIR-SI
FLDATA7,PROT,VISC,AIR-SI
FLDATA7,PROT,COND,AIR-SI
FLDATA8,NOMI,COND,-1
FLDATA7,PROT,SPHT,AIR-SI
!*
ACEL,0,9.81,0                    ! Set flow environment gravity
FLDATA18,METH,PRES,1             ! Set FLOTRAN solver controls
FLDATA19,TDMA,PRES,100
FINISH
/SOLU
SOLVE                            ! Solve model
FINISH
/POST1
SET,LAST
PLNSOL,TEMP                      ! Plot temperature solution
PLNSOL,STRM                      ! Plot streamline contours
/DEVICE,VECTOR,1
!*
PLVECT,V, , , ,VECT,ELEM         ! Plot velocity vectors
wpstyle,0.0005,0.0001,-1,1,0.00005,0,2,,5
EPLOT
TRPOIN,P                         ! The P argument enables graphical picking
                                 ! of trace points via the GUI
PLTRAC,FLUID,TEMP                ! Plot particle traces contoured by
PLTRAC,FLUID,V,SUM               ! temperature and velocity magnitude
FINISH