src/base/clptrg.f90 File Reference

| programmer's documentation

- Main Page
- Related Pages
- Modules
- Data Structures
- Files

- File List
- Globals

All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
clptrg.f90 File Reference

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

## Functions/Subroutines
subroutine | clptrg (nscal, isvhb, icodcl, rtp, propce, rcodcl, velipb, rijipb, coefa, coefb, visvdr, hbord, theipb)

## Detailed Description

Boundary conditions for rought walls (icodcl = 6).

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

The values at a border face stored in the face center of the variable and its diffusive flux are written as:

and

where is the value of the variable at the neighbooring cell.

Warning:

- for a vector field such as the veclocity the boundary conditions may read:

  and

  where and are 3x3 tensor matrix which coupled veclocity components next to a boundary.

## Function/Subroutine Documentation

subroutine clptrg | ( | integer | nscal,
|  | integer | isvhb,
|  | integer, dimension(nfabor,nvarcl) | icodcl,
|  | double precision, dimension(ncelet,*) | rtp,
|  | double precision, dimension(ncelet,*) | propce,
|  | double precision, dimension(nfabor,nvarcl,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) | 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 5 smooth wall and 6 rought wall and 9 free inlet/outlet (input mass flux blocked to 0)
  [in] | rtp | calculated variables at cell centers (at current time step)
  [in] | propce | physical properties at cell 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 for the pressure for a scalar
  [in] | velipb | value of the velocity at of boundary cells
  [in] | rijipb | value of at 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] | theipb | boundary temperature in (more exaclty the energetic variable)