求助：在多组分分析中如何选择某一流体（有实例）

crow0905

求助，在多组分分析中如何选中参与运算的某一流体？

比如在下面的例子中，如何在后处理阶段选择H2 或者O2 亦或N2来分别显示他们的温度？

/filename,spec4
/prep7
!!!!!!!  Define some dimensions - SI units are used
lenin=.3
half=lenin/2.
width=.1
hfwid=width/2.
outlen=1.2
!!!!!!!  Define inlet and outlet rectangles
rect,-lenin,0,-hfwid,hfwid
rect,.2,.3,.25,.25+lenin
rect,.2,.3,-.25,-.25-lenin
rect,.4,.4+outlen,-.15,.15
!!!!!!!  Lines that border the mixing area
l2tan,-3,-8
l2tan,-6,-15
l2tan,-13,-10
l2tan,-12,-1
!!!!!!!  Parameters defined for meshing
nlcurv=16
rlcurv=1
nscurv=9
rscurv=1
nispan=8
rispan=-1.5
nospan=13
rospan=-2
nilen=10
rilen=-2
nolen=24
rolen=4
!!!!!!!
!flst,2,2,4,orde,2
!fitem,2,17
!fitem,2,20
lsel,s,,,17,20,3
lesize,all,,,nlcurv,rlcurv
lsel,s,,,19
!fitem,2,18
!fitem,2,-19
lesize,19,,,nscurv,rscurv
lsel,s,,,18
lesize,18,,,nscurv,rscurv
!flst,2,3,4,orde,3
!fitem,2,2
!fitem,2,5
!fitem,2,11
lsel,s,,,2,5,3
lsel,a,,,11
lesize,all,,,nispan,1
!flst,2,3,4,orde,3
!fitem,2,4
!fitem,2,7
!fitem,2,9
lsel,s,,,4,7,3
lsel,a,,,9
lesize,all,,,nispan,rispan
!flst,2,6,4,orde,6
!fitem,2,1
!fitem,2,3
!fitem,2,6
!fitem,2,8
!fitem,2,10
!fitem,2,12
lsel,s,,,1,3,2
lsel,a,,,6,12,2
lesize,all,,,nilen,rilen
lsel,s,,,13,16
lesize,13,,,nolen,rolen
lesize,15,,,nolen,1./rolen
lesize,16,,,nospan,rospan
lesize,14,,,nospan,rospan
alls
!!!!!!!  Define the mixing area...
!flst,2,8,3
!fitem,2,2
!fitem,2,12
!fitem,2,11
!fitem,2,13
!fitem,2,16
!fitem,2,6
!fitem,2,5
!fitem,2,3
a,2,12,11,13,16,6,5,3
alls
/com ** CONSTRUCTION OF THE MESH
!!!!!!! Put triangles in the mixing area (5)
asel,s,,,5
et,1,141
mshape,1,2d
amesh,5
!!!!!!! Put a mapped mesh in the rectangles (quads)
asel,s,,,1,4
mshape,0,2d
mshkey,1
amesh,all
/com ** BOUNDARY CONDITIONS
!!!!!!! Wall boundary conditions
lsel,s,,,1,3,2
lsel,a,,,6,12,2
lsel,a,,,13,17,2
lsel,a,,,18,20
nsll,,1
d,all,vx,0
d,all,vy,0
!!!!!!! Define velocities, temperature at top, bottom, left
vtop=.1
vbot=.1
vlef=.1
lsel,s,,,7
nsll
d,all,vx,0
d,all,vy,-vtop
nsll,,1
d,all,temp,300
lsel,s,,,4
nsll
d,all,vx,vlef
d,all,vy,0
nsll,,1
d,all,temp,400
lsel,s,,,9
nsll
d,all,vx,0
d,all,vy,vbot
nsll,,1
d,all,temp,300
lsel,s,,,14
nsll,,1
d,all,pres,0
flda,conv,iter,5
save
!finish
alls
/solu
/com **  FLOTRAN INPUT
flda,iter,exec,40
flda,temp,nomi,300      ! Initial temperature field is 300K
flda,prot,dens,air-si   ! Initial solution will be for AIR
flda,vary,dens,true
flda,prot,visc,air-si
flda,vary,visc,true
flda,prot,cond,air-si
flda,vary,cond,true
flda,prot,spht,air-si
save
solve                  ! Solve 40 iterations with air at 300K
finish                 ! Prepare for species transport ...
/prep7
flda,prot,dens,cmix    ! Composite mixture for density
flda,prot,visc,cmix     ! Composite mixture for viscosity
flda,prot,cond,cmix
keyopt,1,1,3         ! Specify that there are 3 species
flda,solu,spec,t       ! Turn on solution for multiple species transport
msdata,2                ! Specify species 2 (N2) as the algebraic species
!!!!!!!  Species Property Input
msspec,1,o2,31.999
msprop,1,DENS,GAS,1.2998,300,1.01325E+5
msvary,1,dens,t
msnomf,1,.3
msprop,1,VISC,CONSTANT,1.2067E-5
msprop,1,mdif,CONSTANT,2.149E-5
msprop,1,cond,CONSTANT,.02674
mscap,1,1
msspec,2,n2,28.018
msprop,2,DENS,GAS,1.1381,300,1.01325E+5
msvary,2,dens,t
msnomf,2,.3
msprop,2,VISC,CONSTANT,1.786E-5
msprop,2,mdif,CONSTANT,1.601E-5
msprop,2,cond,CONSTANT,.02598
mscap,2,1
msspec,3,h2,2.016
msprop,3,DENS,GAS,0.0819,300,1.01325E+5
msvary,3,dens,t
msnomf,3,.4
msprop,3,VISC,CONSTANT,8.94E-6
msprop,3,mdif,CONSTANT,4.964E-5
msprop,3,cond,CONSTANT,.1815
mscap,3,1
msrelx,1,1.0
msrelx,3,1.0
msmeth,1,3
mssolu,1,,,2,1.e-8
msmeth,3,3
mssolu,3,,,2,1.e-8
alls
!!!!!!!  Species boundary conditions
lsel,s,,,7
nsll,,1
d,all,o2,1.0
d,all,n2,0.0
d,all,h2,0.0
lsel,s,,,4
nsll,,1
d,all,o2,0.0
d,all,n2,1.0
d,all,h2,0.0
lsel,s,,,9
nsll,,1
d,all,o2,0.0
d,all,n2,0.0
d,all,h2,1.0
alls
!!!!!!!  End of species boundary conditions
flda,iter,exec,20           ! Ask for 20 global iterations
flda,conv,outp,land         ! Adjust convergence monitor output style
save
finish
/solu
solve                       ! 20 iterations with species activated
!
!Prepare for energy solution
flda,solu,temp,t            ! Achieve a constant flow temperature solution
flda,solu,flow,f            ! Freeze the flow field
flda,meth,temp,3            ! Activate PCCR solver
flda,conv,temp,1.e-10       ! Convergence criterion for PCCR
flda,iter,exec,5            ! Need only a few iterations for solution...
flda,relx,temp,1.0          ! no relaxation on temperature
save
solve
!! Prepare for coupled solution
flda,iter,exec,35            ! Achieve a flow and thermal solution
flda,solu,flow,t
solve
flda,iter,exec,50           ! Run 50 more iterations to refine solution
save
solve

crow0905

上面的例子是ANSYS  HELP中的例子，主要就是三种气体H2，O2 ，N2从3个入口进入一个管中的情况。目前我只会在后处理中看到整个流体的速度，密度，温度等物理量，H2，O2 ，N2 的质量组分也可以看到，但是无法分别观测他们的温度......
希望有知道方法的高手指点指点，在下感激不尽

crow0905

怎么没人回答阿......顶起来