Finite difference method for solving turbulent boundary layer equations

A method is presented for numerically solving laminar, transitional and turbulent compressible boundary layer equations for axisymmetric flows with a three-point implicit finite difference scheme. The turbulent boundary layer is treated with a two layer concept, and at each layer appropriate eddy viscosity models are used to represent Reynolds stress terms. In order to relate the turbulent heat flux term with the Reynolds stress, the turbulent Prandtl number is assumed to be a constant. The system of governing equations is solved with an iteration method. For several situations, satisfactory results are obtained using an improved eddy viscosity model and a smooth method of the eddy viscosity distribution.