My Project
programmer's documentation
Functions/Subroutines
typecl.f90 File Reference

Handle boundary condition type code (itypfb). More...

Functions/Subroutines

subroutine typecl (nvar, nscal, iterns, init, itypfb, itrifb, icodcl, isostd, rcodcl)
 

Detailed Description

Handle boundary condition type code (itypfb).

Function/Subroutine Documentation

◆ typecl()

subroutine typecl ( integer  nvar,
integer  nscal,
integer  iterns,
logical  init,
integer, dimension(ndimfb)  itypfb,
integer, dimension(ndimfb)  itrifb,
integer, dimension(ndimfb,nvar icodcl,
integer, dimension(ndimfb+1)  isostd,
double precision, dimension(ndimfb,nvar,3)  rcodcl 
)
Parameters
[in]nvartotal number of variables
[in]nscaltotal number of scalars
[in]iternsiteration number on Navier-Stokes equations
[in]initpartial treatment (before time loop) if true
[in,out]itypfbboundary face types
[out]itrifbtab d'indirection pour tri des faces
[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
[out]isostdstandard output indicator
  • reference face number
[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} $