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

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]	itypfb	boundary face types
[in]	izfppp	zone number for the boundary face for the specific physic module
[in,out]	icodcl	face 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]	rcodcl	value 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$