Computation of three-dimensional turbulent boundary layers with heat transfer in a plane of symmetry using embedded wall-layer functions

In the calculation of turbulent boundary layers, a large number of mesh points are required to adequately resolve the intense variation in the velocity and enthalpy in the near-wall region. A substantial reduction in computational effort may be realized by representing the velocity and enthalpy profiles in the wall layer by analytical embedded functions. The effectively inviscid flow in the outer part of the boundary layer may then be resolved by employing a relatively coarser mesh. To obtain complete profiles, the outer numerical solution is matched asymptotically to the inner wall-layer analytical solution. To date, this approach has been restricted to two-dimensional flows; in the present study, a method which may be utilized for turbulent boundary layers with heat transfer in a plane of symmetry is developed as a first step in the application of the embedded-function method to full three-dimensional flows. The present method uses only about half as many mesh points as that required in a conventional procedure, which calculates the flow all the way to the wall, but there is no degradation in accuracy of the computed results.