Local variables to be added
integer ifac, ii
integer izone
integer ilelt, nlelt
double precision d2s3
double precision zref, xuref
double precision ustar, rugd, rugt
double precision zent, xuent, xvent
double precision xkent, xeent
integer, allocatable, dimension(:) :: lstelt
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
For boundary faces of color 11, assign an inlet boundary condition prescribed from the meteo profile with automatic choice between inlet/ outlet according to the meteo profile.
izone = 1
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
iprofm(izone) = 1
enddo
Example 2
For boundary faces of color 21, assign an inlet boundary condition prescribed from the meteo profile.
izone = 2
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
iprofm(izone) = 1
enddo
Example 3
For boundary faces of color 31, assign an inlet boundary condition prescribed from the meteo profile except for dynamical variables which are prescribed with a rough log law.
izone = 3
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
iprofm(izone) = 1
ustar=
xkappa*xuref/log((zref+rugd)/rugd)
xuent=ustar/
xkappa*log((zent+rugd)/rugd)
xvent = 0.d0
xeent=ustar**3/
xkappa/(zent+rugd)
rcodcl(ifac,
iu,1) = xuent
rcodcl(ifac,
iv,1) = xvent
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
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iep,1) = xeent
rcodcl(ifac,
iphi,1) = d2s3
rcodcl(ifac,
ifb,1) = 0.d0
rcodcl(ifac,
ik,1) = xkent
rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
rcodcl(ifac,
inusa,1) =
cmu*xkent**2/xeent
endif
enddo
Example 4
Prescribe at boundary faces of color '12' an outlet.
izone = 5
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
endif
enddo
Example 5
Prescribe at boundary faces of color 4 a symmetry.
izone = 6
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
enddo