请教：自然对流问题

sophie_junjun

大家好，我想计算一个容器中液体在自然对流下的流动情况，可是计算结果完全不是我从文献中看到的结果，我实在不知道那错了，希望各位帮忙看看。如下图为加载后的图。上面的线为恒温1550度，左右两条线为恒定热流散热7769，下面那条线为恒定热流散热3778w/m2。中间液体为钢液，初始温度为1500，粘度0。006，导热系数46，比热787。

sophie_junjun

以下为设置求解部分的命令流：
FLDATA4,TIME,STEP,2,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,1.0e06
FLDATA4,TIME,TEND,10,   
FLDATA4,TIME,GLOB,20,   
FLDATA4,TIME,VX,0.01,   
FLDATA4,TIME,VY,0.01,   
FLDATA4,TIME,VZ,0.01,   
FLDATA4,TIME,PRES,1e-006,   
FLDATA4,TIME,TEMP,1e-006,   
FLDATA4,TIME,ENKE,0.01,
FLDATA4,TIME,ENDS,0.01,
FLDATA4A,STEP,OVER,0
FLDATA4,TIME,OVER,0,
FLDATA4A,STEP,APPE,0
FLDATA12,PROP,DENS,2

FLDATA13,VARY,DENS,1   密度可以随温度变化
FLDATA12,PROP,VISC,0
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0

FLDATA7,PROT,DENS,LIQUID
FLDATA8,NOMI,DENS,7000,
FLDATA9,COF1,DENS,300,  
FLDATA10,COF2,DENS,0.0001,  
FLDATA11,COF3,DENS,0,   
FLDATA7,PROT,VISC,CONSTANT  
FLDATA8,NOMI,VISC,0.006,
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,46,   
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,787,  
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0

ACEL,0,-9.81,0,
FLDATA15,PRES,REFE,101350,  
FLDATA16,BULK,BETA,1e+015,  
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,1500,
FLDATA14,TEMP,TTOT,293,
FLDATA14,TEMP,BULK,20,  
TOFFSET,273,

FLDATA18,METH,TEMP,1
FLDATA19,TDMA,TEMP,100,
FLDATA18,METH,PRES,1
FLDATA19,TDMA,PRES,100,
SOLVE

sophie_junjun

以下为设置求解部分的命令流：
FLDATA4,TIME,STEP,2,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,1.0e06
FLDATA4,TIME,TEND,10,   
FLDATA4,TIME,GLOB,20,   
FLDATA4,TIME,VX,0.01,   
FLDATA4,TIME,VY,0.01,   
FLDATA4,TIME,VZ,0.01,   
FLDATA4,TIME,PRES,1e-006,   
FLDATA4,TIME,TEMP,1e-006,   
FLDATA4,TIME,ENKE,0.01,
FLDATA4,TIME,ENDS,0.01,
FLDATA4A,STEP,OVER,0
FLDATA4,TIME,OVER,0,
FLDATA4A,STEP,APPE,0
FLDATA12,PROP,DENS,2

FLDATA13,VARY,DENS,1   密度可以随温度变化
FLDATA12,PROP,VISC,0
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0

FLDATA7,PROT,DENS,LIQUID
FLDATA8,NOMI,DENS,7000,
FLDATA9,COF1,DENS,300,  
FLDATA10,COF2,DENS,0.0001,  
FLDATA11,COF3,DENS,0,   
FLDATA7,PROT,VISC,CONSTANT  
FLDATA8,NOMI,VISC,0.006,
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,46,   
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,787,  
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0

ACEL,0,-9.81,0,
FLDATA15,PRES,REFE,101350,  
FLDATA16,BULK,BETA,1e+015,  
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,1500,
FLDATA14,TEMP,TTOT,293,
FLDATA14,TEMP,BULK,20,  
TOFFSET,273,

FLDATA18,METH,TEMP,1
FLDATA19,TDMA,TEMP,100,
FLDATA18,METH,PRES,1
FLDATA19,TDMA,PRES,100,
SOLVE

sophie_junjun

图形

sophie_junjun

自己顶一下，帮帮忙

sophie_junjun

请问一下，如果光对流体建模，是不是不能直接在流体边界加热通量和对流换热系数的边界条件呢？
我又计算了几遍，还是不对，很着急，帮帮忙行吗