JCHLI2000 发表于 2005-11-12 21:01:44

我的流－热耦合，大家看看错在哪里？

做一个流－热耦合，对两个边界加热流，热流随位置的变化而变化，且不是一个连续的函数，我用一个表格数组来完成，当求解时总出现一下问题。

JCHLI2000 发表于 2005-11-12 21:03:42

文件的命令如下

文件的命令如下：
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,0
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
KEYW,PR_CFD,1   
/GO
!*
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,FLOTRAN CFD
!*
/PREP7
!*
ET,1,FLUID141
/PREP7
mptemp,1,375,600,1200,1485,1505,1540
mpdata,kxx,1,1,40,39,36,36,192,192
mpdata,dens,1,1,7500,7200,7000,6977,6970,6952   
mpdata,c,1,1,420,560,660,710,806,860
!mpdata,enth,1,1,0,1.6e9,5.04e9,6.81e9,9.98e9,1.054e10   
!*   
!*
RECTNG,0,0.15,0,0.9,
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,3   
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,20,-3, , , ,1
!*
FLST,5,2,4,ORDE,2   
FITEM,5,2   
FITEM,5,4   
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,20,4, , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,2   
CM,_Y1,LINE
LSEL, , , ,P51X
*GET,_z1,LINE,,COUNT
*SET,_z2,0
*DO,_z5,1,_z1   
*SET,_z2,LSNEXT(_z2)
*GET,_z3,LINE,_z2,ATTR,NDNX
*GET,_z4,LINE,_z2,ATTR,SPNX
*get,_z6,line,_z2,attr,kynd
*IF,_z3,GT,0,THEN   
*IF,_z4,NE,0,THEN   
LESIZE,_z2,,,_z3,1/_z4,,,,_z6   
*ENDIF
*ENDIF
*ENDDO
CMSEL,S,_Y1
CMDELE,_Y1
!*
MSHAPE,0,2D
MSHKEY,1
!*
CM,_Y,AREA
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y
!*
AMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
alls

FLST,5,2,4,ORDE,2   
FITEM,5,2   
FITEM,5,4   
LSEL,R, , ,P51X
nsll,s,1
NPLOT   
d,all,vx,0
d,all,vy,0
alls
LPLOT   

LSEL,R, , ,       1
nsll,s,1
nplot   
d,all,press,0   
alls
eplot   

FLST,5,3,1,ORDE,2   
FITEM,5,51
FITEM,5,-53
NSEL,R, , ,P51X
NPLOT   
d,all,vy,-0.6   
d,all,temp,1350
alls
NPLOT   
FLST,5,18,1,ORDE,5
FITEM,5,22
FITEM,5,42
FITEM,5,-50
FITEM,5,54
FITEM,5,-61
NSEL,R, , ,P51X
nplot   
d,all,vy,0
alls

alls
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*dim,tt,table,13,1,1,y,Hflux    !建立一个热流随Y位置变化的table
tt(1)=784000,1.8E6,1.814E6,1.608E6,1.34E6,1.065E6,1.017E6,997000!初始值
tt(9)=955000,893000,840000,818000,804000
tt(1,0)=0,0.029,0.064,0.086,0.107,0.16,0.236,0.3455
tt(9,0)=0.451,0.581,0.71,0.822,0.9
*vplot,tt(1,0),tt(1,1)    !显示热流随位置Y变化的曲线
FLST,5,2,4,ORDE,2      !选取定义热流的线
FITEM,5,2   
FITEM,5,4   
LSEL,s, , ,P51X
nsll,s,1
FINISH
/SOL
SF,all,HFLUX, %TT%   !给线上节点加载热流table

alls

/solu

TOFFST,273

FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,TURB,1
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,1
!*
FLDATA,ALGR,SEGR,SIMPLEN
FLDATA,ALGR,HFLM,MATX   
!*
/COM,,Transient Analysis,1
*SET,_z4,-999   
*SET,_z14,1

FLDATA12,PROP,DENS,0
FLDATA13,VARY,DENS,0
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,TABLE
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
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,TABLE
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,TABLE
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0

!*
FLDATA,ADVM,MOME,SUPG   
FLDATA,ADVM,TURB,MSU
FLDATA,ADVM,PRES,MSU
FLDATA,ADVM,TEMP,SUPG   


FLDATA4,TIME,STEP,0.05,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,500
FLDATA4,TIME,TEND,1.0e06
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,   
FLDATA4,TIME,APPE,1.0e6
FLDATA4A,STEP,SUMF,10,
FLDATA4,TIME,SUMF,1.0e6
FLDATA4,TIME,BC,0   
FLDATA4,TIME,TEND,1000000,
ACEL,0,-9.81,0,

solve

JCHLI2000 发表于 2005-11-12 21:06:54

Re:我的流－热耦合，大家看看错在哪里？

感觉热流也加上了，通过list/loads可以看出，热流也加上了，但不知什么原因出现上面的提示？

JCHLI2000 发表于 2005-11-15 06:39:31

Re:我的流－热耦合，大家看看错在哪里？

知道错在哪里了，哈哈

Amary 发表于 2005-11-15 12:46:06

Re:我的流－热耦合，大家看看错在哪里？

呵呵呵，很好玩，你在这里自问自答，我的问题也能这么解决就好了

JCHLI2000 发表于 2005-11-15 14:33:20

Re:我的流－热耦合，大家看看错在哪里？

解决是用函数来定义随位置变化的边界条件是可以了,但用网格数组来定义随位置而变化的边界条件还是会出现上面的情况,也不知怎么了?

epmkkll 发表于 2005-11-30 11:21:24

Re:我的流－热耦合，大家看看错在哪里？

flotran不允许加用表格定义的热边界条件
我想问一下怎么加用函数定义的热流边界条件，万分感谢！！！

JCHLI2000 发表于 2005-12-2 11:43:50

Re:我的流－热耦合，大家看看错在哪里？

在parameters/functions/difine and edit中定义

cookyzhao 发表于 2006-1-15 15:23:01

Re:我的流－热耦合，大家看看错在哪里？

修改后的文件，贴一个让大家瞻仰一下吧

xingkong52700 发表于 2006-4-27 12:17:01

请

麻烦楼主能不能把正确的程序给贴出来，大家瞻仰一下！

还有你的程序是不是中间有用GUI模式的，给出来吧，

这个程序有很多看不明白，不通俗！

谢谢！

maxiuihua 发表于 2006-4-30 22:27:54

页: [1]
查看完整版本: 我的流－热耦合，大家看看错在哪里？