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

Boundary condition code for the ALE module. More...

Functions/Subroutines

subroutine altycl (itypfb, ialtyb, icodcl, impale, init, dt, rcodcl, xyzno0)
 

Detailed Description

Boundary condition code for the ALE module.

Function/Subroutine Documentation

◆ altycl()

subroutine altycl ( integer, dimension(nfabor)  itypfb,
integer, dimension(nfabor)  ialtyb,
integer, dimension(nfabor,nvar icodcl,
integer, dimension(nnod)  impale,
logical  init,
double precision, dimension(ncelet)  dt,
double precision, dimension(nfabor,nvar,3)  rcodcl,
double precision, dimension(3,nnod)  xyzno0 
)
Parameters
[in]itypfbboundary face types
[in,out]ialtybboundary face types for ALE
[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]impaleimposed displacement indicator
[in]initpartial treatment (before time loop) if true
[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} $
[in]xyzno0initial mesh nodes coordinates