VOF中为何加了表面张力效应后结果不对呢？
本想计算一个微管道中的毛细管效应，但算出的结果很奇怪。望大侠们指点一下命令流如下
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 编辑 ] 怎么奇怪？说说你的边界条件和想法？
页:
[1]