Boundary conditions for smooth walls (icodcl = 5). More...

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

subroutine	clptur_scalar (iscal, isvhb, icodcl, rcodcl, byplus, bdplus, buk, hbord, theipb, tetmax, tetmin, tplumx, tplumn)

subroutine	clptur_vector (iscal, isvhb, icodcl, rcodcl, byplus, bdplus, buk, hbord)

Detailed Description

Boundary conditions for smooth walls (icodcl = 5).

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 veclocity $\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 wall boundary conditions section of the theory guide for more informations, as well as the clptur section.

Function/Subroutine Documentation

subroutine clptur	(	integer	nscal,
		integer	isvhb,
		integer, dimension(nfabor,nvar)	icodcl,
		double precision, dimension(nfabor,nvar,3)	rcodcl,
		double precision, dimension(nfabor,ndim)	velipb,
		double precision, dimension(nfabor,6)	rijipb,
		double precision, dimension(ncelet)	visvdr,
		double precision, dimension(nfabor)	hbord,
		double precision, dimension(nfabor)	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	dynamic viscosity after V. Driest damping in boundary cells
[out]	hbord	exchange coefficient at boundary
[in]	theipb	value of thermal scalar at $\centip$ of boundary cells

subroutine clptur_scalar	(	integer	iscal,
		integer	isvhb,
		integer, dimension(nfabor,nvar)	icodcl,
		double precision, dimension(nfabor,nvar,3)	rcodcl,
		double precision, dimension(nfabor)	byplus,
		double precision, dimension(nfabor)	bdplus,
		double precision, dimension(nfabor)	buk,
		double precision, dimension(nfabor)	hbord,
		double precision, dimension(nfabor)	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]	bdplus	dimensionless shift to the wall for scalable wall functions
[in]	buk	dimensionless velocity
[in,out]	hbord	exchange coefficient at boundary
[in]	theipb	value of thermal scalar at $\centip$ of boundary cells
[out]	tetmax	maximum local ustar value
[out]	tetmin	minimum local ustar value
[out]	tplumx	maximum local tplus value
[out]	tplumn	minimum local tplus value

subroutine clptur_vector	(	integer	iscal,
		integer	isvhb,
		integer, dimension(nfabor,nvar)	icodcl,
		double precision, dimension(nfabor,nvar,3)	rcodcl,
		double precision, dimension(nfabor)	byplus,
		double precision, dimension(nfabor)	bdplus,
		double precision, dimension(nfabor)	buk,
		double precision, dimension(nfabor)	hbord
	)

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]	bdplus	dimensionless shift to the wall for scalable wall functions
[in]	buk	dimensionless velocity
[in,out]	hbord	exchange coefficient at boundary