TwoLayer Model
From ThermalFluidsPedia
External Turbulent Flow/Heat Transfer 
The above mixing length theory assumes that the turbulent flow takes place within the entire boundary layer thickness. In reality, the turbulent boundary layer can be further divided into two regions (see figure to the right). The first region is viscous sublayer and it occupies less than 1% of the total turbulent boundary layer thickness. The momentum and heat transfer in this region are dominated by viscous shear and heat conduction, respectively. Outside of the viscous sublayer is a full turbulent region which comprises almost the entirety of the turbulent boundary layer. Figure (a) shows the instantaneous velocity profiles in a turbulent boundary layer on a flat plate taken at 17 instances. The 17 superimposed profiles as well as the average profile are shown in Fig. (b) and (c), respectively. It can be clearly seen from figure(c) that the velocity increases sharply in the sublayer near the wall and its change in the second region is fairly insignificant. While the thickness of the turbulent boundary layer increases with increasing x, the viscous sublayer remains at a fairly constant thickness. As the thickness of the turbulent boundary layer increases, the viscous sublayer inhabits a smaller and smaller portion of the entire turbulent boundary layer. Therefore, the mixing length model based on the assumption that the entire turbulent boundary layer is full turbulent flow should be improved by considering the twolayer structure of the turbulent boundary layer.
In the twolayer structure of the turbulent boundary layer, the velocity is dominated by viscous shear stress in the sublayer and by turbulent mixing in the fully turbulent region. Near the wall, the shear stress is equal to the shear stress at the wall, i.e., . Thus, eq. (4) in Algebraic Models for Eddy Diffusivity near the wall becomes

where the partial derivative has been approximated as because in the region near the wall. For the viscous sublayer, and eq. (1) is simplified to:

Equation (2) can be integrated from the wall to yield
i.e.,

Introducing the wall coordinate

where

is shear velocity or friction velocity, eq. (3) becomes

which is the velocity distribution in the viscous sublayer. For the fully turbulent region, and eq. (1) is simplified to:

which can be rewritten to the following format using the Prandtl’s mixing length model

Substituting eq. (4), eq. (8) can be nondimensionalized into

In order to integrate eq. (9), we must know the thickness of the viscous sublayer. For the case without blowing or suction on the wall and zeropressure gradient, we can choose y^{ + } = 10.8 as the dimensionless thickness of the viscous sublayer. Therefore, eq. (9) can be integrated from y^{ + } = 10.8, i.e.,

which can be integrated to yield

where the von Kármán’s constant is taken as κ = 0.4. Equation (11) is often referred to as the law of the wall. It can also be approximated into the following power law form:

which fits into eq. (11) at least up to y^{ + } = 1500.
The twolayer model can then be summarized as

Figure to the right shows the comparison between the above twolayer model and experimental results obtained using water and air. The Reynolds number shown in the figure is the momentum thickness Reynolds number defined as , where is the momentum thickness. It can be seen that the results from the twolayer model agreed with the experimental results very well except at the outer region of the turbulent boundary layer. In addition, the agreement between eq. (13) and the experimental results is also not very good in the region near y^{ + } = 10.8. Some researchers suggested there is a buffer region (5 < y^{ + } < 30) between the sublayer and the fully turbulent region and the velocity profile in the buffer region is ^{[1]}^{[2]}:
