Mixing Length Model
From ThermalFluidsPedia
Line 1:  Line 1:  
{{Turbulence Category}}  {{Turbulence Category}}  
  [[Image:Fig4.33.pngthumb400 pxalt=Mixing length model   +  [[Image:Fig4.33.pngthumb400 pxalt=Mixing length model  Mixing length model.]] 
  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  +  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 the figure to the right, 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 <math>\bar{u}</math> according to the definition of mixing length. On the other hand, the timeaveraged velocity at point ''B'' is <math>\bar{u}+(\partial \bar{u}/\partial y)l</math> according to the profile of the timeaveraged velocity (see the figure). 
+  
+  Therefore, the packet must have a negative velocity fluctuation equal to <math>(\partial \bar{u}/\partial y)l</math> in order to keep its time averaged velocity unchanged. Thus, the fluctuation of the velocity component in the ''x''direction is:  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 9:  Line 11:  
 width="100%" <center>   width="100%" <center>  
<math>{u}'=l\left( \frac{\partial \bar{u}}{\partial y} \right)</math>  <math>{u}'=l\left( \frac{\partial \bar{u}}{\partial y} \right)</math>  
  
</center>  </center>  
{{EquationRef(1)}}  {{EquationRef(1)}}  
}  }  
  When the velocity component in the xdirection has the above negative fluctuation, the velocity component in the ydirection must have a positive fluctuation, <math>{v}'</math>, with the same scale, i.e.,  +  
+  When the velocity component in the ''x''direction has the above negative fluctuation, the velocity component in the ''y''direction must have a positive fluctuation, <math>{v}'</math>, with the same scale, i.e.,  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 23:  Line 25:  
{{EquationRef(2)}}  {{EquationRef(2)}}  
}  }  
  where C is a local constant. Thus, the timeaverage of the product of the velocity fluctuations, <math>\overline{{u}'{v}'}</math>, must be negative for this case.  +  where ''C'' is a local constant. Thus, the timeaverage of the product of the velocity fluctuations, <math>\overline{{u}'{v}'}</math>, 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 <math>{u}'</math> will be positive and <math>{v}'</math> will be negative. Therefore, <math>\overline{{u}'{v}'}</math> must be negative for any cases. Combining eqs. (  +  Similarly, we can also analyze motion of the fluid packet from point ''B'' to point ''A'', in which case <math>{u}'</math> will be positive and <math>{v}'</math> will be negative. Therefore, <math>\overline{{u}'{v}'}</math> must be negative for any cases. Combining eqs. (1) and (2) yields 
<math>\overline{{u}'{v}'}=Cl^{2}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}</math>  <math>\overline{{u}'{v}'}=Cl^{2}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}</math>  
  Since l is still undetermined, it will be beneficial to absorb C into l and yield  +  Since ''l'' is still undetermined, it will be beneficial to absorb ''C'' into ''l'' and yield 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 38:  Line 40:  
{{EquationRef(3)}}  {{EquationRef(3)}}  
}  }  
  It follows from the definition of the eddy diffusivity  +  It follows from the definition of the eddy diffusivity that 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 48:  Line 50:  
{{EquationRef(4)}}  {{EquationRef(4)}}  
}  }  
  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:  +  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: 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 57:  Line 59:  
{{EquationRef(5)}}  {{EquationRef(5)}}  
}  }  
  where <math>\kappa </math> is an empirical constant with order of 1, and is referred to as von Kármán’s constant. Equation (  +  where <math>\kappa </math> is an empirical constant with order of 1, and is referred to as von Kármán’s constant. Equation (5) is valid only if <math>\kappa </math> is really a constant. Substituting eq. (5) into eq. (4), the eddy diffusivity of momentum becomes 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 67:  Line 69:  
{{EquationRef(6)}}  {{EquationRef(6)}}  
}  }  
  Substituting eq. (  +  Substituting eq. (6) into eq. (4) of [[Algebraic Models for Eddy Diffusivity]], the shear stress in the twodimensional turbulent flow becomes 
{ class="wikitable" border="0"  { class="wikitable" border="0" 
Revision as of 02:07, 21 July 2010
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 the figure to the right, 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 the figure).
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. (1) and (2) 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 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:

where κ is an empirical constant with order of 1, and is referred to as von Kármán’s constant. Equation (5) is valid only if κ is really a constant. Substituting eq. (5) into eq. (4), the eddy diffusivity of momentum becomes

Substituting eq. (6) into eq. (4) of Algebraic Models for Eddy Diffusivity, the shear stress in the twodimensional turbulent flow becomes
