Mixing Length Model
From ThermalFluidsPedia
External Turbulent Flow/Heat Transfer 
The mixing length model proposed by Prandtl is the simplest turbulent model. The distinctive feature of turbulent flow is the existence of eddies and vortices so that the transport in the turbulent flow is dominated by packets of the molecules, instead of the behavior of individual molecules. Considering a turbulent flow near a flat plat as shown in Fig. 1, the mixing length can be defined as the maximum length that a packet can travel vertically while maintaining its time averaged velocity unchanged. The concept of mixing length for turbulent flow is similar to the mean free path for random molecular motion. When a fluid packet located at point A travels to point B by moving upward a distance that equals the mixing length, l, its timeaveraged velocity should be kept at according to the definition of mixing length. On the other hand, the timeaveraged velocity at point B is according to the profile of the timeaveraged velocity (see Fig. 4.33). Therefore, the packet must have a negative velocity fluctuation equal to in order to keep its time averaged velocity unchanged. Thus, the fluctuation of the velocity component in the xdirection is:

When the velocity component in the xdirection has the above negative fluctuation, the velocity component in the ydirection must have a positive fluctuation, v', with the same scale, i.e.,

where C is a local constant. Thus, the timeaverage of the product of the velocity fluctuations, , must be negative for this case. Similarly, we can also analyze motion of the fluid packet from point B to point A, in which case u' will be positive and v' will be negative. Therefore, must be negative for any cases. Combining eqs. (4.410) and (4.411) yields
Since l is still undetermined, it will be beneficial to absorb C into l and yield

It follows from the definition of the eddy diffusivity [see eqs. (4.399) and (4.402)] that

where the absolute value is to ensure a positive eddy diffusivity. While the general rule for determining the mixing length, l, is lacking, the mixing length for turbulent boundary layer cannot exceed the distance to the wall. Therefore, we can assume:
l = κy 
where κ is an empirical constant with order of 1, and is referred to as von Kármán’s constant. Equation (4.414) is valid only if κ is really a constant. Substituting eq. (4.414) into eq. (4.413), the eddy diffusivity of momentum becomes

Substituting eq. (4.415) into eq. (4.402), the shear stress in the twodimensional turbulent flow becomes
