METEO 300
Fundamentals of Atmospheric Science

11.5 Here’s How Reynolds Did Averaging

For any variable, the observed value can be written as a sum of the mean value and a turbulent value:

u= u ¯ +u' MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadwhacqGH9aqpceWG1bGbaebacqGHRaWkcaWG1bGaai4jaaaa@3863@

 

where u ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaanaaabaGaamyDaaaaaaa@33D5@ is the mean, or average value, and u MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiqadwhagaqbaaaa@33D0@ is the turbulent part.

The average can be obtained by averaging u over time or over space or even by doing a number of samples and averaging over the samples.

temporal averaging:   u ¯ = u(t)dt dt = i=0 N1 u( t i ) N ;spatial averaging u ¯ = u(x)dx dx = i=0 N1 u( x i ) N MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaabshacaqGLbGaaeyBaiaabchacaqGVbGaaeOCaiaabggacaqGSbGaaeiiaiaabggacaqG2bGaaeyzaiaabkhacaqGHbGaae4zaiaabMgacaqGUbGaae4zaiaabQdacaqGGaGaaeiiaiqadwhagaqeaiabg2da9maalaaabaWaa8qaaeaacaWG1bGaaiikaiaadshacaGGPaGaamizaiaadshaaSqabeqaniabgUIiYdaakeaadaWdbaqaaiaadsgacaWG0baaleqabeqdcqGHRiI8aaaakiabg2da9maalaaabaWaaabCaeaacaWG1bGaaiikaiaadshadaWgaaWcbaGaamyAaaqabaGccaGGPaaaleaacaWGPbGaeyypa0JaaGimaaqaaiaad6eacqGHsislcaaIXaaaniabggHiLdaakeaacaWGobaaaiaacUdacaaMf8Uaae4CaiaabchacaqGHbGaaeiDaiaabMgacaqGHbGaaeiBaiaabccacaqGHbGaaeODaiaabwgacaqGYbGaaeyyaiaabEgacaqGPbGaaeOBaiaabEgacaaMf8UabmyDayaaraGaeyypa0ZaaSaaaeaadaWdbaqaaiaadwhacaGGOaGaamiEaiaacMcacaWGKbGaamiEaaWcbeqab0Gaey4kIipaaOqaamaapeaabaGaamizaiaadIhaaSqabeqaniabgUIiYdaaaOGaeyypa0ZaaSaaaeaadaaeWbqaaiaadwhacaGGOaGaamiEamaaBaaaleaacaWGPbaabeaakiaacMcaaSqaaiaadMgacqGH9aqpcaaIWaaabaGaamOtaiabgkHiTiaaigdaa0GaeyyeIuoaaOqaaiaad6eaaaaaaa@8EBB@

 

If the turbulence does not change with time and is homogeneous (i.e., the same in all directions and for all time), then these averages equal each other.

In Lesson 10, we developed the equation of motion without really considering the short-term and small-scale variations, except to say that they led to turbulent drag, which acts to resist the mean flow in the upper boundary layer. Now we want to think about how to correctly capture the dynamic effects of turbulent motion. What we want to do is to write down the equations of motion that you learned in Lesson 10; substitute mean and turbulent parts for the variables such as u, v, w; average over all the terms; and then see if we can sort out the terms to create an equation for the mean wind and an equation for the turbulent wind. This type of averaging is called Reynold’s averaging.

But first we need to learn the rules for averaging.

Rules of Averaging

c=constant; u,v=variables u ¯ =0 cu ¯ = c ¯ u ¯ u+v ¯ = u ¯ + v ¯ ( u ¯ v ) ¯ = u ¯ v ¯ ( du dt ) ¯ = d u ¯ dt MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOabaeqabaGaam4yaiabg2da9iaabogacaqGVbGaaeOBaiaabohacaqG0bGaaeyyaiaab6gacaqG0bGaae4oaiaabccacaqG1bGaaeilaiaabAhacaqG9aGaaeODaiaabggacaqGYbGaaeyAaiaabggacaqGIbGaaeiBaiaabwgacaqGZbaabaWaa0aaaeaaceWG1bGbauaaaaGaeyypa0JaaGimaaqaamaanaaabaGaam4yaiaaykW7caWG1baaaiabg2da9maanaaabaGaam4yaaaacaaMc8+aa0aaaeaacaWG1baaaaqaamaanaaabaGaamyDaiabgUcaRiaadAhaaaGaeyypa0JabmyDayaaraGaey4kaSIabmODayaaraaabaWaa0aaaeaadaqadaqaaiqadwhagaqeaiaaykW7caWG2baacaGLOaGaayzkaaaaaiabg2da9iqadwhagaqeaiaaykW7ceWG2bGbaebaaeaadaqdaaqaamaabmaabaWaaSaaaeaacaWGKbGaamyDaaqaaiaadsgacaWG0baaaaGaayjkaiaawMcaaaaacqGH9aqpdaWcaaqaaiaadsgaceWG1bGbaebaaeaacaWGKbGaamiDaaaaaaaa@6F83@

