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

Functions/Subroutines
subroutine	clptur (nvar, nscal, isvhb, icodcl, dt, rtp, rtpa, propce, propfa, propfb, rcodcl, velipb, rijipb, coefa, coefb, visvdr, hbord, thbord)

Detailed Description

Boundary conditions for smooth walls (icodcl = 5).

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

The values at a border 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 neighbooring 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. This is only available when the option ivelco is set to 1.

Function/Subroutine Documentation

◆ clptur()

subroutine clptur	(	integer	nvar,
		integer	nscal,
		integer	isvhb,
		integer, dimension(nfabor,nvar)	icodcl,
		double precision, dimension(ncelet)	dt,
		double precision, dimension(ncelet,*)	rtp,
		double precision, dimension(ncelet,*)	rtpa,
		double precision, dimension(ncelet,*)	propce,
		double precision, dimension(nfac,*)	propfa,
		double precision, dimension(nfabor,*)	propfb,
		double precision, dimension(nfabor,nvar,3)	rcodcl,
		double precision, dimension(nfabor,ndim)	velipb,
		double precision, dimension(nfabor,6)	rijipb,
		double precision, dimension(nfabor,*)	coefa,
		double precision, dimension(nfabor,*)	coefb,
		double precision, dimension(ncelet)	visvdr,
		double precision, dimension(nfabor)	hbord,
		double precision, dimension(nfabor)	thbord
	)

Parameters

[in]	nvar	total number of variables
[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 rought wall and $\vect{u} \cdot \vect{n} = 0$ 9 free inlet/outlet (input mass flux blocked to 0)
[in]	dt	time step (per cell)
[in]	rtp,rtpa	calculated variables at cell centers (at current and previous time steps)
[in]	propce	physical properties at cell centers
[in]	propfa	physical properties at interior face centers
[in]	propfb	physical properties at boundary face centers
[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 roughtness 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]	coefa	explicit boundary condition coefficient
[out]	coefb	implicit boundary condition coefficient
[out]	visvdr	viscosite dynamique ds les cellules de bord apres amortisst de v driest
[out]	hbord	coefficients d'echange aux bords
[out]	thbord	boundary temperature in $\centip$ (more exaclty the energetic variable)

Here is the call graph for this function:

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ clptur()