【讨论】流场求解中温度出现负值的问题(急)

cmjssl

我求解的三维流场问题,温度出现负值,提示说要么重新细剖,要么放松参数,可是我都试了,还是不行,不知为什么?而二维的由于电流升热在流场分析中求解温度,和实际情况不符,而且查看后处理时,出现下面的警告:
The current solution may have been produced using diffrent model or boundary condtion data than is currently stored. Post1 results may be erroneous unless you perform a new solution using the stored data.
不知为什么?  
有高手指点一二,同时希望斑竹关注一下,急呀

jnjl

产生错误的原因是多方面的，把你的题目大体说一下，或者介绍一下你的大概思路，这样大家也就能更好的替你出谋划策了。

cmjssl

非常感谢您的回应.
我的题目是在一个钢外壳中放几根导体,将电流产生的焦耳热读到流场中,属于空气自然对流,外壳外表面加自然对流换热系数,外壳内表面和导体外表面是固壁,对空气施加重力加速度.
我今天又发现用8.0版本就不出现那个警告,用的求解器不同得到的结果也不同,我这题目属于流固耦合,不知用哪个求解器是正确的.
对于三维的模型温度出现负值我还是不知怎么解决?我用的方法与二维的一样,二维的就不出现负值
请各位高手关注一下,多谢了

jnjl

关于用何种求解器的问题，看书中的介绍一般都是用默认的求解器，你看一下帮助中关于不同求解器的介绍，看看哪种更适合你的问题的求解。“不同的求解器得到的结果不同”，是结果差别很大吗？如果那样的话就有问题了。
二维和三维的计算差别很大的，流体分析中需要设置的参数太多，稍有不慎可能就会导致结果的不正确，你再仔细核对一下。

cmjssl

谢谢您的回复.我仔细核对过了,我自己找不出什么毛病,我把命令流贴上去,麻烦高手帮我看看.
！加交流电，求解磁场
/prep7
et,1,plane53,1
MP,MURX,1,1          !空气的相对磁导率
MP,MURX,2,1           !母排的相对磁导率 （铜）        
MP,MURX,3,500         !外壳的相对磁导率（钢）
MP,RSVX,1,1.7e8
MP,RSVX,2,1.7E-8      ！三者的电阻率
MP,RSVX,3,9.7E-8
！建模
RECTNG,0,0.18,0,0.18
RECTNG,0.002,0.178,0.002,0.178
RECTNG,0.086,0.094,0.025,0.155  
lfillt,9,12,0.004
lfillt,9,10,0.004
lfillt,10,11,0.004
lfillt,12,11,0.004
lsel,s,,,13,16,3
lsel,a,,,18,21,3
asbl,3,all
adele,4,7,1,1
numcmp,all
agen,3,3,,,0.02
agen,2,3,,,-0.02
AOVLAP,ALL
NUMCMP,ALL
/PNUM,AREA,1
/REPLOT
！赋属性
ASEL,S,,,6
AATT,1,,1
ASEL,S,,,5
AATT,3,,1
ASEL,S,,,1,4,1
AATT,2,,1
！分网
SMRT,3   
allsel,all
MSHAPE,0,2D  
MSHKEY,0
AMESH,all
  
