VOF中为何加了表面张力效应后结果不对呢？

kstone

本想计算一个微管道中的毛细管效应，但算出的结果很奇怪。望大侠们指点一下
命令流如下
FINISH  
/CLEAR,NOSTART  
/FILNAME,SURFTEN,0  
!mg,ms,micrometer,K
/nopr
LS2    = 3                           ! S/2 length
LW2   = 1                            ! W/2 length
LD   = 0.5                         ! D Height        
LH0   = 0.2                          ! H0

time = 10e3                           ! Total simulation time
nts = 20                             ! Total number of time steps
delta = time/nts                      ! Time step size
appf = delta                          ! Results append frequency
VEL = -1e-6                            ! simulation time with velocity VEL
SURFTENS=29.7e-3                      !Surface tension
GRAVITY=9.8                           !Gravity acceleration                     
RHO   = 1.17e-9                      ! Fluid density
MU    = 3.0                          ! Fluid viscosity
WETANGLE=25                          !wetting angle

*dim,timutab,table,2,,,time          ! left wall displacement history
timutab(1,0)=0,time     
timutab(1,1)=0,VEL*time
/PREP7 $smrt,off
et,1,141                              ! 2-D XY system
keyopt,1,4,1                          ! enable disp DOFs
esha,2                                ! Quad elements
rect,,3,,5
!rect,,LS2,LH0,(LH0+LD)
lesize,all,0.1
amesh,1
EPLOT

lsel,s,,,4                           !
nsll,s,1
d,all,vx
d,all,ux
lsel,s,,,2                           !
nsll,s,1
d,all,vx
d,all,ux

lsel,s,,,1                            ! Stationary right wall                  
nsll,s,1
d,all,vx
d,all,vy      
d,all,ux,0.0
d,all,uy,0.0

alls

nsel,s,loc,y,0,2                  ! select lower half elements
esln,s,1
ice,all,vfrc,1.0                       ! apply vof initial conditions
allsel

/SOLU
FLDATA8,NOMI,SFTS,SURFTENS               !surface tension
FLDATA8,NOMI,WSCA,WETANGLE                !contact angle
FLDA,NOMI,DENS,RHO                    ! Nominal density
FLDA,NOMI,VISC,MU                     ! Nominal viscosity
!acel,,GRAVITY                        ! Apply gravity loading
flda,step,sumf,1
flda,solu,tran,t
flda,time,step,delta
flda,time,numb,nts
flda,time,tend,time
flda,time,glob,50
flda,time,appe,appf
flda,solu,ale,t                       ! Ale on
flda,solu,vof,t                       ! Vof on
flda,solu,sfts,t                      !Surface tension on
flda,advm,mome,supg
!flda,mir,mome,1.0
save
SOLVE
/DSCALE,1,1.0
/POST1  
PLVFRC,0
ANTIME,30,0.5, ,1,0,0,0

[ 本帖最后由 kstone 于 2006-10-31 16:53 编辑 ]

fangdian

怎么奇怪？说说你的边界条件和想法？