7.2
general documentation
Functions/Subroutines
pptycl.f90 File Reference

Boundary conditions for specific physics modules. More...

Functions/Subroutines

subroutine pptycl (nvar, init, icodcl, itypfb, izfppp, dt, rcodcl)
 

Detailed Description

Boundary conditions for specific physics modules.

Function/Subroutine Documentation

◆ pptycl()

subroutine pptycl ( integer  nvar,
logical  init,
integer, dimension(nfabor,nvar)  icodcl,
integer, dimension(nfabor)  itypfb,
integer, dimension(nfabor)  izfppp,
double precision, dimension(ncelet)  dt,
double precision, dimension(nfabor,nvar,3)  rcodcl 
)
Parameters
[in]nvartotal number of variables
[in]initpartial treatment (before time loop) if true
[in,out]icodclface 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 rough wall and $ \vect{u} \cdot \vect{n} = 0 $
  • 9 free inlet/outlet (input mass flux blocked to 0)
  • 13 Dirichlet for the advection operator and Neumann for the diffusion operator
[in,out]itypfbboundary face types
[in,out]izfpppindex of the zone for the boundary faces (for the specific physics)
[in]dttime step (per cell)
[in,out]rcodclboundary 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
    1. for the velocity $ (\mu+\mu_T) \gradv \vect{u} \cdot \vect{n} $
    2. for the pressure $ \Delta t \grad P \cdot \vect{n} $
    3. for a scalar $ cp \left( K + \dfrac{K_T}{\sigma_T} \right) \grad T \cdot \vect{n} $