flotran 旋转壁面困扰的问题
本帖最后由 battlestar 于 2012-9-28 15:34 编辑在用flotran做含有旋转壁面的流场分析时，我将fluid142单元的坐标系换成柱坐标系，在下图的moving wall中加载角速度vy=188rad/sec，外面的的wall三个方向速度均为零，下端设置成入口vz=1.5，上部为出口压力pres=0；然后用k-e求解，考虑重力加速度，所得流场速度比moving wall的线速度还大，不知道哪里出了问题，望高手指点！！
速度场向量图
下面是我的命令流文件，希望大家给于解答啊~
fini
/clea
s_r=0.45
s_h=0.5
f_r=0.408
f_h=0.43
a_r=0.146
a_h=0.7
as_r=0.155
/prep7
TOFFST,273,
csys,4
et,,142!流体3D单元
et,,226!耦合场单元
cylin,f_r,,f_h
wpof,,,-(s_h-f_h)/2
cylin,s_r,,s_h
cylin,a_r,,s_h
wpof,,,-(a_h-s_h)/2
cylin,a_r,,a_h
cylin,as_r,,a_h
vovl,all
wpro,,90
vsbw,all
wpro,,,90
vsbw,all
kwpa,1
wpro,,90
vsbw,all
kwpa,7
vsbw,all
FLST,5,28,6,ORDE,6
FITEM,5,6
FITEM,5,-17
FITEM,5,37
FITEM,5,-48
FITEM,5,57
FITEM,5,-60
VSEL,s, , ,P51X
cm,solid,volu
vatt,,,2
vsel,inve
cm,fluid,volume
FLST,5,4,5,ORDE,4
FITEM,5,188
FITEM,5,191
FITEM,5,194
FITEM,5,196
ASEL,S, , ,P51X
cm,input,area
FLST,5,4,5,ORDE,4
FITEM,5,198
FITEM,5,201
FITEM,5,204
FITEM,5,206
ASEL,S, , ,P51X
cm,output,area
cmse,s,fluid
aslv
cmse,u,input
cmse,u,output
csys,0
FLST,5,36,5,ORDE,30
FITEM,5,13
FITEM,5,36
FITEM,5,43
FITEM,5,51
FITEM,5,94
FITEM,5,98
FITEM,5,101
FITEM,5,106
FITEM,5,110
FITEM,5,114
FITEM,5,116
FITEM,5,120
FITEM,5,123
FITEM,5,126
FITEM,5,128
FITEM,5,131
FITEM,5,140
FITEM,5,142
FITEM,5,144
FITEM,5,146
FITEM,5,149
FITEM,5,-150
FITEM,5,152
FITEM,5,-153
FITEM,5,174
FITEM,5,-176
FITEM,5,178
FITEM,5,-180
FITEM,5,182
FITEM,5,-187
ASEL,s, , ,P51X
cm,interface,area
FLST,5,28,5,ORDE,28
FITEM,5,2
FITEM,5,6
FITEM,5,11
FITEM,5,16
FITEM,5,30
FITEM,5,33
FITEM,5,38
FITEM,5,41
FITEM,5,45
FITEM,5,49
FITEM,5,54
FITEM,5,58
FITEM,5,189
FITEM,5,192
FITEM,5,195
FITEM,5,197
FITEM,5,199
FITEM,5,202
FITEM,5,205
FITEM,5,207
FITEM,5,209
FITEM,5,-210
FITEM,5,213
FITEM,5,-214
FITEM,5,217
FITEM,5,-218
FITEM,5,220
FITEM,5,-221
ASEL,s, , ,P51X
cm,wall,area
alls
cmse,s,fluid
SMRT,1
vswe,all
cmse,s,solid
SMRT,6
vswe,all
alls
!流场
et,1,142
KEYOPT,1,3,3
et,2,0
!FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,TEMP,1 !考虑温度的流体计算
FLDATA1,SOLU,TURB,1 !湍流模型
FLDATA2,ITER,EXEC,50
FLDATA5,OUTP,TAUW,1 !考虑避免剪切力
FLDATA7,PROT,DENS,CONSTANT!流体密度
FLDATA8,NOMI,DENS,712,
FLDATA7,PROT,VISC,CONSTANT!流体粘度
FLDATA8,NOMI,VISC,9.12e-5,
FLDATA12,PROP,IVIS
FLDATA7,PROT,COND,CONSTANT!流体热导率
FLDATA8,NOMI,COND,53.96e-2,
FLDATA7,PROT,SPHT,CONSTANT!流体比热容
FLDATA8,NOMI,SPHT,5.376,
FLDATA4A,STEP,APPE,20,
FLDATA4,TIME,STEP,0.01,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,200
FLDATA,ALGR,SEGR,SIMPLEN
!流体分析环境边界条件
Da,wall,VX,0,1
Da,wall,VY,0,1
Da,wall,VZ,0,1
Da,interface,VX,0,1
Da,interface,VY,188,1
Da,interface,VZ,0,1
!DA,interface,ENKE,-1,1
CMSE,S,INPUT
nsla,,1
D,ALL,TEMP,300
Da,all,VX,0,1
Da,all,Vy,0,1
Da,all,VZ,1.5,1
CMSE,S,OUTPUT
nsla,,1
D,ALL,PRES,17.5e6,1
alls
dtran
/SOLU
ACEL,0,0,9.8,
solv
页:
[1]