【分享】A FLOTRSN example

lanweiming

/batch
/prep7
!!  Note that for consistency, the energy unit must
!!  be in centiergs.  E+9 centiergs = 1 Joule
!!  Length in mm, mass in grams, temperature in C.
!!
!!   FLOTRAN solids can only rotate in a swirl problem
!!   Inner domain consists of only the inner solid
!!   Outer domain is the gap plus the outer solid
!!
!!  In this case, the velocity of interest is VZ, the swirl  
!!  velocity.
!!
!! Make the fluid domain
mat,1
idshaft = 20.0
odshaft = 40.
odouter = 60.0
gap = .015
yy = gap/5
slice = 0.2
velp = 20.
et,1,141,,,1
rect,odshaft/2,gap+odshaft/2,0, yy
nrs = 50
rs1 = -150
ng = 30
rg = -3
ny = 4
lsel,s,,,1,3,2
lesize,all,,,ng,rg
lsel,s,,,2,4,2
lesize,all,,,ny
alls
amesh,1
!! Make the Solid domain
!! Inner Solid is steel - area 2
!!   Use carbon steel properties
rect,idshaft/2.,odshaft/2.,0.,yy
mp,kxx,2,4.33E+7
mp,dens,2,7.8E-3
mp,c,2,1.0 ! (Arbitrary)
!! Outer Solid is lead bronze - area3
mp,kxx,3,4.0E+8
mp,dens,3,8.6E-3
mp,c,3,1.0 ! (Arbitrary)
rect,gap+odshaft/2.,odouter/2.,0,yy
nummrg,all
lsel,s,,,5,11,2
lesize,all,,,nrs,rs1
lsel,s,,,2,4,2
lsel,a,,,8,10,2
lesize,all,,,ny
alls
mat,2
amesh,2
mat,3
amesh,3
rpm = 250000
pi = 3.1415
thetav =(2* pi/60.)*rpm
!! Boundary conditions
dl,4,,enke,-1,1
dl,2,,vx,0,1
dl,2,,vy,0,1
dl,2,,vz,0,1
dl,10,,pres,0
!! Make the entire inner solid rotate
asel,s,,,2
nsla,,1
*get,MVN,node,0,count
*do,i,1,MVN
*get,thisnode,node,0,num,min
velz = thetav*NX(thisnode)
d,thisnode,vx,0
d,thisnode,vy,0
d,thisnode,vz,velz
nsel,u,,,thisnode
*enddo
alls
!! Temperature on the inside
dl,8,,temp,200,1
!! heat transfer coefficient on the outside
sfl,10,conv,2.E+6,,20
!!periodic boundary
lsel,s,,,1
lsel,a,,,5,9,4
dl,all,,vy,0,1
nsll,,1
peri,0,yy,0
alls
fini
/solu
!! Flotran commands
flda,iter,exec,50
flda,prot,dens,air-mm
flda,prot,visc,air-mm
flda,prot,cond,constant
flda,nomi,cond,2.6294E+4
flda,prot,spht,constant
flda,nomi,spht,1.004E+9
flda,temp,nomi,20.
toff,273
flda,solu,temp,t
flda,solu,ivsh,t
flda,solu,swrl,t
flda,relx,temp,1.0
flda,meth,temp,5
flda,meth,pres,5
flda,meth,vz,5
flda,advm,mome,supg
flda,advm,temp,supg
save
solve

jnjl

谢谢楼主的倾情奉献！
希望大家以后能多提供一些ｆｌｏｔｒａｎ的例子！

jnjl

将其存到一文本文档中，然后改名为*.inp文件
用/input,*,inp输入
你最好在问问题前仔细搜一下论坛中的贴子，这个问题已经有n多人提出来过了。。。
不过看在过节的份上，就再发一遍吧！