programmer's documentation
Functions/Subroutines
cfxtcl.f90 File Reference

Initialisation of the variables if the compressible flow model is enabled. More...

Functions/Subroutines

subroutine cfxtcl (nvar, icodcl, itypfb, dt, rcodcl)
 

Detailed Description

Initialisation of the variables if the compressible flow model is enabled.

Handle boundary condition type code (itypfb) when the compressible model is enabled.

This subroutine is called at the beginning of a computation (or when a computation is resumed) before the start of the time loop.

It allows to initialise or modify (for resumed computations) the variables and the time step values.

Before this subroutine call, the density and the molecular viscosity have been initialised at ro0 and viscl0 respectively or they have been read in a checkpoint file in the case of a resumed computation. If the scalar diffusivities (visls) and the isobaric specific heat (cp) were defined (i.e. variable), their values are here at hand only if a computation is resumed.

Any modification of a physical property (density, molecular viscosity, scalar diffusivity, isobaric specific heat) shall be performed in the ppphyv subroutine and never here.

Function/Subroutine Documentation

subroutine cfxtcl ( integer  nvar,
integer, dimension(nfabor,nvarcl)  icodcl,
integer, dimension(nfabor)  itypfb,
double precision, dimension(ncelet)  dt,
double precision, dimension(nfabor,nvarcl,3)  rcodcl 
)
Parameters
[in]nvartotal number of variables
[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]itypfbboundary face types
[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} $