The examples provided here are only concern inlet boundary conditions. Other boundary condition types are not specific to the pulverized coal model, so generic examples apply (and as usual, the GUI is recommended).
Note that turbulence variables BC's are usually computed automatically:
- If icalke(izone) = 1: hydraulic diameter and reference velocity (similar to turbulence_bc_inlet_hyd_diam).
- If icalke(izone) = 2: hydraulic diameter, reference velocity and turbulence intensity (similar to turbulence_bc_inlet_turb_intensity).
- If icalke(izone) = 0: non-automatic definition. Turbulence BC's must be defined here for all turbulence variables (k, epsilon, omega, Rij, ... depending on active model).
Local variables to be added
integer ifac, ii
integer izone
integer icha, iclapc
integer ilelt, nlelt
integer, allocatable, dimension(:) :: lstelt
This firs example
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
Secondary or tertiary air inlet, for zone 'inlet'.
call getfbr(
'inlet', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 1
izfppp(ifac) = izone
rcodcl(ifac,
isca(ii),1) = 1.d0
enddo
endif
enddo
Example 2
Primary air inlet, for zone 'primary_inlet'.
call getfbr(
'primary_inlet', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 2
izfppp(ifac) = izone
enddo
enddo
icha = 1
enddo
enddo
Example 3
Secondary or tertiary air inlet, for zone 'inlet', using the pulverized coal with moisture model (with Lagrangian particle tracking).
call getfbr(
'inlet', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izone = 1
izfppp(ifac) = izone
rcodcl(ifac,
isca(ii),1) = 1.d0
enddo
endif
enddo