Mixing length model

{{Short description|Mathematical model in fluid dynamics}}

Image:Mixing length.jpg will keep its original characteristics before dispersing them into the surrounding fluid. Here, the bar on the left side of the figure is the mixing length.]]

In fluid dynamics, the mixing length model is a method attempting to describe momentum transfer by turbulence Reynolds stresses within a Newtonian fluid boundary layer by means of an eddy viscosity. The model was developed by Ludwig Prandtl in the early 20th century.{{cite book |last=Holton |first=James R. |title=Dynamic Meteorology |edition=4th |series=International Geophysics Series |volume=88 |year=2004 |publisher=Elsevier Academic Press |location=Burlington, MA |pages=124–127 |chapter=Chapter 5 – The Planetary Boundary Layer }} Prandtl himself had reservations about the model,{{cite journal | last = Prandtl | first = L. | title = 7. Bericht über Untersuchungen zur ausgebildeten Turbulenz | journal = Z. Angew. Math. Mech. | volume = 5 | issue = 1 | pages = 136–139 | date = 1925 | doi = 10.1002/zamm.19250050212 | bibcode = 1925ZaMM....5..136P }} describing it as, "only a rough approximation,"{{cite journal | doi = 10.1038/249135b0 | last = Bradshaw | first = P. | title = Possible origin of Prandt's mixing-length theory | journal = Nature | volume = 249 | issue = 6 | pages = 135–136 | date = 1974 | bibcode = 1974Natur.249..135B | s2cid = 4218601 }}

but it has been used in numerous fields ever since, including atmospheric science, oceanography and stellar structure.{{cite journal | doi = 10.1126/science.235.4787.465 | last = Chan | first = Kwing |author2=Sabatino Sofia | title = Validity Tests of the Mixing-Length Theory of Deep Convection | journal = Science | volume = 235 | issue = 4787 | pages = 465–467 | date = 1987 | pmid = 17810341|bibcode = 1987Sci...235..465C | s2cid = 21960234 }} Also, Ali and Dey{{cite journal | last1 = Ali | first1 = S.Z. | last2 = Dey | first2 = S. | year = 2020 | title = The law of the wall: A new perspective |journal=Physics of Fluids | volume = 36| pages = 121401| doi = 10.1063/5.0036387 }} hypothesized an advanced concept of mixing instability.

Physical intuition

The mixing length is conceptually analogous to the concept of mean free path in thermodynamics: a fluid parcel will conserve its properties for a characteristic length, \ \xi' , before mixing with the surrounding fluid. Prandtl described that the mixing length,{{cite book | last = Prandtl | first = L. | title = Proc. Second Intl. Congr. Appl. Mech. | location = Zürich | date = 1926 }}

{{cquote|may be considered as the diameter of the masses of fluid moving as a whole in each individual case; or again, as the distance traversed by a mass of this type before it becomes blended in with neighbouring masses...}}

In the figure above, temperature, \ T, is conserved for a certain distance as a parcel moves across a temperature gradient. The fluctuation in temperature that the parcel experienced throughout the process is \ T'. So \ T' can be seen as the temperature deviation from its surrounding environment after it has moved over this mixing length \ \xi'.

Mathematical formulation

To begin, we must first be able to express quantities as the sums of their slowly varying components and fluctuating components.

=Reynolds decomposition=

This process is known as Reynolds decomposition. Temperature can be expressed as:{{cite news|url=https://web1.eng.famu.fsu.edu/~dommelen/courses/flm/flm00/topics/turb/node2.html |title=Reynolds Decomposition | publisher=Florida State University | date=6 December 2008 | access-date=2008-12-06}}

T = \overline{T} + T',

where \overline{T}, is the slowly varying component and T' is the fluctuating component.

In the above picture, T' can be expressed in terms of the mixing length considering a fluid parcel moving in the z-direction:

T' = -\xi' \frac{\partial \overline{T}}{\partial z}.

The fluctuating components of velocity, u', v', and w', can also be expressed in a similar fashion:

u' = -\xi' \frac{\partial \overline{u}}{\partial z}, \qquad \ v' = -\xi' \frac{\partial \overline{v}}{\partial z}, \qquad \ w' = -\xi' \frac{\partial \overline{w}}{\partial z}.

although the theoretical justification for doing so is weaker, as the pressure gradient force can significantly alter the fluctuating components. Moreover, for the case of vertical velocity, w' must be in a neutrally stratified fluid.

Taking the product of horizontal and vertical fluctuations gives us:

\overline{u' w'} = \overline{\xi' ^2} \left | \frac{\partial \overline{w}}{\partial z}\right| \frac{\partial \overline{u}}{\partial z}.

The eddy viscosity is defined from the equation above as:

K_m=\overline{\xi'^2} \left| \frac{\partial \overline{w}}{\partial z}\right|,

so we have the eddy viscosity, K_m expressed in terms of the mixing length, \xi'.

See also

References