Algorithm for rapid integration of turbulence model equations on parabolic regions

While developing a three-dimensional boundary layer program using a standard parabolic matching scheme, the author has found computing time with the Wilcox-Rubesin (1979) two-equation turbulence model to be very lengthy. The long computing time occurs because converged solutions are possible only when very small streamwise steps are taken. The proposed remedy reduces computing time by increasing the maximum permissible step size.