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

Automatic boundary conditions for 3 PTHEM gas diffusion flame model. More...

Functions/Subroutines

subroutine d3ptcl (itypfb, izfppp, icodcl, rcodcl)
 

Detailed Description

Automatic boundary conditions for 3 PTHEM gas diffusion flame model.

Function/Subroutine Documentation

◆ d3ptcl()

subroutine d3ptcl ( integer, dimension(nfabor)  itypfb,
integer, dimension(nfabor)  izfppp,
integer, dimension(nfabor,nvar icodcl,
double precision, dimension(nfabor,nvar,3)  rcodcl 
)
Parameters
[in]nvartotal number of variables
[in]itypfbboundary face types
[out]izfpppboundary face zone number
[out]icodclboundary 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)
[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) \gradt \, \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} $