3dec流热耦合问题请教

wangcbkxjh

请问这个热源命令有问题吗，为什么一打开set convection on热源温度就透不过裂缝
new
config  thermal gw
set th_numerical on
set atol 0.0008
poly reg 0 0.75,0 0.30,0 1.0
poly reg -0.0015 -0.0005,0 0.3,0.4 0.8
jset dip 90 dd 90 origin 0.284,0,1 id=1
jset dip 90 dd 90 origin 0.484,0,1 id=2
jset dip 15 dd 90 origin 0,0,0.4 id=3
jset dip 15 dd 90 origin 0,0,0.8 id=4
joint on
prop jmat=1 azero 2e-3 ares 2e-3 amax 2e-3
prop jmat=2 azero 1.5e-3 ares 1.5e-3 amax 1.5e-3
change jmat=1 range dip 90 dd 90
change jmat=2 range dip 15 dd 90
fluid density 1000 bulk 2e9  visc 1e-3
fluid fspec_heat 4180.0 fth_cond 0.6  fht_coe 10
gen edge 0.4
prop mat=1  dens 2700  k 6e9 g 4e9
prop mat=1 spec_heat 750  cond 2.26 diffus 3.2e-5
prop jmat=1 jkn 1e9 jks 1e9
prop jmat=2 jkn 1e9 jks 1e9
insitu pp 1e6
fix range x -0.0015 1.0
set flow off therm on mech off
cycle 100
bou disch 1.2e-3 range x 0 0.75 y 0 0.3 z 0.99 1.01
bou disch -1.2e-3 range x 0 0.75 y 0 0.3 z -0.01 0.01
bou disch 9e-4 range x -0.0005 0 y 0 0.3 z 0.79 0.81
bou disch 0 range x -0.0005 0 y 0 0.3 z 0.39 0.41
bou disch 9e-4 range x -0.01 0.01 y 0 0.3 z 0 1
bou disch -9e-4 range x 0.749 0.751 y 0 0.3 z 0 1
set flow on therm off mech off
cycle 100
bound flux=-5.24 range x -0.01 0.01 y 0 0.3 z 0 1
bound flux=-5.24 range x 0.749 0.751 y 0 0.3 z 0 1
bound flux=-5.24 range x 0 0.75 y -0.01 0.01 z 0 1
bound flux=-5.24 range x 0 0.75 y 0.299 0.301 z 0 1
bound flux=-5.24 range x 0 0.75 y 0 0.3 z -0.01 0.01
bound flux=-5.24 range x 0 0.75 y 0 0.3 z 0.999 1.001
bound fluidtemp 26 range  x 0 0.75 y 0 0.3 z 0.99 1.01
bound fluidtemp 26 range x -0.0005 y 0 0.3 z 0.79 0.81
bound fluidtemp 26 range x -0.01 0.01 y 0 0.3 z 0 1
fluid fht_coe 10
apply thermal convection 50 1e99 range x -0.0015 -0.0005
hist thtime
set therm on mech off flow off
set th_implicit on
set convection on
set thdt 1
step 10000
plot contour temp
plot fpcontour temp
plot fpcontour pp
return