【讨论】这个问题的动力来源于哪里呢？

femsun

这是ANSYS的VOF中的第三个例题。
LX   = 30                                  !   x length in micrometer (um)
LY   = 30                                   !   y length in micrometer (um)
RadF = 20                               !   radius of the final sphere in micrometer (um)
aRad = 18                                !   x radius in micrometer (um)
bRad = RadF**3/aRad**2     !   y radius in micrometer (um)
DELX = 0.1*RadF                  !   x spacing in micrometer (um)
DELY = 0.1*RadF                   !   y spacing in micrometer (um)
DELT = 2.0e-4                        !   time increment in millisecond (ms)
APPT = DELT*5                      !   append time in millisecond (ms)
NDT  = 1000                            !   total number of time steps
Sigm = 0.073                          !   SFTS coefficient in unit of (mg)/(ms)^2
  
rho  = 1.0e-9                           !   density in unit of (mg)/(um)^3
mu   = 1.0e-6                            !   viscosity in unit of (mg)/(um)(ms)
pexact=2*Sigm/RadF            !   theoretical pressure in equilibrium
  
/prep7   
et,1,141,,,1                               !   RY geometry
rect,0,LX,-LY,LY
lsel,s,,,2,4,2
lesize,all,DELY
lsel,s,,,1,3,2
lesize,all,DELX
amesh,all
  
nsel,s,loc,x,0                                     ! symmetry condition
d,all,vx,0.0
alls
icvf,elpt,0.0,0.0,aRad,bRad                 ! initial VFRC
ic,all,vx,0                                              ! initial velocity
ic,all,vy,0
save
fini
  
/solu
flda,solu,tran,t
flda,solu,vof,t
flda,solu,sfts,t
flda,nomi,sfts,sigm
flda,nomi,dens,rho
flda,nomi,visc,mu
flda,time,glob,10
flda,time,numb,NDT
flda,time,step,DELT
  
flda,advm,mome,supg
flda,mir,mome,1.0
flda,relx,pres,1.0
flda,time,appe,APPT
solve   

femsun

这是某个子步的速度求解结果，在上边的命令流中，设置了走遍的边界条件d,all,vx,0.0  
       alls  
设置了流体初始存在区域和所有分析区域速度为0
     icvf,elpt,0.0,0.0,aRad,bRad ! initial VFRC  
     i c,all,vx,0 ! initial velocity  
     ic,all,vy,0  
但是动力来自哪里呢？