My Project
programmer's documentation
Functions/Subroutines
cs_ctwr_bcond.f90 File Reference

Automatic boundary condition for cooling towers. More...

Functions/Subroutines

subroutine cs_ctwr_bcond (itypfb, izfppp, icodcl, rcodcl)
 

Detailed Description

Automatic boundary condition for cooling towers.

Function/Subroutine Documentation

◆ cs_ctwr_bcond()

subroutine cs_ctwr_bcond ( integer, dimension(nfabor)  itypfb,
integer, dimension(nfabor)  izfppp,
integer, dimension(nfabor,nvar icodcl,
double precision, dimension(nfabor,nvar,3)  rcodcl 
)
Parameters
[in]itypfbboundary face types
[in]izfpppzone number for the boundary face for the specific physic module
[in,out]icodclface boundary condition code:
  • 1 Dirichlet
  • 2 Radiative outlet
  • 3 Neumann
  • 4 sliding and $ \vect{u} \cdot \vect{n} = 0 $
  • 5 smooth wall and $ \vect{u} \cdot \vect{n} = 0 $
  • 6 rough wall and $ \vect{u} \cdot \vect{n} = 0 $
  • 9 free inlet/outlet (input mass flux blocked to 0)
  • 13 Dirichlet for the advection operator and Neumann for the diffusion operator
[in,out]rcodclvalue of the boundary conditions to edge faces
                         boundary condition values:
                          - rcodcl(1) value of the dirichlet
                          - rcodcl(2) value of the exterior exchange
                          -  coefficient (infinite if no exchange)
                          -  rcodcl(3) value flux density
                          -  (negative if gain) \f$w.m^{-2} \f$ or
                          -  roughness in \f$m\f$ if  icodcl=6
                           -# for velocity:
                                      \f$(\mu+\mu_T)\gradv \vect{u}\f$
                           -# for pressure: \f$ \Delta \grad P
                                            \cdot \vect{n} \f$
                           -# for scalar:   \f$ C_p \left ( K +
                                            \dfrac{K_T}{\sigma_T} \right)
                                            \grad T \cdot \vect{n} \f$