Solution of a hyperbolic system of turbulence-model equations by the “box” scheme

Abstract The momentum and continuity equations for a two-dimensional boundary layer, together with an empirical first-order partial differential equation for shear-stress transport, form a hyperbolic system of equations for u , v and τ. Here the Bradshaw-Ferriss-Atwell version of this turbulence model is solved by the Keller-Cebeci “box” scheme, which is particularly suited to systems of equations that are individually of first order. Computing time is about equal to that taken by a method-of-characteristics program if the same number of grid points are used across the layer and in the streamwise direction. However, the box scheme allows larger x -steps to be taken in the streamwise direction leading to smaller computing times.