热电场耦合问题 求解答

yutingting1986

我做的是两电极间空气的加热过程，但是两电极间施加电压后，空气相当于电阻，求解结束后中间空气的温度并没有变化，不知是什么原因。我的几何模型外边界为阳极，内边界为阴极边界，最上边最大的圆柱为被加热后流出的空气。
下面为我的命令流：
/batch,list
/show
/filename,65
/title,thermal-electric
/prep7
/units,si
et,1,solid69
mptemp,1,2000,2500,3000,3500,4000,4500
mptemp,7,5000,5500,6000,6500,7000,7500
mptemp,13,8000,8500,9000,9500,10000,10500
mptemp,19,11000,11500
mpdata,kxx,1,1,0.121, 0.252,0.383,0.44,0.5,0.69
mpdata,kxx,1,7,0.881,1.4,2.13,2.83,3.52,2.67
mpdata,kxx,1,13,1.97,1.56,1.12,1.12,1.12,1.26
mpdata,kxx,1,19,1.41,1.56
mpdata,rsvx,1,1,1000,0.56,4.67e-02,3.18e-03,4.11e-04,1.2e-04
mpdata,rsvx,1,7,3.4e-05,2.1e-06,8.85e-06,4.12e-06,2.56e-06,1.13e-06
mpdata,rsvx,1,13,9.61e-07,7.81e-07,5.26e-07,4.85e-07,3.66e-07,3.34e-07
mpdata,rsvx,1,19,2.92e-07,2.65e-07
cyl4,,,,,0.0075,360,0.0042
cone,0.0075,0.0025,0.0042,0.0073,0,360
local,11,1,,,.0073
wpave
cyl4,,,,,.0025,360,0.003
csys,0
wpave
cone,0.0025,0.005,0.0103,0.0128,,360
local,12,1,,,0.0128
wpave
cyl4,,,,,0.015,360,0.0372
csys,0
wpave
cyl4,,,,,0.003,360,0.002
cone,0.003,0.0005,0.002,0.0073,,360
csys,11
wpave
cyl4,,,,,0.0015,360,0.001




用GUI手动将体1,2,3减去6,7,8
cyl4,,,,,0.0015,360,0.001
csys,0
wpave
csdele,11
csdele,12
Vglue,all
numcmp,volu
numcmp,area
vsel,s,volu,,1,6
vatt,1,1,1

mshape,1,3d
esize,0.0006
vmesh,1
vmesh,3,6
esize,0.001
vmesh,2
/solu
antype,static
TOFFST,273
nsubst,20,100,2
neqit,40
*ABBR,NSL,NSLA,,1
asel,s,,,24
nsl
d,all,temp,893
asel,s,,,27,28
nsl
d,all,temp,1573
asel,s,,,19,20
nsl
d,all,temp,1873
asel,s,,, 16,17
nsl
d,all,temp,2073
asel,s,,,7,8
nsl
d,all,temp,1873
asel,s,,,13
nsl
d,all,temp,873
asel,s,,,9,11
nsl
d,all,temp,693
asel,s,,,21,22
nsl
d,all,temp,2473
asel,s,,,25,26
nsl
d,all,temp,2073
asel,s,,,14,15
nsl
d,all,temp,1773
asel,s,,,4
nsl
d,all,temp,3373

allsel


asel,s,,,7,8
nsl
d,all,volt,0
asel,s,,,16,17
nsl
d,all,volt,0
asel,s,,,19,20
nsl
d,all,volt,0
allsel


asel,s,,,9,11
bfa,all,ef,,0
asel,s,,,13,15
bfa,all,ef,0
asel,s,,,24
bfa,all,ef,,,0
asel,s,,,27,28
bfa,all,ef,0
asel,s,,,4
nsl
cp,1,volt,all
ni=ndnext(0)
f,ni,amps,6

allsel

sbctran
/solu
allsel
solve