We have developed a code to model the axisymmetric incompressible turbulent Navier-Stokes equations for chemical reactor problems. A conventional co-located finite volume method is used to discretize the conservative form of the continuum equations. The coupling between the pressure and velocity equations is achieved using the SIMPLE algorithm. The resulting iterative equations are solved using either the standard Tri-Diagonal Matrix Algorithm (TDMA) for the momentum and turbulence model equations or preconditioned conjugate gradient solvers (PCCG) for the pressure correction equation. The k-e and Reynolds stress (RSM) models of turbulence were used to achieve a closure of governing equations. A parallel version of the code was developed within the Communicating Sequential Process (CSP) programming model implemented within a message passing harness portable to both shared memory and distributed memory systems. The total parallel efficiency of the code is presented for both turbulence models implemented in the present code.
[1]
T. Gatski,et al.
Modelling the pressure–strain correlation of turbulence: an invariant dynamical systems approach
,
1991,
Journal of Fluid Mechanics.
[2]
J. H. Whitelaw,et al.
Velocity Characteristics of a Confined Coaxial Jet
,
1979
.
[3]
H. V. D. Vorst,et al.
Conjugate gradient type methods and preconditioning
,
1988
.
[4]
D. Kershaw.
The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations
,
1978
.
[5]
S. Patankar.
Numerical Heat Transfer and Fluid Flow
,
2018,
Lecture Notes in Mechanical Engineering.