【讨论】一个多相问题

femsun

一个多 相问题，说明如下：是ANSYS帮助中一个例题，在FLUIT THERMAL分析中的 ，下部和上边各有一种流体，中间传热体，图形如下：

femsun

它的 命令流：
!!!!模型尺寸参数
LX    = 10                           
NDX   = 10                          
LY1   = 2                           
LY2   = 0.1                        
LY3   = 1.0                          
NDY1  = 10
NDY2  = 3   
NDY3  = 10      
YOFFSET = 0               
  
    
/prep7  
et,1,141                              ! 2-D XY system
  
keyopt,1,1,2                        !  设置两个物相
esha,2                                  !Quad elements
!!!!
rect,,LX,YOFFSET,YOFFSET+LY1                  
rect,,LX,YOFFSET+LY1,YOFFSET+LY1+LY2          
rect,,LX,YOFFSET+LY1+LY2,YOFFSET+LY1+LY2+LY3
nummrg,all
numcmp,all
  
！网格设置
lsel,s,,,3,9,3
lsel,a,,,1
lesi,all,,,NDX,      
lsel,s,,,2,4,2
lesi,all,,,NDY1,-5.0   
lsel,s,,,5,7,2
lesi,all,,,NDY2,-5.0   
lsel,s,,,8,10,2
lesi,all,,,NDY3,-5.0   
allsel
  
！材料属性分配
mat,1             !for fluids material must be 1
amesh,1           !mesh fluid 1 region
amesh,3           !mesh fluid 2 region
mat,2             !for solid set material to 2
amesh,2           !mesh solid region
  
! 对下边区域流体入口边界条件设置
lsel,s,,,4
nsll,s,1                     ！选择全部节点和 选择内部的 节点的 区别是什么
d,all,pres,10            !inlet pressure & temp specified for fluid 1
d,all,temp,100
  
!对下边区域流体出口边界条件设置
lsel,s,,,2
nsll,s,1
d,all,pres,0.       !outlet pressure specified for fluid 1
  
!对下边区域流体底边边界条件设置
lsel,s,,,1
nsll,s,1
d,all,vy,0        !symmetry surface for fluid 1 only
  
!对上边区域入口边界条件设置
lsel,s,,,8
nsll,s,1
d,all,pres,10.           !inlet pressure & temp specified for fluid 2
d,all,temp,400
  
!对上边区域出口边界条件设置
lsel,s,,,10
nsll,s,1
d,all,pres,0.     !outlet pressure specified for fluid 2
  
!对上边流体区域上边边界设置
lsel,s,,,9  
nsll,s,1
d,all,vx,0                       !wall boundary conditions  
d,all,vy,0   
allsel
  
!对中间隔离体属性设置，设置了密度传导系数和比热
mp,dens,2,5.                 
mp,kxx,2,100
mp,c,2,13
  
/SOLU
!对各个物相的初始条件进行了设置
nsel,s,loc,x,,lx       !select nodes for fluid 1
nsel,s,loc,y,o,ly1
ic,all,sp01,1.0         !set mass fraction for fluid 1
ic,all,sp02,0.0
nsel,s,loc,x,,lx
nsel,s,loc,y,ly1+ly2,ly1+ly2+ly3    !select nodes for fluid 2
ic,all,sp01,0.0                                !set mass fraction for fluid 2
ic,all,sp02,1.0
allsel
  
！定义流体类型
FLDA,ITER,EXEC,200              ！设置总体迭代步数
FLDA,PROT,DENS,CMIX         ! Fluid density
FLDA,PROT,VISC,CMIX           ! Fluid viscosity  
FLDA,PROT,COND,CMIX         ! Fluid conductivity
FLDA,PROT,SPHT,CMIX          ! Fluid specific heat
  
!设置所有的 流体属性为可变
flda,vary,dens,t                             !all properties MUST be variable
flda,vary,visc,t
flda,vary,cond,t
flda,vary,spht,t
  
！定义流体性质方程的 名义项参数
FLDA,NOMI,DENS,1.0          ! initial density for all fluid region
FLDA,NOMI,VISC,2.0            ! initial viscosity for all fluid region
FLDA,NOMI,COND,3.0          ! initial conductivity for all fluid region
FLDA,NOMI,SPHT,4.0           ! initial specific heat for all fluid region
！
！设置各个物相的属性
msprop,1,spht,constant,1       !sp heat for fluid 1
msprop,2,spht,constant,2       !sp heat for fluid 2
msprop,1,dens,constant,1      !density for fluid 1
msprop,2,dens,constant,2      !density for fluid 2
msprop,1,visc,constant,1.0    !viscosity for fluid 1
msprop,2,visc,constant,2.0     !viscosity for fluid 2
msprop,1,cond,constant,.1      !conductivity for fluid 1
msprop,2,cond,constant,.2      !conductivity for fluid 2
  
！收敛设置
FLDA,CONV,PRES,1.0E-10       ! PCCR convergence criterion
FLDA,TERM,PRES,1.E-09                  
  
FLDA,OUTP,SP01,T
FLDA,OUTP,SP02,T
save
SOLVE                       !solve for flow only

femsun

这个问题无法求解，系统报错。我是个新手，我 考虑了 一下，找了 一些问题，相信很多的新手都渴望和大家共讨论
  
问题1：flda,vary,既然设置了所有的流体的 属性为可变，为和在msprop
  
中又定义各个物相属性为constant？
  
问题2：esha,2  这个命令为设置单元形状，有物GUI菜单？
  
问题3：FLDA,CONV,PRES,1.0E-10       ! PCCR convergence criterion
  查了 帮助FLDATA6控制收敛监测量的输出,没有这个PRES项目，也 不
  
是后办解释的  PCCR收敛准则？
  
问题4：同样的 问题是FLDA,OUTP,SP01,T，FLDATA5是设置输出及保
  
存文件控制，它的第二个选项可以自己定义吗 ？迷惑！上边几个 问题也许就是部 能求解的 原因。
  
问题5：也 是我 最想知道的 一个 问题。关于材料分配。
！材料属性分配
mat,1             !for fluids material must be 1
amesh,1           !mesh fluid 1 region
amesh,3           !mesh fluid 2 region
mat,2             !for solid set material to 2
amesh,2           !mesh solid region
  难道在流体分析 中不要定义MAT吗？可以直接使用吗？MAT,1在前未见定义啊？
  **如果我 想在区域2中定义一种气体，自然存在，是静止的，不要计算时射入，区域1中后来射入，该如何 定义，在结构分析中可以定义两种单元两种属性分别划分网格，这里怎么 处理呢 ？

femsun

多相问题是个比较棘手的 问题，渴望各位斑竹和各位高手多多赐教，一这个题目为引子，抛砖引玉，给生手们讲解多相问题，既提高了大家的 水平也繁荣了论坛，谢谢，如果有说明不详细的 地方请指出来，谢谢！