/solu
antype,3
harfrq,0,50
kbc,0
LSEL,S,,,1,4,1
NSLL,S,1
D,ALL,AZ,0,0    ！加平行边界条件
allsel,all
ASEL,S,,,1
nsla,s,1
cp,1,volt,all
*get,n1,node,,num,min
f,n1,amps,2600,0
asel,s,,,2
nsla,s,1
cp,2,volt,all
*get,n2,node,,num,min
f,n2,amps,1300,-2251.6
asel,s,,,4
nsla,s,1
cp,3,volt,all
*get,n3,node,,num,min
f,n3,amps,-1300,-2251.6
ALLSEL,ALL
solve
finish
/post1
set,last
PLNSOL,B,SUM
FINISH
！求解流场
/PREP7
et,1,FLUID141
mp,kxx,2,399
mp,kxx,3,36.7
mp,c,2,385
mp,c,3,460
mp,dens,2,8933
mp,dens,3,7850
lsel,s,,,5,32,1
dl,all,6,VX,0,1   !定义无滑移边界条件
dl,all,6,VY,0,1  
/solu
allsel,all
lsel,s,,,1,4,1          ！定义外壳外表面对流换热系数
sfl,all,conv,5,,40
LDREAD,HGEN,,, ,0,'mxc','rst',' '  !读焦尔热  
!定义流体性质 (由LOG文件得到的）
FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,4
FLDATA13,VARY,SPHT,0
!*   
FLDATA7,PROT,DENS,AIR-SI
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0  
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,AIR-SI
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0  
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS   
FLDATA7,PROT,COND,AIR-SI
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0  
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,AIR-SI
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0  
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
  
!定义参考压力，名义温度，环境温度
  
FLDATA15,PRES,REFE,101350,   
FLDATA14,TEMP,NOMI,40,  
FLDATA14,TEMP,BULK,40,   
TOFFSET,273,
ACEL,0,9.8,0,   ！定义重力加速度
  
!设置为稳态层流，不可压缩，求解温度、流场方程
FLDATA1,SOLU,TRAN,0  
FLDATA1,SOLU,FLOW,1  
FLDATA1,SOLU,TEMP,1  
FLDATA1,SOLU,TURB,0  
FLDATA1,SOLU,COMP,0  
FLDATA1,SOLU,VOF,0   
FLDATA1,SOLU,SFTS,0  
FLDATA1,SOLU,IVSH,0  
FLDATA1,SOLU,SWRL,0  
FLDATA1,SOLU,SPEC,0  
FLDATA1,SOLU,ALE,0   
FLDATA1,SOLU,RDSF,0
!设置求解器为三对角求解器
!FLDATA18,METH,PRES,1
!FLDATA19,TDMA,PRES,100,  
!FLDATA18,METH,TEMP,2
!FLDATA19,TDMA,TEMP,100,  
FLDATA18,METH,TEMP,2
FLDATA22,MAXI,TEMP,1000,
FLDATA20,SRCH,TEMP,2,   
FLDATA21,CONV,TEMP,1e-012,   
FLDATA23,DELT,TEMP,1e-010,   
!迭代100次（其中一些设置是默认的）
/COM,,Steady State Analysis,0   
FLDATA2,ITER,EXEC,300,
FLDATA2,ITER,OVER,30,   
FLDATA2,ITER,APPE,30,   
FLDATA3,TERM,VX,0.01,   
FLDATA3,TERM,VY,0.01,   
FLDATA3,TERM,VZ,0.01,   
FLDATA3,TERM,PRES,1e-008,   
FLDATA3,TERM,TEMP,1e-008,   
FLDATA3,TERM,ENKE,0.01,  
FLDATA3,TERM,ENDS,0.01,  
FLDATA5,OUTP,SUMF,10,   
allsel,all
SOLVE  
多谢了

cmjssl

帖子发了这么多天了，为什么没有回应呀，高手都哪里去了，各位大侠多多指教呀

jnjl

最近这两天机子被别人占着，一直不能进行调试，你也不要着急，自己慢慢改改，大家也都帮你看呢。
你定义外表面换热系数为什么定义为5呢？我觉得是不是有些小呢？当然这是题外话，对你的结果可能不会有太大的影响。
lsel,s,,,1,4,1 ！定义外壳外表面对流换热系数  
sfl,all,conv,5,,40
其他的问题还没有看出来。。。

cmjssl

谢谢您的回复,我再研究研究,你也帮我仔细看看,真的非常感谢您,希望我们能成为朋友

cmjssl

我现在求解的三维流场温度不出现负值了,但我选择和二维同样的求解器不可以,压力和温度求解器都选择三对角求解器(第一种)求出的结果趋势和二维求出的不一样(用的是第二种),有高手知道二维和三维流场在选择求解器方面的不同吗,希望指点一二,多谢了!急呀

callousman

      按我的理解：求解器的不同选择可能会导致结果的差异，但不至于出现结果的变化趋势不同。同意jnjl的观点。应该是设置的问题。
      我原来求解流场时，有些区域也出现了负值。但是，在我将原来的常温边界条件（一个是常温，另一个是常热通量）改为常热通量时，就没有负值了，具体什么原因请教高手！

cmjssl

高手快来呀，我急呀

jiajinghua

希望可以和楼主交流下， QQ1213908066。谢谢！！