我做的一个热流耦合的算例，请大家看看是否合理

mwei

计算的是流体绕流固体的情况
A1：固体，导热系数为10
A2：流体，导热系数、比热、粘性、密度均为5
边界条件：
L1：热流密度=100
L5，L6：VY=0
L9：Vx=50，Vy=0，T=10
L7：P=0
L8：自由边界

mwei

图

mwei

温度场结果：

mwei

速度场:

mwei

log文件,还没有时间整理出来

/BATCH  
/COM,ANSYS RELEASE  8.0    UP20030930       09:11:02    04/09/2005
!*  
/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   
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,2,,10
K,1,0,0,,   
K,2,0,0.04,,
K,3,0,0.06,,
K,4,0,0.1,,
K,5,0.1,0,,
K,6,0.04,0.02,,
K,7,0.06,0.02,,
K,7,0.06,0.02,,
kdel,2,3
K,2,0.02,0.04,,
K,2,0.04,0.02,,
K,3,0.06,0.02,,
K,8,0.1,0.1,,   
/replot
kdel,2,3
K,2,0.04,0,,
K,3,0.06,0,,
kplot   
FLST,2,4,3  
FITEM,2,2   
FITEM,2,3   
FITEM,2,7   
FITEM,2,6   
A,P51X  
FLST,2,8,3  
FITEM,2,1   
FITEM,2,2   
FITEM,2,6   
FITEM,2,7   
FITEM,2,3   
FITEM,2,5   
FITEM,2,8   
FITEM,2,4   
A,P51X  
AATT,       2, ,   1,       0,  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,1,,  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,1,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,VISC,1,,
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
SMRT,6  
MSHAPE,0,2D
MSHKEY,0
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,23,1,ORDE,12
FITEM,5,1   
FITEM,5,-4  
FITEM,5,9   
FITEM,5,-12
FITEM,5,31  
FITEM,5,-36
FITEM,5,41  
FITEM,5,-42
FITEM,5,45  
FITEM,5,-48
FITEM,5,60  
FITEM,5,-62
CM,_Y,NODE  
NSEL, , , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
NREFINE,_Y1, , ,1,1,1,1
CMDELE,_Y1  
!*  
elist   
FINISH  
/SOL
lplot   
FLST,2,1,4,ORDE,1   
FITEM,2,1   
/GO
!*  
SFL,P51X,HFLUX,100,100  
FLST,2,1,4,ORDE,1   
FITEM,2,9   
!*  
/GO
DL,P51X, ,TEMP,10,0
FLST,5,1,4,ORDE,1   
FITEM,5,9   
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
!*  
DL,_Y1, ,VX,50,1
DL,_Y1, ,VY,5,1
!*  
CMDELE,_Y1  
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,5   
FITEM,5,-6  
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
!*  
DL,_Y1, ,VY,0,1
!*  
CMDELE,_Y1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,9   
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
!*  
DL,_Y1, ,VX,50,1
DL,_Y1, ,VY,0,1
!*  
CMDELE,_Y1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,7   
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
/GO
!*  
DL,_Y1, ,PRES,0,1   
!*  
CMDELE,_Y1  
!*  
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,1
!*  
FLDATA,ALGR,SEGR,SIMPLEF
FLDATA,ALGR,HFLM,MATX   
!*  
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,CONSTANT  
FLDATA8,NOMI,DENS,5,
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,CONSTANT  
FLDATA8,NOMI,VISC,5,
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,5,
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,5,
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
!*  
FLDATA18,METH,VX,4  
FLDATA22,MAXI,VX,100,   
FLDATA20,SRCH,VX,12,
FLDATA20A,PGMR,FILL,6,  
FLDATA21,CONV,VX,1e-005,
FLDATA23,DELT,VX,1e-010,
!*  
FLDATA18,METH,VX,1  
FLDATA18,METH,TEMP,4
FLDATA22,MAXI,TEMP,1000,
FLDATA20,SRCH,TEMP,12,  
FLDATA20A,PGMR,FILL,6,  
FLDATA20A,PGMR,MODP,0,  
FLDATA21,CONV,TEMP,1e-012,  
FLDATA23,DELT,TEMP,1e-010,  
!*  

SOLVE   
SOLVE   
SOLVE   
FINISH  
/POST1  
/EFACE,1
!*  
PLNSOL,TEMP, ,0,
SET,LAST
/EFACE,1
!*  
PLNSOL,TEMP, ,0,
FINISH  
/SOL
FLDATA24,TURB,MODL,1
FLDATA24,TURB,CMU,0.09,
FLDATA24,TURB,C1,1.44,  
FLDATA24,TURB,C2,1.92,  
FLDATA24,TURB,SCTK,1,   
FLDATA24,TURB,SCTD,1.3,
!*  
SOLVE   
SOLVE   
/POST1  
/EFACE,1
!*  
PLNSOL,TEMP, ,0,
!*  
/VSCALE,1,1,0   
!   
!*  
PLVECT,V, , , ,VECT,ELEM,ON,0   
aplot   
/pnum,area,1
aplot   
lplot   
/pnum,line,1
lplot   
FINISH  
! /EXIT,NOSAV

mwei

看看有什么错误，欢迎大家回帖，
如果没有问题的话，解决共轭传热这类问题就可以这个干了
另，请问斑竹，我怎么可以挣到两个积分？刚好差2个积分：）

zwfhit

我想问一下，为什么你后来要运行下面的命令呢，你是怎么确定这些参数的？
/SOL
FLDATA24,TURB,MODL,1
FLDATA24,TURB,CMU,0.09,
FLDATA24,TURB,C1,1.44,
FLDATA24,TURB,C2,1.92,
FLDATA24,TURB,SCTK,1,
FLDATA24,TURB,SCTD,1.3,

mwei

先算层流再算湍流，用默认的参数试的

zwfhit

不过从你的命令来看，好像只是设置了湍流的有关参数，但没有打开湍流运算： FLDATA1,SOLU,TURB,1