The Flow/chart can be described in the following schematic diagrams (Figure 36.1)

MAIN: Main Program that calls all the modules and subroutines.
CONTI : Module for solving continuity equation
BCC: Boundary conditions for the confining surfaces
BCD: Boundary condition for obstacte
CEQCP: Continuity equation for constant property cases
DIV: Velocity divergence in each cell
E PSI: Pre-defined small numerical value
ITI: Iterative counter for the continuity equation
BCNS: Boundary conditions for the Navier-Stokes equations
VELALT: Changes the converged velocity arrays, i.e., makes the converged values at (n+1)th in time step as nth time step to start the next time step.
TICORR: Calculates δt
DTMAX: Maximum discrepancy of the dependent variables between two time steps
STAT: Predefined small numerical value
NSEQCP: Discretized Navier-Stokes equations for constant property
BCOV: Special Boundary conditions for obstacles
TIGRAD: Calculation of change in dependent variables between two time steps
ITA: Iterative counter for calculating the Navier-Stokes equations
ENERGY: Module for calculation of ENERGY equation
BCT: Boundary conditions for temperatures
TEQCP: Discretized Energy equation for constant property cases
DT: Maximum discrepancy of temperature in any cell between two iterative steps
TSTAT: Predefined small numerical value
TEMALT: Changes the temperature arrays from (m+1) iterative level m th iterative level to start the next iteration