[11.1]
 

Now let’s apply these rules to a variable with a mean and a turbulent part. For example, consider the product uv, which is just the advection of horizontal wind in one direction by the wind in the other.

So, using the rules:

u ¯ = u ¯ +u' ¯ = u ¯ ¯ + u' ¯ , but u' ¯ =0,by definition so,  uv ¯ = ( u ¯ +u')( v ¯ +v') ¯ =...= u ¯ v ¯ + u ¯ v' ¯ + u' v ¯ ¯ + u'v' ¯ and u ¯ v' ¯ = u ¯ v' ¯ = u ¯ 0=0; u' v ¯ ¯ = u' ¯ v ¯ ¯ =0 v ¯ =0 so uv ¯ = u ¯ v ¯ + u'v' ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOabaeqabaWaa0aaaeaacaWG1baaaiabg2da9maanaaabaGabmyDayaaraGaey4kaSIaamyDaiaacEcaaaGaeyypa0Zaa0aaaeaaceWG1bGbaebaaaGaey4kaSYaa0aaaeaacaWG1bGaai4jaaaacaGGSaaabaGaamOyaiaadwhacaWG0bGaaGzbVpaanaaabaGaamyDaiaacEcaaaGaeyypa0JaaGimaiaacYcacaaMf8UaaeOyaiaabMhacaqGGaGaaeizaiaabwgacaqGMbGaaeyAaiaab6gacaqGPbGaaeiDaiaabMgacaqGVbGaaeOBaaqaaiaabohacaqGVbGaaeilaiaabccaaeaadaqdaaqaaiaadwhacaaMc8UaamODaaaacqGH9aqpdaqdaaqaaiaacIcaceWG1bGbaebacqGHRaWkcaWG1bGaai4jaiaacMcacaGGOaGabmODayaaraGaey4kaSIaamODaiaacEcacaGGPaaaaiabg2da9iaac6cacaGGUaGaaiOlaiabg2da9iqadwhagaqeaiaayIW7ceWG2bGbaebacaaMc8Uaey4kaSYaa0aaaeaaceWG1bGbaebacaWG2bGaai4jaaaacqGHRaWkdaqdaaqaaiaadwhacaGGNaGabmODayaaraaaaiabgUcaRmaanaaabaGaamyDaiaacEcacaWG2bGaai4jaaaaaeaacaWGHbGaamOBaiaadsgaaeaadaqdaaqaaiqadwhagaqeaiaaykW7caWG2bGaai4jaaaacqGH9aqpceWG1bGbaebacaaMc8+aa0aaaeaacaWG2bGaai4jaaaacqGH9aqpceWG1bGbaebacqGHflY1caaIWaGaeyypa0JaaGimaiaacUdacaaMf8+aa0aaaeaacaWG1bGaai4jaiqadAhagaqeaaaacqGH9aqpdaqdaaqaaiaadwhacaGGNaaaamaanaaabaGabmODayaaraaaaiabg2da9iaaicdacqGHflY1daqdaaqaaiaadAhaaaGaeyypa0JaaGimaaqaaiaadohacaWGVbaabaWaa0aaaeaacaWG1bGaaGPaVlaadAhaaaGaeyypa0JabmyDayaaraGaaGjcVlqadAhagaqeaiaaykW7cqGHRaWkdaqdaaqaaiaadwhacaGGNaGaamODaiaacEcaaaaaaaa@ADC3@

 

This second term, a product of two turbulent terms is not necessarily zero! Whereas the average of one turbulent term is zero by definition, the average of two turbulent terms is not necessarily zero.

The following video (3.11) further describes Reynold's averaging:

Click here for a transcript for Reynold's Averaging Example Video

Example

Consider two random numbers varying between -0.5 and 0.5, called u and v. The figure below shows u, v, the mean of uv and the mean of u’v’. Of course, the mean of u’v’ might be zero, but is it not necessarily zero, as shown in this figure.

graph of The product of two random numbers whose means are zero
The product of two random numbers whose means are zero. Note that the mean of u’v’ is not zero!
Credit: W. Brune

The same thinking applies to u2.

u 2 ¯ = ( u ¯ +u')( u ¯ +u') ¯ =...= u ¯ 2 + u ' 2 ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaanaaabaGaamyDaiaaykW7daahaaWcbeqaaiaaikdaaaaaaOGaeyypa0Zaa0aaaeaacaGGOaGabmyDayaaraGaey4kaSIaamyDaiaacEcacaGGPaGaaiikaiqadwhagaqeaiabgUcaRiaadwhacaGGNaGaaiykaaaacqGH9aqpcaGGUaGaaiOlaiaac6cacqGH9aqpceWG1bGbaebacaaMi8+aaWbaaSqabeaacaaIYaaaaOGaaGPaVlabgUcaRmaanaaabaGaamyDaiaacEcadaahaaWcbeqaaiaaikdaaaaaaaaa@4E12@

[11.2]
 

Remember your statistics and the concept of variance:

σ u 2 = 1 N1 i=0 N1 ( u i u ¯ ) 2 1 N i=0 N1 ( u i u ¯ ) 2 = 1 N i=0 N1 ( u i ' ) 2 = u ' 2 ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBaaaleaacaWG1baabeaakmaaCaaaleqabaGaaGOmaaaakiabg2da9maalaaabaGaaGymaaqaaiaad6eacqGHsislcaaIXaaaamaaqahabaWaaeWaaeaacaWG1bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IabmyDayaaraaacaGLOaGaayzkaaaaleaacaWGPbGaeyypa0JaaGimaaqaaiaad6eacqGHsislcaaIXaaaniabggHiLdGcdaahaaWcbeqaaiaaikdaaaGccqGHijYUdaWcaaqaaiaaigdaaeaacaWGobaaamaaqahabaWaaeWaaeaacaWG1bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IabmyDayaaraaacaGLOaGaayzkaaaaleaacaWGPbGaeyypa0JaaGimaaqaaiaad6eacqGHsislcaaIXaaaniabggHiLdGcdaahaaWcbeqaaiaaikdaaaGccqGH9aqpdaWcaaqaaiaaigdaaeaacaWGobaaamaaqahabaWaaeWaaeaacaWG1bWaaSbaaSqaaiaadMgaaeqaaOGaai4jaaGaayjkaiaawMcaaaWcbaGaamyAaiabg2da9iaaicdaaeaacaWGobGaeyOeI0IaaGymaaqdcqGHris5aOWaaWbaaSqabeaacaaIYaaaaOGaeyypa0Zaa0aaaeaacaWG1bGaai4jamaaCaaaleqabaGaaGOmaaaaaaaaaa@6D6A@

[11.3]
 

So, the variance is the same as the mean value for the square of the turbulent part of the variable.

The covariance of u and v is given by the equation:

covariance(u,v)= u'v' ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaabogacaqGVbGaaeODaiaabggacaqGYbGaaeyAaiaabggacaqGUbGaae4yaiaabwgacaGGOaGaamyDaiaacYcacaWG2bGaaiykaiabg2da9maanaaabaGaamyDaiaacEcacaWG2bGaai4jaaaaaaa@4463@

[11.4]
 

The standard deviation is just:

σ u = ( u ' 2 ¯ ) 1/2 MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBaaaleaacaWG1baabeaakiabg2da9maabmaabaWaa0aaaeaacaWG1bGaai4jamaaCaaaleqabaGaaGOmaaaaaaaakiaawIcacaGLPaaadaahaaWcbeqaaiaaigdacaGGVaGaaGOmaaaaaaa@3D4C@

[11.5]
 

We can get a better sense of how large this variance is by dividing by the mean value:

I =  σ u u ¯ MathType@MTEF@5@5@+=faaagCart1ev2aaaKnaaaaWenf2ys9wBH5garuavP1wzZbItLDhis9wBH5garmWu51MyVXgaruWqVvNCPvMCG4uz3bqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9vqqj=hHeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpi0dc9GqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaabMeacaqGGaGaaeypaiaabccadaWcaaqaaiabeo8aZnaaBaaaleaacaWG1baabeaaaOqaaiqadwhagaqeaaaaaaa@39B1@

[11.6]