Local variables to be added
integer ifac, izone, ii
integer ilelt, nlelt
double precision uref2, xkent, xeent, d2s3
integer, allocatable, dimension(:) :: lstelt
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
Definition of a burned gas inlet (pilot flame) for each face of colour 11
call getfbr(
'11', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 1
izfppp(ifac) = izone
rcodcl(ifac,
iw,1) = 21.47d0
uref2 = rcodcl(ifac,
iu,1)**2 &
uref2 = max(uref2,1.d-12)
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
elseif (
iturb.eq.70)
then
rcodcl(ifac,
inusa,1) =
cmu*xkent**2/xeent
endif
endif
enddo
Example 2
Definition of an unburned gas inlet for each face of colour 12
call getfbr(
'12', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 2
izfppp(ifac) = izone
rcodcl(ifac,
iu,1) = 60.d0
enddo
Example 3
Definition of a wall for each face of colour 51 and 5
call getfbr(
'51 or 5', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 4
izfppp(ifac) = izone
enddo
Example 4
Definition of an exit for each face of colour 91 and 9
call getfbr(
'91 or 9', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 5
izfppp(ifac) = izone
enddo
Example 5
Definition of symmetric boundary conditions for each face of colour 41 and 4.
call getfbr(
'41 or 4', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 6
izfppp(ifac) = izone
enddo