Local variables to be added
integer ifac, ii
integer izone
integer iclafu
integer ilelt, nlelt
double precision uref2, d2s3
double precision xkent, xeent
integer, allocatable, dimension(:) :: lstelt
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
The 12 color is a pure air inlet
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 1
ientat(izone) = 1
inmoxy(izone) = 1
qimpat(izone) = 1.46d-03
timpat(izone) = 400.d0 +
tkelvi
qimpfl(izone) = 0.d0
uref2 = rcodcl(ifac,
iu,1)**2 &
uref2 = max(uref2,1.d-12)
xkent, xeent )
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iep,1) = xeent
rcodcl(ifac,
ir11,1) = d2s3*xkent
rcodcl(ifac,
ir22,1) = d2s3*xkent
rcodcl(ifac,
ir33,1) = d2s3*xkent
rcodcl(ifac,
ir12,1) = 0.d0
rcodcl(ifac,
ir13,1) = 0.d0
rcodcl(ifac,
ir23,1) = 0.d0
rcodcl(ifac,
iep,1) = xeent
elseif (
iturb.eq.50)
then
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iep,1) = xeent
rcodcl(ifac,
iphi,1) = d2s3
rcodcl(ifac,
ifb,1) = 0.d0
elseif (
iturb.eq.60)
then
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
endif
endif
rcodcl(ifac,
isca(ii),1) = 1.d0
enddo
endif
enddo
Example 2
Inlet of both primary Air and Fuel
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 2
ientfl(izone) = 1
inmoxy(izone) = 1
qimpat(izone) = 1.46d-03
timpat(izone) = 800.d0 +
tkelvi
qimpfl(izone) = 1.46d-04/360.d0
iclafu = 1
distfu(izone,iclafu) = 100.d0
timpfl(izone) = 100.d0 +
tkelvi
uref2 = rcodcl(ifac,
iu,1)**2 &
uref2 = max(uref2,1.d-12)
xkent, xeent )
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iep,1) = xeent
rcodcl(ifac,
ir11,1) = d2s3*xkent
rcodcl(ifac,
ir22,1) = d2s3*xkent
rcodcl(ifac,
ir33,1) = d2s3*xkent
rcodcl(ifac,
ir12,1) = 0.d0
rcodcl(ifac,
ir13,1) = 0.d0
rcodcl(ifac,
ir23,1) = 0.d0
rcodcl(ifac,
iep,1) = xeent
elseif (
iturb.eq.50)
then
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iep,1) = xeent
rcodcl(ifac,
iphi,1) = d2s3
rcodcl(ifac,
ifb,1) = 0.d0
elseif (
iturb.eq.60)
then
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
endif
endif
enddo
Example 3
Color 15 is a wall
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 3
enddo
Example 4
Color 19 is an outlet
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 4
enddo
Example 5
14 and 4 are symmetry
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 5
enddo