Translation of the boundary conditions given by cs_user_boundary_conditions in a form that fits to the solver. More...

Functions/Subroutines
subroutine	condli

subroutine	set_dirichlet_scalar

subroutine	set_dirichlet_vector

subroutine	set_neumann_scalar

subroutine	set_neumann_vector

subroutine	set_convective_outlet_scalar

subroutine	set_convective_outlet_vector

Detailed Description

Translation of the boundary conditions given by cs_user_boundary_conditions in a form that fits to the solver.

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:

if we consider an increment of a variable, the boundary conditions read:
$\delta P_\centf = B_P^g \delta P_\centi$
and
$\delta Q_\centf = B_P^f \delta P_\centi$
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

subroutine condli ( )

Parameters

[in]	nvar	total number of variables
[in]	nscal	total number of scalars
[in]	isvhb	indicator to save exchange coeffient at the walls
[in]	iterns	iteration number on Navier-Stokes equations
[in]	isvtb	indicator to save the temperature at the walls
[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 rought wall and $\vect{u} \cdot \vect{n} = 0$ 9 free inlet/outlet (input mass flux blocked to 0)
[in,out]	isostd	indicator for standard outlet and reference face index
[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}$
[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)
[in]	frcxt	external force responsible for the hydrostatic pressure

Here is the call graph for this function:

subroutine set_convective_outlet_scalar ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	cofaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	cofbf	implicit BC coefficient for diffusive flux
[in]	pimp	Flux value to impose
[in]	cfl	Local Courant number used to convect
[in]	hint	Internal exchange coefficient

subroutine set_convective_outlet_vector ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	cofaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	cofbf	implicit BC coefficient for diffusive flux
[in]	pimpv	Dirichlet value to impose
[in]	cflv	Local Courant number used to convect
[in]	hint	Internal exchange coefficient

subroutine set_dirichlet_scalar ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	coefaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	coefbf	implicit BC coefficient for diffusive flux
[in]	pimp	Dirichlet value to impose
[in]	hint	Internal exchange coefficient
[in]	hext	External exchange coefficient (10^30 by default)

subroutine set_dirichlet_vector ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	coefaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	coefbf	implicit BC coefficient for diffusive flux
[in]	pimpv	Dirichlet value to impose
[in]	hint	Internal exchange coefficient
[in]	hextv	External exchange coefficient (10^30 by default)

subroutine set_neumann_scalar ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	coefaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	coefbf	implicit BC coefficient for diffusive flux
[in]	qimp	Flux value to impose
[in]	hint	Internal exchange coefficient

subroutine set_neumann_vector ( )

Parameters

[out]	coefa	explicit BC coefficient for gradients
[out]	coefaf	explicit BC coefficient for diffusive flux
[out]	coefb	implicit BC coefficient for gradients
[out]	coefbf	implicit BC coefficient for diffusive flux
[in]	qimpv	Flux value to impose
[in]	hint	Internal exchange coefficient

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation