Boundary conditions for rough walls (icodcl = 6). More...

Functions/Subroutines
subroutine	clptrg (nscal, isvhb, icodcl, rcodcl, velipb, rijipb, visvdr, hbord, theipb)

subroutine	clptrg_scalar (iscal, isvhb, icodcl, rcodcl, byplus, buk, buet, bcfnns, hbord, theipb, tetmax, tetmin, tplumx, tplumn)

Detailed Description

Boundary conditions for rough walls (icodcl = 6).

The wall functions may change the value of the diffusive flux.

The values at a boundary face $\fib$ stored in the face center $\centf$ of the variable $ P $ and its diffusive flux $ Q $ are written as:

$P_\centf = A_P^g + B_P^g P_\centi$

and

$Q_\centf = A_P^f + B_P^f P_\centi$

where $P_\centi$ is the value of the variable $ P $ at the neighboring cell.

Warning:

for a vector field such as the velocity $\vect{u}$ the boundary conditions may read:
$\vect{u}_\centf = \vect{A}_u^g + \tens{B}_u^g \vect{u}_\centi$
and
$\vect{Q}_\centf = \vect{A}_u^f + \tens{B}_u^f \vect{u}_\centi$
where $\tens{B}_u^g$ and $\tens{B}_u^f$ are 3x3 tensor matrix which coupled veclocity components next to a boundary.

Please refer to the clptrg section of the theory guide for more informations.

Function/Subroutine Documentation

subroutine clptrg	(	integer	nscal,
		integer	isvhb,
		integer, dimension(:,:), pointer	icodcl,
		double precision, dimension(:,:,:), pointer	rcodcl,
		double precision, dimension(:,:)	velipb,
		double precision, dimension(:,:), pointer	rijipb,
		double precision, dimension(:), pointer	visvdr,
		double precision, dimension(:), pointer	hbord,
		double precision, dimension(:), pointer	theipb
	)

Parameters

[in]	nscal	total number of scalars
[in]	isvhb	indicator to save exchange coeffient
[in,out]	icodcl	face boundary condition code: 1 Dirichlet 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)
[in,out]	rcodcl	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) in w/m2 or roughness in m if icodcl=6 for the velocity $(\mu+\mu_T) \gradv \vect{u} \cdot \vect{n}$ for the pressure $\Delta t \grad P \cdot \vect{n}$ for a scalar $cp \left( K + \dfrac{K_T}{\sigma_T} \right) \grad T \cdot \vect{n}$
[in]	velipb	value of the velocity at $\centip$ of boundary cells
[in]	rijipb	value of $R_{ij}$ at $\centip$ of boundary cells
[out]	visvdr	viscosite dynamique ds les cellules de bord apres amortisst de v driest
[out]	hbord	coefficients d'echange aux bords
[in]	theipb	boundary temperature in $\centip$ (more exaclty the energetic variable)

subroutine clptrg_scalar	(	integer	iscal,
		integer	isvhb,
		integer, dimension(:,:), pointer	icodcl,
		double precision, dimension(:,:,:), pointer	rcodcl,
		double precision, dimension(:)	byplus,
		double precision, dimension(:)	buk,
		double precision, dimension(:)	buet,
		double precision, dimension(:)	bcfnns,
		double precision, dimension(:), pointer	hbord,
		double precision, dimension(:), pointer	theipb,
		double precision	tetmax,
		double precision	tetmin,
		double precision	tplumx,
		double precision	tplumn
	)

Parameters

[in]	iscal	scalar id
[in]	isvhb	indicator to save exchange coeffient
[in,out]	icodcl	face boundary condition code: 1 Dirichlet 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)
[in,out]	rcodcl	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) in w/m2 or roughness in m if icodcl=6 for the velocity $(\mu+\mu_T) \gradv \vect{u} \cdot \vect{n}$ for the pressure $\Delta t \grad P \cdot \vect{n}$ for a scalar $cp \left( K + \dfrac{K_T}{\sigma_T} \right) \grad T \cdot \vect{n}$
[in]	byplus	dimensionless distance to the wall
[in]	buk	dimensionless velocity
[in]	buet	boundary ustar value
[in]	bcfnns	boundary correction factor
[in,out]	hbord	exchange coefficient at boundary
[in]	theipb	boundary temperature in $\centip$ (more exaclty the energetic variable)
[out]	tetmax	maximum local ustar value
[out]	tetmin	minimum local ustar value
[out]	tplumx	maximum local tplus value
[out]	tplumn	minimum local tplus value