Introduction to turbulence/Free turbulent shear flows
From CFD-Wiki
(→The shear layer equations) |
m (fixed typo) |
||
(174 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
+ | {{Introduction to turbulence menu}} | ||
== Introduction == | == Introduction == | ||
- | Free shear flows are | + | Free shear flows are inhomogeneous flows with mean velocity gradients that develop in the absence of boundaries. Turbulent free shear flows are commonly found in natural and engineering environments. The jet of of air issuing from one's nostrils or mouth upon exhaling, the turbulent plume from a smoldering cigarette, and the buoyant jet issuing from an erupting volcano - all illustrate both the omnipresence of free turbulent shear flows and the range of scales of such flows in the natural environment. Examples of the multitude of engineering free shear flows are the wakes behind moving bodies and the exhausts from jet engines. Most combustion processes and many mixing processes involve turbulent free shear flows. |
Free shear flows in the real world are most often turbulent. Even if generated as laminar flows, they tend to become turbulent much more rapidly than the wall-bounded flows which we will discuss later. This is because the three-dimensional vorticity necessary for the transition to turbulence can develop much more rapidly in the absence of walls that inhibit the qrowth velocity components normal to them. | Free shear flows in the real world are most often turbulent. Even if generated as laminar flows, they tend to become turbulent much more rapidly than the wall-bounded flows which we will discuss later. This is because the three-dimensional vorticity necessary for the transition to turbulence can develop much more rapidly in the absence of walls that inhibit the qrowth velocity components normal to them. | ||
Line 68: | Line 69: | ||
</td><td width="5%">(3)</td></tr></table> | </td><td width="5%">(3)</td></tr></table> | ||
- | In addition, we have the two-dimensional mean | + | In addition, we have the two-dimensional mean continuity equation which reduces to: |
<table width="70%"><tr><td> | <table width="70%"><tr><td> | ||
Line 75: | Line 76: | ||
</math> | </math> | ||
</td><td width="5%">(4)</td></tr></table> | </td><td width="5%">(4)</td></tr></table> | ||
+ | |||
+ | === Order of magnitude estimates === | ||
+ | |||
+ | Now let's make an order of magnitude estimate for each of the terms. This procedure may seem trivial to some, or hopeless hand-waving to others. The fact is that if you fall either of these groups, you have missed something important. Learning to make good order-of-magnitude arguments and knowing when to use them (and when not to use them) are two of the most important skills in fluid mechanics, and especially turbulence. To do this right we will have to be very careful to make sure our estimates accurately characterize the terms we are making them for. | ||
+ | |||
+ | Naturally we should not expect changes of anything in the <math> x </math> - direction to scale the same as changes in the <math> y </math>-direction, especially in view of the above. (This, after all, is the whole idea of "thin" shear flow.) So let's agree that we will pick a length scale, say <math> L </math>, characteristic of changes or mean quantities in the <math> x </math>-direction; i.e. | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial}{\partial x} \sim \frac{1}{L} | ||
+ | </math> | ||
+ | </td><td width="5%">(5)</td></tr></table> | ||
+ | |||
+ | where for now '<math> \sim </math>' means "of order of magnitude of". And we can do the same thing for changes of mean quantities in the <math> y </math>-direction by defining a second length scale, say <math> \delta </math>, to mean: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial}{\partial y } \sim \frac{1}{\delta } | ||
+ | </math> | ||
+ | </td><td width="5%">(6)</td></tr></table> | ||
+ | |||
+ | Note that both these scales will vary with the streamwise position where we evaluate them. A good choice for <math> \delta </math> might be proportional to the local lateral extent of the flow (or its "width"), while <math> L </math> is related to the distance from the source. | ||
+ | |||
+ | Consider the mean velocity in equation 10 in section "Wall-bounded turbulent flows". It occurs in five different terms: <math> U </math> alone; twice with <math> x </math> -derivatives, <math> \partial U/ \partial x </math> and <math> \partial^{2} U / \partial x^2 </math>; and twice with <math> y </math> -derivatives,<math> \partial U/ \partial y </math> and <math> \partial^{2} U / \partial y^2 </math>. Now it would be tempting to simply pick a scale velocity for <math> U </math> say <math> U_s </math>, and use it to estimate all five terms, say as: <math> U_s </math>, <math> U_s/L </math>, <math> U_s/L^2 </math>, <math> U_s/ \delta </math>, and <math> U_s/ \delta^2 </math>. But this is much too naive, and fails to appreciate the true role of the terms we are evaluating. | ||
+ | |||
+ | Look at the examples shown in Figure ??. Our simple approach would provide an appropriate estimate for a jet if we took our velocity scale equal to the mean centerline velocity at a given streamwise location; i.e., <math> U_s \left( x \right) = U_{cl} \left( x \right) </math>. This is because both the mean velocity and the ''changes'' in the mean velocity across the flow are characterized by the centerline velocity. But by contrast, look at the wake. Even relatively close to the wake generator, the wake deficit,<math> \Delta U_{cl} = U_{0} - U_{cl} </math> is small compared to the free stream velocity (i.e. <math> \Delta U_{cl} << U_{0} </math> ). So the obvious choice to scale <math> U \left( x,y \right) </math> is <math> U_{0}</math> On the other hand, an estimate for the velocity gradient across the flow of <math> U_{0}/\delta </math> would be much too big, again because the deficit is so small. Obviously a better choice would be to use the centerline mean velocity ''deficit'' <math> \Delta U_{cl} </math>; i.e., | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial U}{\partial y} \sim \frac{\Delta U_{cl}}{\delta} | ||
+ | </math> | ||
+ | </td><td width="5%">(7)</td></tr></table> | ||
+ | |||
+ | In the order of magnitude analysis below, we shall try to keep the discussion as general as possible by using <math> U_{s} </math> to characterize the mean velocity when it appears by itself, and <math> \Delta U_{s}</math> to represent changes in the mean velocity. For the jet example of the preceding paragraph, both <math> U_{s} </math> and <math> \Delta U_{s} </math> are the same; i.e. <math> U_{s} = U_{cl} </math> and <math> \Delta U_{s} = U_{cl} </math>. But for the wake they are different because of the external stream; i.e., <math> U_{s} = U_{0} </math> and <math> \Delta U_{s} = \Delta U_{cl} </math>. If you can keep in mind ''why'' these differences exist among the various flows, it will be a lot easier to both understand the results and not confuse them. | ||
+ | |||
+ | |||
+ | Now we could distinguish changes of velocity in the <math> x </math> -direction from those in the <math> y </math> - direction. But this level of complexity is not necessary (at least for the examples considered here), especially since we have left the precise defenition of <math> L </math> rather nebulous. What we can do is to use the same estimate for changes in the velocity as for the <math> y </math> - direction, and ''define'' our length scale <math> L </math> to make the estimate based on both correct; i.e. <math> \partial U / \partial x \sim \Delta U_{s} / L </math>. To see why this make sence physically and can be reasoned (as opposed to guessed), let's look at the wake. Pick a spot outside, near the edge of the wake fairly close the generator (say point A). Now proceed at constant <math> y </math> far downstream in <math> x </math>. Eventually the wake will have spread past you and you will be close enough to the centerline so the local mean velocity will be closer to <math> U_{cl} </math> than <math> U_{0} </math>. Obviously we have simply traveled far enough at constant <math> y </math> to ensure that <math> \Delta U_{cl} </math> is the proper scale for ''the changes'' in velocity in the <math> x </math> - direction. If the distance downstream over which this change occured is taken as <math> L </math>, then proper estimate is easily seen to be: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial U}{\partial x} \sim \frac{\Delta U_{cl}}{L} | ||
+ | </math> | ||
+ | </td><td width="5%">(8)</td></tr></table> | ||
+ | |||
+ | But this is exactly what we would have gotten by taking <math> U_{s} = U_{cl} </math> as we agreed above. We simply have absorbed any differences into our choice of <math> L </math>. When considering specific problems and applying similarity techniques, the seemingly arbitrary choices here become quite precise constraints (as we shall see). | ||
+ | |||
+ | We still haven't talked about how to estimate the velocity scale for <math> V </math>, the cross-stream mean velocity component. From the continuity equation we know that: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial V}{\partial y} = - \frac{\partial U}{\partial x} | ||
+ | </math> | ||
+ | </td><td width="5%">(9)</td></tr></table> | ||
+ | |||
+ | From our considerations above, we know that: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial U}{\partial x} \sim \frac{\Delta U_{s}}{L} | ||
+ | </math> | ||
+ | </td><td width="5%">(10)</td></tr></table> | ||
+ | |||
+ | If there is no mean cross flow in the external stream, then the scale for | ||
+ | <math> V </math> is the same as the scale for changes in <math> V </math>. Therefore, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial V}{\partial y} \sim \frac{V_{s}}{\delta} | ||
+ | </math> | ||
+ | </td><td width="5%">(11)</td></tr></table> | ||
+ | |||
+ | It follows immediately that the order of magnitude of the cross-stream velocity is: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | V_{s} \sim \Delta U_{s} \frac{\delta}{L} | ||
+ | </math> | ||
+ | </td><td width="5%">(12)</td></tr></table> | ||
+ | |||
+ | We might have expected something like this if we had thought about it. If the <math> V </math>-velocity were of the same order as the <math> U </math> - velocity, how could the flow in any sense be a "thin shear flow". On the other hand, it also makes sense that <math> V_{s}/U_{s} \propto \delta/L </math> since both are some measure of how the flow spreads/ Note that equation 12 would not be correct estimate for <math> V_{s} </math> if there were an imposed cross-flow, since then we would have to consider <math> V </math> and changes in <math> V </math> separately (exactly as for <math> U </math> ). | ||
+ | |||
+ | The mean pressure gradient term is always a problem to estimate at the outset. Therefore it is better to simply leave this problem alone, and see what is left at the end. In the estimates below you will see a question mark, which simply means we are postponing judgement until we have more information. Sometimes we will have to keep the term simply because we don't know enough to throw it away. Other times it will be obvious that it must remain because there is only one term left that must be balanced by something. | ||
+ | |||
+ | Now we have figured out how to estimate the order of magnitude of all the terms except the turbulence terms. For most problems this turns out to be pretty straight-forward if you remember our discussion of the pressure strain-rate terms. They so effectively distribute the energy that even that the three components of velocity are usually about the same order of magnitude. So if we choose a turbulence scale as simply <math> u </math>, then <math> \left\langle u^{2} \right\rangle \sim u^{2}</math>, <math> \left\langle v^{2} \right\rangle \sim u^{2}</math>, <math> \left\langle w^{2} \right\rangle \sim u^{2}</math>. But what about the Reynolds shear stress components like <math> \left\langle uv \right\rangle </math> ? When acting against the mean shear to produce energy, they tend to be well-correlated and the maximum value of <math> \left\langle uv \right\rangle / u_{rms}v_{rms} \approx 0.4 </math> . Obviously the right choice though since some cross-stress like <math> \left\langle uw \right\rangle </math> are identically zero in plane flows because of homogeniety. Allowing for the possibility of a Reynolds number dependent correlation coefficient will be sent to be especially important when the terms in the turbulence kinetic energy and dissipation equations are considered. | ||
+ | |||
+ | === Streamwise momentum equation === | ||
+ | |||
+ | Let's look now at the <math> x </math>-component of the mean momentum equation and write below each term its order of the magnitude. | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | U \frac{\partial U}{\partial x} & + & V \frac{\partial U}{\partial y} \\ | ||
+ | U_{s} \frac{\Delta U_{s}}{L} & & \left( \Delta U_{s} \frac{\delta}{L} \right) \frac{\Delta U_{s}}{ \delta} \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | = -\frac{1}{\rho} \frac{\partial P}{ \partial x}& - & \frac{ \partial \left\langle U^{2} \right\rangle}{ \partial x} & - & \frac{ \partial \left\langle uv \right\rangle}{ \partial y} & + & \nu \frac{\partial^{2}U}{\partial x^2} & + & \nu \frac{\partial^{2}U}{\partial y^2} \\ | ||
+ | ? & & \frac{u^{2}}{L} & & \frac{u^{2}}{\delta} & & \nu \frac{\Delta U_{s}}{L^{2}}& & \nu \frac{\Delta U_{s}}{\delta^{2}} | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Now the survival of at least one of the terms on the left-hand side of the equation is the essence of free shear flow. Since the largest of these is the first, let's re-scale the others by deviding all the estimates by <math> U_{s} \Delta U_{s} / L </math> . Now we have: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | U \frac{\partial U}{\partial x} & + & V \frac{\partial U}{\partial y} \\ | ||
+ | 1 & & \frac{\Delta U_{s}}{U_{s}} \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | = -\frac{1}{\rho} \frac{\partial P}{ \partial x}& - & \frac{ \partial \left\langle U^{2} \right\rangle}{ \partial x} & - & \frac{ \partial \left\langle uv \right\rangle}{ \partial y} & + & \nu \frac{\partial^{2}U}{\partial x^2} & + & \nu \frac{\partial^{2}U}{\partial y^2} \\ | ||
+ | ? & & \frac{u^{2}}{U_{s} \Delta U_{s}} & & \frac{u^{2}}{U_{s} \Delta U_{s}} \left( \frac{L}{\delta} \right) & & \frac{\nu}{U_{s} L} & & \frac{\nu}{U_{s} L} \left( \frac{L}{\delta} \right) | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Just in case you need to be reminded, this is a turbulence course. Therefore almost all the interesting flows are at high Reynolds number. But what does this mean? Now we can say. It could mean the viscous terms in the mean momentum equation are negligible, which in tuurn means that: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{U_{s}L}{\nu} >> 1 | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | and | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{U_{s} \delta}{\nu} >> \frac{L}{\delta} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Obviously the second criterion is much more stringent, since <math> L / \delta </math> is typically of order 10. When <math> U_{s} \delta / \nu > 1000 </math>, the contributions of the viscous stresses are certainly negligible, at least as far as the <math> x </math>-component of the mean mometum equation is concerned. But if <math> U_{s} \delta / \nu \sim 100 | ||
+ | </math> only, this is a pretty marginal assumption, and you might want to retain the last term in your analysis. Such is unfortunately the case in many experiments that are often at quite low Reynolds numbers. | ||
+ | |||
+ | So, what then, tou ask, do we do about the turbulence terms? To repeat this is a TURBULENCE course! Which means there is no way we are going to throw away all the turbulence terms on the right-hand side of the equation. There is no magic or hocus-pocus about this, you simply can't have turbulence without turbulence terms. And if there is no turbulence, we are really aren't too interested. So, what does this tell us? It tells us about <math>\delta </math>! Surprised? I bet you thought it would tell us about <math> u </math>, right? Nope. Look at the right-hand side of the second equation of this subsection and the orders of magnitude below it. Almost always <math> u < U_{s} </math>; sometimes much less and certainly no bigger. Then the biggest turbulence term is one involving the Reynolds shear stress <math> \partial \left\langle -uv \right\rangle / \partial y </math> which we have estimated as <math> \left[ u^{2} / \left( U_{s} \Delta U_{s} \right) \right] \left( L / \delta \right) </math>. Hence there can be no turbulence terms at all unless: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\delta}{L} \sim \frac{u^{2}}{U_{s} \Delta U_{s}} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Wow! Look what we have learned about turbulence without ever solving a single equation. We know how the growth of our free shear flow relates to the turbulence intensity - at least in terms of order of magnitude. The simlarity theories described below will be able to actually deduce the <math> x </math> - ''dependence'' of <math>\delta </math>, again without solving the equations. | ||
+ | |||
+ | Does the argument above mean the other turbulence term is negligible since it is only of order <math> u^{2} / U_{s} \Delta U_{s} </math>? Well, no question it is usually small, typically less than 5%. But unlike the viscous terms, this term often does not get smaller the father we go in the streamwise direction. So the answer depends on the question we are asking and how accurately we want our answer to be. If we are asking only first order questions and are willing to accept errors in our solutions (or experiments) to within a 5-10 percent, then this normal stress is certainly negligible. Before you are ready to forget about this, consider the following: since <math> \delta / L < 0.1 </math> typically, ''to first order free shear flows do not grow at all!'' So if you want to talk about whether <math> d \delta / dx </math> for a particular shear flow is 0.095 instead of 0.09 (turbulence wars have been waged over less), then you better be dealing with a full desk - meaning these second order terms had better be retained in any calculation or experimental verification (like momentum or energy balances). | ||
+ | |||
+ | Much of the confusion in the literature about free shear flows comes from the failure to use second order accuracy to make second order statements. And even more unbelievably, most attempts to do so-called code validation of second-order turbulence models are based on mesurements which are only first order accurate. Now I know you think I am kidding - but try to do a simple momentum balance on the data to within order <math> u^{2} / U_{s} \Delta U_{s} </math>. Did I hear you say it was ridiculous to use second-order equations (like the Reynolds stress equations or even the kinetic energyequation) and determine constants to three decimal places using first order accurate data? Shhhh!!!! Your collegues might be annoyed unless you keeep this a secret. | ||
+ | |||
+ | Hopefully you have followed all this. If not, go forward and "faith" will come to you. I'll summarize how we do this: | ||
+ | * You have to look very carefully at the particular flow you wish to analyze. | ||
+ | * Then, based on physical reasoning and data (if you have it), you make estimates you think are appropriate. | ||
+ | * Then you use these estimates to decide which terms you have to keep, and which you can ignore. | ||
+ | * And when you are completely finished, you carefully look back to make sure that all the terms you kept really are important and remain important as the flow develops. | ||
+ | |||
+ | Finally, if you are really good at this, you look carefully to make sure that all the terms you kept do not all vanish at the same place (by going from positive to negative, for example). If this happens, then the terms you neglected may in fact be the only terms left in the equation thet are not exactly zero - oops! This can make fo some very interesting situations/ (The famous "critical layer" of boundary layer stability theory is an example of this.) | ||
+ | |||
+ | === The transverse momentum equation === | ||
+ | |||
+ | Now we must also cosider the transverse or <math> y </math>-momentum equation. The most important thing we have to remember is that the <math> y </math>- equation can NOT be considered indepedently from the <math> x </math>-momentum equation. We would never consider trying to simplify a vector like <math> (a,b,c) </math> by dividing only one component by <math> S </math> say to produce <math> (Sa,b,c) </math>, since this would destroy the whole idea of a vector as having a direction. Similarly we must not do this for a vector equation. We have already decided thet the first term on the left-hand side of the <math> x </math>-momentum equation was the term we had to keep, and we divided by its order of magnitude, <math> U_{s} \Delta U_{s} / L </math> , to make sure it was of order one. Thus we have already decided how we are going to scale all the components of the vector equation, and so we must do exactly the same thing here. but first we must estimate the order of magnitude of each term, exactly as before. Using our previous results, here's what we get: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | U\frac{\partial V}{ \partial x} & + & V \frac{\partial V}{\partial y} \\ | ||
+ | U_{s} \frac{\Delta U_{s} \delta / L}{L} & &\left(\Delta U_{s} \frac{\delta}{L} \right) \frac{\Delta U_{s} \delta / L}{\delta} \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | = - & \frac{1}{\rho} \frac{\partial P}{ \partial y} & - & \frac{\partial \left\langle uv \right\rangle}{ \partial x} & - & \frac{ \partial \left\langle v^{2} \right\rangle}{ \partial y} & + & \nu \frac{\partial^{2} V}{ \partial x^{2}} & + & \nu \frac{\partial^{2} V}{ \partial y^{2}} \\ | ||
+ | & ? & & \frac{ u^{2}}{L}& & \frac{ u^{2}}{\delta} & &\nu \frac{\Delta U_{s} \delta / L}{L^{2}} & & \nu \frac{\Delta U_{s} \delta / L}{\delta^{2}} \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Now dividing each term by <math> U_{s} \Delta U_{s} / L </math>, exactly as before, we obtain: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | U\frac{\partial V}{ \partial x} & + & V \frac{\partial V}{\partial y} \\ | ||
+ | \frac{\delta}{L} & & \frac{\Delta U_{s}}{U_{s}} \left( \frac{\delta}{L} \right) \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \begin{matrix} | ||
+ | = - & \frac{1}{\rho} \frac{\partial P}{ \partial y} & - & \frac{\partial \left\langle uv \right\rangle}{ \partial x} & - & \frac{ \partial \left\langle v^{2} \right\rangle}{ \partial y} & + & \nu \frac{\partial^{2} V}{ \partial x^{2}} & + & \nu \frac{\partial^{2} V}{ \partial y^{2}} \\ | ||
+ | & ? & & \frac{u^{2}}{U_{s} \Delta U_{s} } & & \frac{u^{2}}{U_{s} \Delta U_{s} } \left( \frac{L}{\delta} \right) & & \frac{\nu}{U_{s} L} \left( \frac{\delta}{L}\right) & & \frac{\nu}{U_{s} \delta} \\ | ||
+ | \end{matrix} | ||
+ | </math> | ||
+ | </td><td width="5%"></td></tr></table> | ||
+ | |||
+ | Unless you have seen this all before, the left-hand side is probably a surprise: none of the mean convection terms are of order one! In fact the only estimated term that is of order one in the whole equation is <math> \partial \left\langle v^{2} \right\rangle / \partial y </math> and only because we have already agreed that <math> \left( u^{2} / U_{s} \Delta U_{s} \right) \left( L / \delta \right) </math> had to be of order one to keep a turbulence term in the <math> x </math>-momentum equation. Of course, there cannot be an equation with only a single term equal to zero, unless we have badly over-estimated its order of magnitude. Fortunately there is the pressure term left to balance it, so to first order in <math> \delta / L \sim u^{2} / U_{s} \Delta U_{s} </math>, the <math> y </math>-momentum equation reduces to simply: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | 0 \approx - \frac{1}{\rho} \frac{\partial P}{\partial y} - \frac{\partial \left\langle v^{2} \right\rangle}{\partial y} | ||
+ | </math> | ||
+ | </td><td width="5%">(13)</td></tr></table> | ||
+ | |||
+ | This simple equation has an equally simple interpretation. It says thet the change in the mean pressure is only due to the radisl gradient of the transverse component of the Reynolds normal stress, <math> \left\langle v^{2} \right\rangle </math>. | ||
+ | |||
+ | We can integrate equation 13 ''across'' the shear layer from a given value of <math> y </math> to infinity (or anywhere else for thet matter) to obtain: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | P \left( x, y \right) = P \left( x, \infty \right) - \rho \left\langle v^{2} \right\rangle | ||
+ | </math> | ||
+ | </td><td width="5%">(14)</td></tr></table> | ||
+ | |||
+ | If the free stream is at constant (or zero) mean velocity, then <math> P \left( x, \infty \right) = P_{\infty} = constant </math>, so we can write simply: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | P \left( x, y \right) = P_{\infty} - \rho \left\langle v^{2} \right\rangle | ||
+ | </math> | ||
+ | </td><td width="5%">(15)</td></tr></table> | ||
+ | |||
+ | Either of these can be substituted into <math> x </math>-momentum equation to eliminate the pressure entirely, as we shall show below. | ||
+ | |||
+ | Now since <math> \left\langle v^{2} \right\rangle << \Delta U^{2}_{s} </math> typically <math> P \approx P_{\infty} </math>, at least to first order in <math> u^{2} / \Delta U^{2}_{s} </math>. Therefore one might be tempted to conclude that the small pressure changes across the flow are not important. And they are not, of course, if only the streamwise momentun is considered. But without this small mean pressure gradient ''across'' the flow, there would be no entrainment and no growth of the shear layer. In other words, the flow would remain parallel. clearly whether an effect is negligible or not depends on which question is being asked. | ||
+ | |||
+ | === The free shear layer equations === | ||
+ | |||
+ | If we assume the Reynolds number is always large enough that the viscous terms can be neglected, then as noted above, the pressure term in the terms of <math> P_{\infty} </math> and <math> \left\langle v^{2} \right\rangle </math>. Thus, ''to second-order in'' <math> u^{2} / U_{s} \Delta U_{s} </math> ''or'' <math> \delta / L </math>, the momentum equations for a free shear flow reduce to a single equation: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U \frac{\partial U}{ \partial x} + \left\{ V \frac{\partial U}{ \partial y} \right\} = - \frac{d P_{\infty}}{dx} - \frac{\partial }{ \partial y} \left\langle uv \right\rangle - \left\{ \frac{\partial }{ \partial x} \left[ \left\langle u^{2} \right\rangle - \left\langle v^{2} \right\rangle \right] \right\} | ||
+ | </math> | ||
+ | </td><td width="5%">(14)</td></tr></table> | ||
+ | |||
+ | As we have seen above, the second term on the left-hand side may or may not be important, depending on whether <math> U_{s} = \Delta U_{s} </math> or <math> \Delta U_{s}/ U_{s} << 1 </math>. Clearly this depends on the velocity deficit (or excess) relative to the free stream. The term in brackets on the right-hand side is second-order (in <math> u^{2} / U_{s} \Delta U_{s} \sim \delta / L </math> ) compared to the others, and so could rigthgly be neglected. It has been retained for now since in some flows it does not vanish with increasing distance from the source. Moreover, it can be quite important when considering integrals of the momentum equation, since the profiles of <math> \left\langle u^{2} \right\rangle </math> and <math> \left\langle v^{2} \right\rangle </math> do not vanish as rapidly with increasing <math> y </math> as do the other terms. | ||
+ | |||
+ | == Two-dimensional Turbulent Jets == | ||
+ | |||
+ | Turbulent jets are generated by a concentrated source of momentum issuing into an ambient environment. The state of the environment and the external boundary conditions are very important in determining how the turbulent flow evolves. The simplest case to consider is that in which the environment is at rest and unbounded, so all the boundary conditions are homogeneous. The jet itself can be assumed to isssue from a line source or a slot as shown in figure '''??'''. The jet then evolves and spreads downstream by entraining mass from the surroundings which are at most in irrotational motion induced by the vortical fluid in the jet. But no new momentum is added to the flow, and it is this fact that distinguishes the jet from all other flows. As we shall see below the rate at which momentum crosses any <math> x </math> - plane is not quite constant at the source value, however, due to the small streamwise pressure gradient arising from the turbulence normal stresses. | ||
+ | |||
+ | The avereged continuity equation can be integrated from the centerline (where it is zero by symmetry) to obtain <math> V </math>, i.e., | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | V = - \int^{y}_{0} \frac{\partial U}{ \partial x} d \tilde{y} | ||
+ | </math> | ||
+ | </td><td width="5%">(17)</td></tr></table> | ||
+ | |||
+ | It immediately follows that the <math> V </math> velocity at <math> + / - \infty </math> is given by: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | V_{\infty} = - V_{-\infty} = - \frac{ d }{ d x} \int^{\infty}_{0} U \left( x , y \right) dy | ||
+ | </math> | ||
+ | </td><td width="5%">(18)</td></tr></table> | ||
+ | |||
+ | Twice the integral represents the total volume flow rate crossing a given <math> x </math> - plane. Therefore it makes sense that the rate of increase of this integral is the enrainment velocity, <math> V_{\infty} </math>. Obviously, <math> V_{\infty} </math> cannot be zero if the jet is to spread, at least in two-dimensional flow. It is easy to see that the entrainment velocity calculated from the integral is consistent with our order of magnitude estimate above (i.e., <math> V \sim \Delta U_{s} \delta / L </math> ). Note that the integral makes no sense if the mean velocity <math> U </math> does ''not'' go to zero with increasing <math> \left| y \right| </math>. | ||
+ | |||
+ | If the free stream is assumed to have zero streamwise velocity, then the scales for the velocity and gradienits of it are the same, so the equations reduce to simply equation 16 in '''the next section'''(need to edit) and the mean continuity equation. The latter can be multiplied by <math> U </math> to yield: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U\left\{ \frac{\partial U}{\partial x} + \frac{\partial V}{ \partial y} \right\} = 0 | ||
+ | </math> | ||
+ | </td><td width="5%">(19)</td></tr></table> | ||
+ | |||
+ | When this is added to equation 16 from '''the next section'''(need to edit) the terms can be combined to obtain: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{\partial U^{2}}{\partial x} + \frac{\partial UV}{\partial y} + \frac{\partial \left\langle uv \right\rangle}{ \partial y} + \frac{ \partial }{ \partial x} \left[ \left\langle u^{2} \right\rangle - \left\langle v^{2} \right\rangle \right] = 0 | ||
+ | </math> | ||
+ | </td><td width="5%">(20)</td></tr></table> | ||
+ | |||
+ | This can be integrated across the flow for any given value of <math> x </math> to obtain: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{d}{dx} \int^{\infty}_{-\infty} \left[ U^{2} + \left( \left\langle u^{2} \right\rangle - \left\langle v^{2} \right\rangle \right) \right] dy = 0 | ||
+ | </math> | ||
+ | </td><td width="5%">(21)</td></tr></table> | ||
+ | |||
+ | where we have ''assumed'' that <math> U </math>, <math> \left\langle u^{2} \right\rangle </math>, and <math> \left\langle v^{2} \right\rangle </math> vanish as <math> \left|y \right| \rightarrow \infty </math>. (Remember these assumptions the next time someone tries to tell you that a small co-flowing stream and modest background turbulence level in the external flow are not important.) | ||
+ | |||
+ | Equation (21) can in turn be integrated from the source, say <math> x=0 </math>, to obtain: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | M_{o} = \int^{\infty}_{-\infty} \left[ U^{2} + \left( \left\langle u^{2} \right\rangle - \left\langle v^{2} \right\rangle \right) \right] dy | ||
+ | </math> | ||
+ | </td><td width="5%">(22)</td></tr></table> | ||
+ | |||
+ | where <math> M_{o} </math> is the rate at which momentum per unit mass per unit length is added at the source. The first two terms of the itegral are the flux of streamwise momentum due to the mean flow and the turbulence. The last term is due to the streamwise pressure gradient obtained by integrating the <math> y </math>-momentum equation. | ||
+ | |||
+ | This is a very powerful constraint, and together with the reduced governing | ||
+ | equations allows determination of a similarity solution for this problem. Equation 22 also is of great value to experimentalists, since it can be used to confirm whether their measurements are valid and whether the flow indeed approximates a plane jet in an infinite environment at rest. Be fore-warned, most do not. And the explanations as to why range from hilarious to pathetic. Consider this, if you measure the velocity moments to within say 1%, then shouldn’t you be able to estimate the integral even more accurately, since the random errors from the integration should be reduced by one over the square root of the number of points you used in the integration? Be very wary of measurements which have not been tested against the momentum integral — for this and all flows. There is usually | ||
+ | a reason the authors fail to mention it. | ||
+ | |||
+ | === Similarity analysis of the plane jet === | ||
+ | |||
+ | '''Foreword''' The results of this section are VERY MUCH at odds with those in all texts, including Pope 2000. There is an increasing body of literature to support the position taken, however, especially the multiplicity of solutions and the role of upstream or initial conditions. Since these ideas have a profound effect on how we think about turbulence, it is important that you consider both the old and the new, and be aware of why the conclusions are different. | ||
+ | |||
+ | The message you should take from the previous section is that regardless of how the plane jet begins, the momentum ''integral'' should be constant at <math> M_{o} </math>. Note that Schneider (1985) has shown that the nature of the entrained flow in the neighborhood of the source can modify this value, especially for plane jets. Also the term <math> \left\langle u^{2} \right\rangle - \left\langle v^{2} \right\rangle </math>, while important for accounting for the momentum balance in the analysis below with no loss in generality. | ||
+ | |||
+ | The search for similarity solutions to the averaged equations of motion is exactly like the approach utilized in laminar flows (eg. Batchelor 1967) except that here there are more unknows than equations, i.e. the averaged equations are not closed. Thus solutions are sought which are of the form | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U = U_{s} \left( x \right) f \left( \overline{y} \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(23)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | - \left\langle uv \right\rangle = R_{s} \left( x \right) g \left( \overline{y} \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(24)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | missed | ||
+ | </math> | ||
+ | </td><td width="5%">(25)</td></tr></table> | ||
+ | |||
+ | where | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \overline{y} = y / \delta \left( x \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(26)</td></tr></table> | ||
+ | |||
+ | The scale functions are functions of <math> x </math> only, i.e. <math> U_{s} = U_{s} \left( x \right) \ R_{s} = R_{s} \left( x \right) </math> and <math> \delta = \delta \left( x \right) </math> only. Note that the point of departure from the traditional turbulence analyses was passed when it was not arbitrarily decided that <math> R_{s} = {U}^{2}_{s} </math>, etc. (c.f. Tennekes and Lumley 1972, Townsend 1976, Pope 2000). | ||
+ | |||
+ | Now let you think there is something weird about the way we have written the form of the solutions we say we are looking for, consider this. Suppose you were looking for some <math> x </math>- dependent quantities you could normalize the profiles by to make them all collepse together - like any experimentalist almost always does. Such a search is called the search for "''scales''" , meaning "''scales''" which collapse the data. Now precisely what does "collapse the data" really mean? Exactly this: "collapse the data" means that the form of the solution must be an <math> x </math>-independent solution to the governing equations. But this is exactly what we are looking for in equations 23 and 24. Obviously we need to plug them into the averaged equations and see if such solution is possible. We already know which terms in these equations are important. What we want to do now is to see if these govering equations admit to solutions for which these terms ''remain'' important throughout the jet development. These are called ''equilibrium similarity solutions''. Note that for any other form of solution, one or more terms could die off, so that others dominate. | ||
+ | |||
+ | Differentiating equations 23 and 24 using the chain-rule, substituting into the averaged momentum equation, and clearing the terms yields: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \left[ \frac{\delta}{U_{s}} \frac{ dU_{s}}{dx} \right] f^{2} - \left[ \frac{ \delta }{U_{s}} \frac{ dU_{s}}{dx} + \frac{d \delta}{dx} \right] f' \int^{\overline{y}}_{0} f\left( \overline{y}' \right) d \overline{y}' = \left[ \frac{R_{s}}{U^{2}_{s}} \right] g' | ||
+ | </math> | ||
+ | </td><td width="5%">(27)</td></tr></table> | ||
+ | |||
+ | Now, all of the <math> x </math>-dependence is in the square bracketed terms, so equilibrium similarity solutions are possible only if all the bracketed terms have the same <math> x </math>-dependence, i.e., | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{d \delta}{d x} \propto \frac{ \delta}{ U_{s}} \frac{ d U_{s}}{ dx } \propto \frac{ R_{s}}{ U_{s}} | ||
+ | </math> | ||
+ | </td><td width="5%">(28)</td></tr></table> | ||
+ | |||
+ | The coefficients of proportionality, say <math> n </math> and <math> B </math> respectively, can at most depend on the details of how the flow began; i.e., the mistery argument "*". | ||
+ | |||
+ | Substitution into the momentum integral yields to first order, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | A U^{2}_{s} \delta = M_{o} | ||
+ | </math> | ||
+ | </td><td width="5%">(29)</td></tr></table> | ||
+ | |||
+ | where | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | A \equiv \int^{\infty}_{0} f^{2} \left( \overline{y} \right) d \overline{y} | ||
+ | </math> | ||
+ | </td><td width="5%">(30)</td></tr></table> | ||
+ | |||
+ | <math> A </math> and <math> B </math> are parameters which also can at most depend on how the flow began (since every other dependence has been eliminated). Moreover, since <math> M_{o} </math> is independent of <math> x </math>, then <math> n=1/2 </math> and it follows that: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U_{s} \propto M^{1/2}_{o} \delta^{1/2} | ||
+ | </math> | ||
+ | </td><td width="5%">(31)</td></tr></table> | ||
+ | |||
+ | It's easy to see that all of the relations involving <math> U_{s} </math> are proportional to <math> d \delta / dx </math>, leaving only | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | - \frac{1}{2} \left\{ f^{2} - f' \int^{\overline{y}}_{0} f \left( \overline{y}' \right) d \overline{y}' \right\} = \left[ \frac{R_{s}}{ U^{2}_{s} d \delta / d x} \right] g' | ||
+ | </math> | ||
+ | </td><td width="5%">(32)</td></tr></table> | ||
+ | |||
+ | Thus the only remaining necessary condition for similarity is | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | R_{s} \propto U^{2}_{s} \frac{d \delta}{ d x } | ||
+ | </math> | ||
+ | </td><td width="5%">(33)</td></tr></table> | ||
+ | |||
+ | Note that, unlike the similarity solutions encountered in laminar flows, it's possible to have a jet which is similar ''without'' having some form of power law behavior. In fact, the <math> x </math>-dependence of the flow may not be known at all because of the closure problem. Nonetheless, the profiles will collapse with the local length scale. In fact, it is easy to show that the profiles for all plane jets will be alike, if normalized properly, even if the growth rates are quite different. The scale velocity <math> U_{s} </math> can be defined to be the centerline velocity <math> U_{cl} </math> by absorbing an appropriate factor into the profile function <math> f \left( \overline{y} \right) </math>. Also, the entire factor in brackets on the right-hand side of equation 32 can be absorbed into the Reynolds stress profile, <math> g \left( \overline{y} \right) </math>. Finally, if the length scale is chosen the same for all flows under consideration (e.g., the half-width, say <math> \delta_{1/2} </math>, defined as the distance between the center and where <math> U = U_{cl}/2 </math>), then the similarity equation govering all jets reduces to, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | f^{2} + f' \int^{\overline{y}}_{0} f \left( \overline{y}' \right) d \overline{y}' = - 2 \overline{g} ' | ||
+ | </math> | ||
+ | </td><td width="5%">(34)</td></tr></table> | ||
+ | |||
+ | where both <math> f </math> has been modified so that <math> f(0)=1 </math> and <math> \overline{g} \left( \overline{y} \right) </math> is just equal to <math> g </math> with the factor <math> R_{s} / \left( U^{2}_{s} d \delta / dx \right) </math> absorbed into it. Thus all jets, regardless of initial conditions, will produce the same similarity velocity profile, when scaled with <math> U_{cl} </math> and <math> \delta_{1/2} </math>. This will be not true for the Reynolds stress, however, unless the coefficient of proportionality between <math> R_{s} </math> and <math> U^{2}_{s} d \delta / dx </math> is the same for ''all'' plane jets. In general, it is not; so here is where the dependence on source conditions shows up. | ||
+ | |||
+ | It is immediately obvious from equation 32 that usual arbitrary choice of <math>R_{s}={U}^{2}_{s}</math>considerably restricts the possible solutions to those for which <math> d \delta / dx = constant </math>, or plane jets which grow linearly. However, eben if the growth were linear under some conditions, there is nothing in the theory to this point which demands that the constant be universal and independent of how the jet began. | ||
+ | |||
+ | The idea that jets might all be described be the same asymptotic growth rate stems from the idea of a jet formed by a point source of momentum only, say <math> \rho M_{o} </math>. Such a source must be infinitesimal in size, since any finite size will also provide a source of mass flow. In the absence of scaling parameters other than distance from the source, <math> x </math>, the only possibilities are <math> U_{s} \propto \left( M_{o} / x \right)^{1 / 2} </math>, <math> R_{s} \propto M_{o} / x </math>, <math> \delta \propto x </math>. Obviously <math> R_{s}= {U}^{2}_{s} </math>, and the constants of proportionality can be assumed "universal" since there is only one way to create a point source jet, and this is also contarary to the usual ''assumption''. | ||
+ | |||
+ | The problem with the point source idea is not in the idea itself, but that it has been accepted for so long by so many without justification as the asymptotic state of finite source jets. The usual line of argument (eg. Monin and Yaglom 1972) if one is presented at all is as follows. The jet entrains mass and the rate at which this occurs increases with distance downstream, so that the rate of mass entrainment ultimately overwhelms that added at the source. It is then inferred that the latter can be neglected. And what is the experimental proof for all this? Precisely that the profiles of <math> U/U_{cl} </math> collapse to a single profile for all types of initial conditions. But as was demonstrated above, the mean momentum equation tells us that the mean velocity profiles will collapse to a single profile ''no matter what''! If there is a dependence on the initial conditions it can only show up in the other turbulence moments, in profiles of the Reynolds shear stress normalized by <math> {U}^{2}_{cl} </math> (instead of <math> d \delta / dx </math>), and, of course, in <math> d \delta / dx </math> itself. And guess what quantities no two experiments ever seem to agree on? — All of the above! Could all those experiments plus this nice theory with virtually no assumptions be wrong? Many think so. It is truly amazing how much good science people will throw away in order to hang-on to a couple of ideas that were never derived in the first place but assumed; namely, <math> R_{s} = {U}^{2}_{s} </math> and the asymptotic independence of initial conditions. But then there are still many people who believe the world is flat. | ||
+ | |||
+ | George (1989) demonstrated for an axisymmetric jet that even simple dimensional analysis suggests that the classical theory is wrong, and the same arguments apply here. Suppose that in addition to momentum, mass is added at the source at a rate of <math> \rho m_{o} </math>. Now there is an additional parameter to be considered, and as a consequence, an additional length scale given by <math> L = m^{2}_{0} / M_{o} </math> . Thus the most that can be inferred from dimensional analysis is that <math> \delta / x </math>, <math> U_{s} x^{1/2} / M^{1/2}_{o} </math> and <math> R_{s} x / M_{o} </math> are functions of <math> x/L </math>, with no insight offered into what the functions might be. | ||
+ | |||
+ | ===Implications of the Reynolds stress equations=== | ||
+ | |||
+ | George 1989 further argued that some insight into the growth rate of the jet can be obtained by considering the conditions for similarity solutions of the higher moment equations, in particular the kinetic energy equation. The choice of the kinetic energy equation for further analysis was unfortunate since it implicitly assumed that the three components of kinetic energy all scaled in the same manner. This is, in fact, true only if <math> d \delta / dx = constant </math>, which is certainly not true ''a priory'' for turbulent shear flows. Therefore, here the individual component equations of the Reynolds stress will be considered (as should have been done in George 1989). | ||
+ | |||
+ | For the plane jet the equation for <math> \left\langle u^{2} \right\rangle </math> can be written to first oder (Tennekes and Lumley 1972) as | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U \frac{\partial \left\langle u^{2} \right\rangle}{ \partial x} + V \frac{ \partial \left\langle u^{2} \right\rangle}{ \partial y} = | ||
+ | </math> | ||
+ | </td><td width="5%">(35)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | = \frac{2}{\rho} \left\langle p \frac{\partial u}{\partial x} \right\rangle + \frac{\partial}{\partial y} \left\{ -\left\langle u^{2}v \right\rangle \right\} - 2 \left\langle uv \right\rangle \frac{\partial U}{\partial y} - 2 \epsilon_{u} | ||
+ | </math> | ||
+ | </td><td width="5%">(35)</td></tr></table> | ||
+ | |||
+ | where <math> \epsilon_{u} </math> is the energy dissipation rate for <math> \left\langle u^{2} \right\rangle </math>. | ||
+ | |||
+ | By considering similarity forms for the new moments like | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \ \frac{1}{2} \left\langle u^{2}\right\rangle =K_{u}\left(x\right)k\left(\eta\right) | ||
+ | </math> | ||
+ | </td><td width="5%">(36)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \left\langle p \frac{\partial u}{\partial x} \right\rangle = P_{u} \left( x \right) p_{u} \left( \eta \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(37)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | -\frac{1}{2} \left\langle u^{2} v \right\rangle = T_{u^{2}v} \left(x \right) {t} \left( \eta \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(38)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \epsilon_{u} = D_{u} \left( x \right) {d} \left( \epsilon \right) | ||
+ | </math> | ||
+ | </td><td width="5%">(39)</td></tr></table> | ||
+ | |||
+ | and using <math> R_{s} = U^{2}_{s} d \delta / dx </math>, it's easy to show that similarity of the <math> \left\langle u^{2} \right\rangle </math>-equation is possible only if | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | K_{u} \propto U^{2}_{s} | ||
+ | </math> | ||
+ | </td><td width="5%">(40)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | P_{u} \propto \frac{U^{3}_{s}}{\delta} \frac{d \delta }{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(41)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | T_{u^{2}v} \propto U^{3}_{s} \frac{d \delta }{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(42)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | D_{u} \propto \frac{U^{3}_{s}}{\delta} \frac{d \delta }{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(43)</td></tr></table> | ||
+ | |||
+ | All of these are somewhat surprising: The first (even though a second moment like the Reynolds stress) because the factor of <math> d \delta / dx </math> is absent; the second, third and forth because it is present. | ||
+ | |||
+ | Similar equations can be written for the <math>\left\langle {v}^2 \right\rangle </math> ,<math>\left\langle w^2 \right\rangle </math>, and <math>\left\langle -uv \right\rangle </math>-equations; i.e. | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U \frac{\partial \left\langle v^{2} \right\rangle}{ \partial x} + V \frac{\partial \left\langle v^{2} \right\rangle}{ \partial y} = | ||
+ | </math> | ||
+ | </td><td width="5%">(44)</td></tr></table> | ||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | = 2 \left\langle p \frac{ \partial v}{ \partial y} \right\rangle + \frac{ \partial }{ \partial y} \left\{ - \left\langle v^{3} \right\rangle -2 \left\langle pv \right\rangle \right\} - 2 \epsilon_{v} | ||
+ | </math> | ||
+ | </td><td width="5%">(44)</td></tr></table> | ||
+ | |||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U \frac{\partial \left\langle w^{2} \right\rangle}{ \partial x} + V \frac{\partial \left\langle w^{2} \right\rangle}{ \partial y} = | ||
+ | </math> | ||
+ | </td><td width="5%">(45)</td></tr></table> | ||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | = 2 \left\langle p \frac{ \partial w}{ \partial z} \right\rangle + \frac{ \partial }{ \partial y} \left\{ - \left\langle w^{2}v \right\rangle \right\} - 2 \epsilon_{w} | ||
+ | </math> | ||
+ | </td><td width="5%">(45)</td></tr></table> | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | U \frac{\partial \left\langle uv \right\rangle}{ \partial x} + V \frac{\partial \left\langle uv \right\rangle}{ \partial y} = | ||
+ | </math> | ||
+ | </td><td width="5%">(46)</td></tr></table> | ||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | = \left\langle p \left( \frac{ \partial u}{ \partial y} + \frac{ \partial v}{ \partial x} \right) \right\rangle + \frac{ \partial }{ \partial y} \left\{ - \left\langle u v^{2} \right\rangle \right\} - \left\langle v^{2} \right\rangle \frac{ \partial U}{ \partial y} | ||
+ | </math> | ||
+ | </td><td width="5%">(46)</td></tr></table> | ||
+ | |||
+ | When each of the terms in these equations is expressed in similarity variables, the resulting similarity conditions are: | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \ D_{v} \propto P_{v} \propto \frac{U_{s} K_{v}}{ \delta} \frac{ d \delta}{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(47)</td></tr></table> | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | D_{w} \propto P_{w} \propto \frac{U_{s} K_{w}}{ \delta} \frac{ d \delta}{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(48)</td></tr></table> | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | T_{v^{3}} \propto \frac{U_{s} K_{v}}{\delta} \frac{ d \delta}{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(49)</td></tr></table> | ||
+ | |||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | T_{w^{2}v} \propto \frac{U_{s} K_{w}}{\delta} \frac{ d \delta}{dx} | ||
+ | </math> | ||
+ | </td><td width="5%">(50)</td></tr></table> | ||
+ | |||
+ | and the real surprise, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | K_{v} \propto U^{2}_{s} \left( \frac{ d \delta}{dx} \right)^{2} | ||
+ | </math> | ||
+ | </td><td width="5%">(51)</td></tr></table> | ||
+ | |||
+ | There is an additional equation which must be accounted for: namely that the sum of the pressure strain-rate terms in the component energy equations be zero (from ontinuity). Thus, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \left\langle p \frac{ \partial u}{ \partial x} \right\rangle + \left\langle p \frac{ \partial v}{ \partial y} \right\rangle + \left\langle p \frac{ \partial w}{ \partial z} \right\rangle =0 | ||
+ | </math> | ||
+ | </td><td width="5%">(52)</td></tr></table> | ||
+ | |||
+ | or in similarity variables, | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | P_{u} \left( x \right) p_{u} \left( \eta \right) + P_{v} \left( x \right) p_{v} \left( \eta \right) + P_{w} \left( x \right) p_{w} \left( \eta \right) = 0 | ||
+ | </math> | ||
+ | </td><td width="5%">(53)</td></tr></table> | ||
+ | |||
+ | This can be true for all <math> \eta </math> only if | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | P_{u} \propto P_{v} \propto P_{w} | ||
+ | </math> | ||
+ | </td><td width="5%">(54)</td></tr></table> | ||
+ | |||
+ | An immediate consequence is that | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | D_{u} \propto D_{v} \propto D_{w} | ||
+ | </math> | ||
+ | </td><td width="5%">(55)</td></tr></table> | ||
+ | |||
+ | From equations 41,47 and 48 it also follows that the constraint imposed by 54 can be satisfied only if | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | K_{u} \propto K_{v} \propto K_{w} | ||
+ | </math> | ||
+ | </td><td width="5%">(51)</td></tr></table> | ||
+ | |||
+ | But from equation 51, this can be true only if | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | \frac{d \delta}{dx} = constant | ||
+ | </math> | ||
+ | </td><td width="5%">(57)</td></tr></table> | ||
+ | |||
+ | The relations given by equations 56 and 55 were assumed without proof in the George 1989 analysis. The additional constraint imposed by equation 57 was | ||
+ | not derived, however, and arises from the additional information provided by the pressure strain-rate terms. | ||
+ | |||
+ | Hence, similarity solutions of the Reynolds stress equations are possible only if | ||
+ | |||
+ | <table width="70%"><tr><td> | ||
+ | :<math> | ||
+ | D_{s} \left( x \right) \propto \frac{U^{3}_{s}}{\delta} | ||
+ | </math> | ||
+ | </td><td width="5%">(58)</td></tr></table> | ||
+ | |||
+ | It is an immediate consequence of the earlier discussion on the nature of the dissipation that there are only two possibilities for this to occur: | ||
+ | * Either the local Reynolds number of the flow is constant so that the effect of | ||
+ | the dissipation on the energy containing eddies (and those producing the Reynolds stress as well) does not vary with downstream distance; or | ||
+ | |||
+ | * The local turbulence Reynolds number is high enough so that the relation <math>\epsilon \propto q^3 / L </math> is approximately valid (for a physical length <math>L \propto \delta</math>). | ||
+ | |||
+ | |||
+ | Unlike some flows (like the axisymmetric jet or plane wake) where the local Reynolds number is constant, for the plane jet it varies with downstream distance. Therefore the only possibility for similarity at the level of the Reynolds stresses is (''ii''). This can occur only when the turbulence Reynolds number is large enough, typically <math>10^4</math>. Since the local Reynolds number for the plane jet continues to increase with increasing downstream distance, this state will eventually be reached. The higher the source Reynolds number, the closer to the exit plane the similarity of the moments will be realized. | ||
+ | |||
+ | == Other free shear flows == | ||
+ | |||
+ | The plane jet is but one of the flows which can be analyzed in the manner described above. A few of the possibilities which have already been analyzed are the axisymmetric jet, plane and axisymmetric wakes (George 1989, 1995 Hussein | ||
+ | et al. 1994, Moser et al. 1997). Other possibilities include free shear layers, thermal plumes and the self-propelled wake to mention but a few. All of these fall into the two categories described above: Flows which evolve at constant Reynolds number, and flows which do not. The axisymmetric jet and the plane wake are of the former type, and hence when properly scaled (using the techniques described above) will yield Reynolds number and source dependent solutions. These have been already discussed in detail in the cited papers and will not be discussed further here. The second type of flows (those for which the local Reynolds number varies with streamwise distances) also fall into two types: Those for which the local Reynolds number increases downstream (like the plane jet, plume or the shear layer), and those for which it decreases (like the axisymmetric wake). | ||
+ | |||
+ | When the Reynolds number is increasing with <math>x</math>, the flow will eventually reach | ||
+ | the state of full similarity where all of the mean and turbulence quantities collapse. This state will be characterized by the infinite Reynolds number dissipation relation <math>\epsilon \propto q^3 / L</math> which will be manifested in the approach of <math>d \delta /dx</math> to its asymptotic value. (This has been shown above to be constant for the plane jet, but will be different for wakes, for example.) Generally this approach will coincide with a turbulent Reynolds number of <math>q^4 / \epsilon \nu \sim 10^4</math> and the emergence in the spectra of the <math>{k}^{− 5 / 3}</math> range. Before this, the lack of collapse will be most evident in those quantities which depend directly on <math>d \delta /dx</math>, like <math> \left\langle -uv \right\rangle</math>, <math>\left\langle v^2 \right\rangle</math>, etc. Other quantities like the mean flow will collapse much earlier, and as noted above will collapse to profiles independent of source conditions. The latter will not be the case for the second moment quantities since a dependence in the asymptotic value of <math>d \delta /dx</math> will result in differences in the Reynolds stress equations themselves. | ||
+ | |||
+ | |||
+ | Perhaps the most troubling (and for that reason the most interesting) flows | ||
+ | are those where the local Reynolds number is decreasing — like the axisymmetric | ||
+ | wake. In these cases, the mean velocity profiles will collapse, assuming the | ||
+ | Reynolds number to be large enough that the viscous terms in the mean momentum | ||
+ | equation are negligible. The asymptotic growth rate (corresponding to <math>\epsilon \propto q^3 / \delta</math> or in the case of the axisymmetric wake, <math>\delta \propto x^{1/3}</math>) will only be achieved as long as the local Reynolds number is high enough (again <math>q^4 / \epsilon \nu \propto 10^4</math>). As soon as it drops below this value, the growth rate and scale parameters will begin to deviate (perhaps substantially) from the asymptotic power law forms. The turbulence quantities will begin to reflect this in the lack of collapse — again first noticable in quantities with <math>v^2</math>, etc. which have a direct dependence on <math>d \delta /dx</math>. | ||
+ | |||
+ | |||
+ | The mean velocity profile, however, will continue to collapse when scaled in local variables. The same will be true for flows in which the source Reynolds number is not high enough for the flow to ever achieve the requisite turbulence Reynolds number — the mean velocity will collapse even though the <math>x</math>-dependence will be all wrong — at least if asymptotic behavior is expected. |
Latest revision as of 13:42, 19 April 2016
Nature of turbulence |
Statistical analysis |
Reynolds averaged equation |
Turbulence kinetic energy |
Stationarity and homogeneity |
Homogeneous turbulence |
Free turbulent shear flows |
Wall bounded turbulent flows |
Study questions
... template not finished yet! |
Contents |
Introduction
Free shear flows are inhomogeneous flows with mean velocity gradients that develop in the absence of boundaries. Turbulent free shear flows are commonly found in natural and engineering environments. The jet of of air issuing from one's nostrils or mouth upon exhaling, the turbulent plume from a smoldering cigarette, and the buoyant jet issuing from an erupting volcano - all illustrate both the omnipresence of free turbulent shear flows and the range of scales of such flows in the natural environment. Examples of the multitude of engineering free shear flows are the wakes behind moving bodies and the exhausts from jet engines. Most combustion processes and many mixing processes involve turbulent free shear flows.
Free shear flows in the real world are most often turbulent. Even if generated as laminar flows, they tend to become turbulent much more rapidly than the wall-bounded flows which we will discuss later. This is because the three-dimensional vorticity necessary for the transition to turbulence can develop much more rapidly in the absence of walls that inhibit the qrowth velocity components normal to them.
The tendency of free shear flows to become and remain turbulent can be greatly modified by the presence of density gradients in the flow, especially if gravitational effects are also important. Why this is the case can easily be seen by examining the vorticity equation for such flows in the absence of viscosity,
| (1) |
The last term can act to either increase or decrease vorticity production but only in non-barotropic flows. (Recall that a barotropic flow is one in which the gradients of density and pressure are co-linear, because the density is a function of the pressure only). For example, in the vertically-oriented buoyant plume generated by exhausting a lighter fluid into heavier one, the principal density gradient is across the flow and thus perpendicular to the gravitational force which is the principal contributor to the pressure gradient. As a consequence the turbulent buoyant plume develops much more quickly than its uniform density counterpart, the jet. On the other hand, horizontal free shear flows in a stably stratified environment (fluid density decreases with height) can be quickly suppressed since the density and pressure gradients are in opposite directions.
Free turbulent shear flows are distinctly different from the homogeneous shear flows. In a free turbulent shear flow, the vortical fluid is patially confined and is separated from the surrounding fluid by an interface, the turbulent-nonturbulent interface (also known as the ”Corrsin superlayer” after itself discoverer). The turbulent/non-turbulent interface has a thickness which is characterized by the Kolmogorov microscale, thus its characterization as an interface is appropriate. The actual shape of the interface is random and it is severely distorted by the energetic turbulent processes which take place below it, with the result that at any given location the turbulence can be highly intermittent. This means that at a given location, it is sometimes turbulent, sometimes not.
It should not be inferred from the above that the non-turbulent fluid outside the superlayer is quiescent. Quite the opposite is true since the motion of the fluid at the interface produces motions in the surrounding stream just as would the motions of a solid wall. Alternately, the flow outside the interface can be viewed as being “induced” by the vortical motions beneath it. It is easy to show that these induced motions are irrotational. Thus since these random motions of the outer flow have no vorticity, they can not be considered turbulent.
Figure 7.1 shows records of the velocity versus time at a number of locations in the mixing layer of a round jet. When turbulent fluid passes the probes, the velocity signals are characterized by bursts of activity. The smooth undulations between the bursts are the irrotational fluctuations induced by the turbulent vorticity on the other side of the interface. Note that near the center of the mixing layer where the shear is a maximum, the flow is nearly always turbulent while it becomes increasingly intermittent as one proceeds away from the region of maximum production of turbulence energy. This increasing intermittency toward the outer edge is a characteristic of all free shear flows, and is an indication of the fact that the turbulent/non-turbulent interface is constantly changing its position.
One of the most important features of free shear flows is that the amount of fluid which is turbulent is continuously increased by a process known as entrainment. No matter how little fluid is in the flow initially, the turbulent part of the flow will continue to capture new fluid by entrainment as it evolves. The photograph of an air jet in Figure 1.2 illustrates this phenomenon dramatically. The mass flow of the jet increases at each cross-section due to entrainment. Entrainment is not unique to turbulent flows, but is also an important haracteristic of laminar flow, even though the actual mechanism of entrainment is quite different.
There are several consequences of entrainment. The first and most obvious is that free shear flows continue to spread throughout their lifetime. (That such is the case for the air jet of Figure 7.1 is obvious). A second consequence of entrainment is that the fluid in the flow is being continuously diluted by the addition of fluid from outside it. This is the basis of many mixing processes, and without such entrainment our lives would be quite different. A third consequence is that it will never be possible to neglect the turbulent transport terms in the dynamical equations, at least in the directions in which the flow is spreading. This is because the dilution process has ensured that the flow can never reach homogeneity since it will continue to entrain and spread through its lifetime (Recall that the transport terms were identically zero in homogeneous flows). Thus in dealing with free shear flows, all of the types of terms encountered in the turbulence kinetic energy equation of Chapter 4 must be dealt with — advection, dissipation, production, and turbulent transport.
Turbulent free shear flows have another distinctive feature in that they very often give rise to easily recognizable large scale structures or eddies. Figure 1.2 also illustrates this phenomenon, and coherent patterns of a scale equal to the lateral extent of the flow are clearly visible. These large eddies appear to control the shape of the turbulent/non-turbulent interface and play an important role in the entrainment process. They may also be important to the processes by which the turbulence gains and distributes energy from the mean flow.
A feature which free shear flows have in common with the homogeneous flows discussed in Chapter 6 is that their scales continue to grow as long as the flow remains turbulent. The dynamical equations and boundary conditions for many free shear flows can be shown to admit to similarity solutions in which the number of independent variables is reduced by one. According to the equilibrium similarity principle set forth in this chapter, such flows might be expected to asymptotically achieve such a state; and this, in fact, occurs. In the limit of infinite Reynolds number, some such flows can even be characterized by a single time and length scale, thus satisfying the conditions under which the simplest closure models might be expected to work. Care must be taken not to infer too much from the ability of a given closure model to predict such a flow, since any model which has the proper scale relations should work.
Finally there is the important question of whether free shear flows become asymptotically independent of their initial conditions (or source conditions). The conventional wisdom until very recently has been that they do. If correct, this means that there is nothing that could be done to alter the far downstream flow. There is recent theoretical and experimental evidence, however, that this traditional view may well be wrong. If so, this opens up previously un-imagined possibilities for flow control at the source.
In the remainder of this chapter, the averaged equations of motion will be simplified, and similarity solutions for several ideal shear flows will be derived and discussed in detail. The role of the large eddies will be discussed, and mechanisms for turbulent entrainment will be examined. The energy balance of several turbulent free shear flows will be studied in some detail. Finally, the effects of confinement and non-homogeneous external boundary conditions will be considered.
The averaged equations
The shear layer equations
One of the most important ideas in the history of fluid mechanics is that of the boundary layer approximation. These approximations to the Navier-Stokes equations were originally proposed by Prandtl in his famous theory for wall boundary layers. By introducing a different length scale for changes perpendicular to the wall than for changes in the flow direction, he was able to explain how viscous stresses could survive near the wall at high Reynolds number. These allowed the no-slip condition at a surface to be satisfied, and resolved the paradox of how there could be drag in the limit of zero viscosity.
It may seem strange to be talking about Prandtl’s boundary layer idea in a section about free shear flows, but as we shall see below, the basic approximations can be applied to all “thin” (or slowly growing) shear flows with or without a surface. In this section, we shall show that free shear flows, for the most part, satisfy the conditions for these “boundary layer approximations”. Hence they belong to the general class of flows referred as “boundary layer flows”.
One important difference will lie in whether momentum is being added to the flow at the source (as in jets) or taken out (by drag, as in wakes). A related influence is the presence (or absence) of a free stream in which our free shear flow is imbedded. We shall see that stationary free shear flows fall into two general classes, those with external flow and those without. One easy way to see why this makes a difference is to remember that these flows all spread by entraining mass from the surrounding fluid. You don’t have to think very hard to see that the entrained mass is carrying its own momentum into the shear flow. You should expect (and find) that even a small free stream velocity can make a significant difference, since the momentum carrried in is mixed in with that of the fluid particles which are already part of the turbulence. The longer the flow develops (or the farther downstream one looks), the more these simple differences can make a difference in how the flow spreads. In view of this, it should be no surprise that the presence or absence of an external stream plays a major role in petermining which mean convection terms which must be retained in the governing equations.
We will consider only flows which are plane (or two-dimensional) in the mean (although similar considerations can be applied to flows that are axisymmetric in the mean). In effect, this is exactly the same as assuming the flow is homogeneous in the third direction. Also we shall restrict our attention to flows which are statistically stationary, so that time derivatives of averaged quantities can be neglected. And, of course, we have already agreed to confine our attention to Newtonian flows at constant density.
It will be easier to abandon tensor notation for the moment, and use the symbols for the streamwise direction, mean and fluctuating velicities respectively, and for cross-stream. Given all this, the mean momentum equations reduce to:
x-component:
| (2) |
y-component:
| (3) |
In addition, we have the two-dimensional mean continuity equation which reduces to:
| (4) |
Order of magnitude estimates
Now let's make an order of magnitude estimate for each of the terms. This procedure may seem trivial to some, or hopeless hand-waving to others. The fact is that if you fall either of these groups, you have missed something important. Learning to make good order-of-magnitude arguments and knowing when to use them (and when not to use them) are two of the most important skills in fluid mechanics, and especially turbulence. To do this right we will have to be very careful to make sure our estimates accurately characterize the terms we are making them for.
Naturally we should not expect changes of anything in the - direction to scale the same as changes in the -direction, especially in view of the above. (This, after all, is the whole idea of "thin" shear flow.) So let's agree that we will pick a length scale, say , characteristic of changes or mean quantities in the -direction; i.e.
| (5) |
where for now '' means "of order of magnitude of". And we can do the same thing for changes of mean quantities in the -direction by defining a second length scale, say , to mean:
| (6) |
Note that both these scales will vary with the streamwise position where we evaluate them. A good choice for might be proportional to the local lateral extent of the flow (or its "width"), while is related to the distance from the source.
Consider the mean velocity in equation 10 in section "Wall-bounded turbulent flows". It occurs in five different terms: alone; twice with -derivatives, and ; and twice with -derivatives, and . Now it would be tempting to simply pick a scale velocity for say , and use it to estimate all five terms, say as: , , , , and . But this is much too naive, and fails to appreciate the true role of the terms we are evaluating.
Look at the examples shown in Figure ??. Our simple approach would provide an appropriate estimate for a jet if we took our velocity scale equal to the mean centerline velocity at a given streamwise location; i.e., . This is because both the mean velocity and the changes in the mean velocity across the flow are characterized by the centerline velocity. But by contrast, look at the wake. Even relatively close to the wake generator, the wake deficit, is small compared to the free stream velocity (i.e. ). So the obvious choice to scale is On the other hand, an estimate for the velocity gradient across the flow of would be much too big, again because the deficit is so small. Obviously a better choice would be to use the centerline mean velocity deficit ; i.e.,
| (7) |
In the order of magnitude analysis below, we shall try to keep the discussion as general as possible by using to characterize the mean velocity when it appears by itself, and to represent changes in the mean velocity. For the jet example of the preceding paragraph, both and are the same; i.e. and . But for the wake they are different because of the external stream; i.e., and . If you can keep in mind why these differences exist among the various flows, it will be a lot easier to both understand the results and not confuse them.
Now we could distinguish changes of velocity in the -direction from those in the - direction. But this level of complexity is not necessary (at least for the examples considered here), especially since we have left the precise defenition of rather nebulous. What we can do is to use the same estimate for changes in the velocity as for the - direction, and define our length scale to make the estimate based on both correct; i.e. . To see why this make sence physically and can be reasoned (as opposed to guessed), let's look at the wake. Pick a spot outside, near the edge of the wake fairly close the generator (say point A). Now proceed at constant far downstream in . Eventually the wake will have spread past you and you will be close enough to the centerline so the local mean velocity will be closer to than . Obviously we have simply traveled far enough at constant to ensure that is the proper scale for the changes in velocity in the - direction. If the distance downstream over which this change occured is taken as , then proper estimate is easily seen to be:
| (8) |
But this is exactly what we would have gotten by taking as we agreed above. We simply have absorbed any differences into our choice of . When considering specific problems and applying similarity techniques, the seemingly arbitrary choices here become quite precise constraints (as we shall see).
We still haven't talked about how to estimate the velocity scale for , the cross-stream mean velocity component. From the continuity equation we know that:
| (9) |
From our considerations above, we know that:
| (10) |
If there is no mean cross flow in the external stream, then the scale for is the same as the scale for changes in . Therefore,
| (11) |
It follows immediately that the order of magnitude of the cross-stream velocity is:
| (12) |
We might have expected something like this if we had thought about it. If the -velocity were of the same order as the - velocity, how could the flow in any sense be a "thin shear flow". On the other hand, it also makes sense that since both are some measure of how the flow spreads/ Note that equation 12 would not be correct estimate for if there were an imposed cross-flow, since then we would have to consider and changes in separately (exactly as for ).
The mean pressure gradient term is always a problem to estimate at the outset. Therefore it is better to simply leave this problem alone, and see what is left at the end. In the estimates below you will see a question mark, which simply means we are postponing judgement until we have more information. Sometimes we will have to keep the term simply because we don't know enough to throw it away. Other times it will be obvious that it must remain because there is only one term left that must be balanced by something.
Now we have figured out how to estimate the order of magnitude of all the terms except the turbulence terms. For most problems this turns out to be pretty straight-forward if you remember our discussion of the pressure strain-rate terms. They so effectively distribute the energy that even that the three components of velocity are usually about the same order of magnitude. So if we choose a turbulence scale as simply , then , , . But what about the Reynolds shear stress components like ? When acting against the mean shear to produce energy, they tend to be well-correlated and the maximum value of . Obviously the right choice though since some cross-stress like are identically zero in plane flows because of homogeniety. Allowing for the possibility of a Reynolds number dependent correlation coefficient will be sent to be especially important when the terms in the turbulence kinetic energy and dissipation equations are considered.
Streamwise momentum equation
Let's look now at the -component of the mean momentum equation and write below each term its order of the magnitude.
|
|
Now the survival of at least one of the terms on the left-hand side of the equation is the essence of free shear flow. Since the largest of these is the first, let's re-scale the others by deviding all the estimates by . Now we have:
|
|
Just in case you need to be reminded, this is a turbulence course. Therefore almost all the interesting flows are at high Reynolds number. But what does this mean? Now we can say. It could mean the viscous terms in the mean momentum equation are negligible, which in tuurn means that:
|
and
|
Obviously the second criterion is much more stringent, since is typically of order 10. When , the contributions of the viscous stresses are certainly negligible, at least as far as the -component of the mean mometum equation is concerned. But if only, this is a pretty marginal assumption, and you might want to retain the last term in your analysis. Such is unfortunately the case in many experiments that are often at quite low Reynolds numbers.
So, what then, tou ask, do we do about the turbulence terms? To repeat this is a TURBULENCE course! Which means there is no way we are going to throw away all the turbulence terms on the right-hand side of the equation. There is no magic or hocus-pocus about this, you simply can't have turbulence without turbulence terms. And if there is no turbulence, we are really aren't too interested. So, what does this tell us? It tells us about ! Surprised? I bet you thought it would tell us about , right? Nope. Look at the right-hand side of the second equation of this subsection and the orders of magnitude below it. Almost always ; sometimes much less and certainly no bigger. Then the biggest turbulence term is one involving the Reynolds shear stress which we have estimated as . Hence there can be no turbulence terms at all unless:
|
Wow! Look what we have learned about turbulence without ever solving a single equation. We know how the growth of our free shear flow relates to the turbulence intensity - at least in terms of order of magnitude. The simlarity theories described below will be able to actually deduce the - dependence of , again without solving the equations.
Does the argument above mean the other turbulence term is negligible since it is only of order ? Well, no question it is usually small, typically less than 5%. But unlike the viscous terms, this term often does not get smaller the father we go in the streamwise direction. So the answer depends on the question we are asking and how accurately we want our answer to be. If we are asking only first order questions and are willing to accept errors in our solutions (or experiments) to within a 5-10 percent, then this normal stress is certainly negligible. Before you are ready to forget about this, consider the following: since typically, to first order free shear flows do not grow at all! So if you want to talk about whether for a particular shear flow is 0.095 instead of 0.09 (turbulence wars have been waged over less), then you better be dealing with a full desk - meaning these second order terms had better be retained in any calculation or experimental verification (like momentum or energy balances).
Much of the confusion in the literature about free shear flows comes from the failure to use second order accuracy to make second order statements. And even more unbelievably, most attempts to do so-called code validation of second-order turbulence models are based on mesurements which are only first order accurate. Now I know you think I am kidding - but try to do a simple momentum balance on the data to within order . Did I hear you say it was ridiculous to use second-order equations (like the Reynolds stress equations or even the kinetic energyequation) and determine constants to three decimal places using first order accurate data? Shhhh!!!! Your collegues might be annoyed unless you keeep this a secret.
Hopefully you have followed all this. If not, go forward and "faith" will come to you. I'll summarize how we do this:
- You have to look very carefully at the particular flow you wish to analyze.
- Then, based on physical reasoning and data (if you have it), you make estimates you think are appropriate.
- Then you use these estimates to decide which terms you have to keep, and which you can ignore.
- And when you are completely finished, you carefully look back to make sure that all the terms you kept really are important and remain important as the flow develops.
Finally, if you are really good at this, you look carefully to make sure that all the terms you kept do not all vanish at the same place (by going from positive to negative, for example). If this happens, then the terms you neglected may in fact be the only terms left in the equation thet are not exactly zero - oops! This can make fo some very interesting situations/ (The famous "critical layer" of boundary layer stability theory is an example of this.)
The transverse momentum equation
Now we must also cosider the transverse or -momentum equation. The most important thing we have to remember is that the - equation can NOT be considered indepedently from the -momentum equation. We would never consider trying to simplify a vector like by dividing only one component by say to produce , since this would destroy the whole idea of a vector as having a direction. Similarly we must not do this for a vector equation. We have already decided thet the first term on the left-hand side of the -momentum equation was the term we had to keep, and we divided by its order of magnitude, , to make sure it was of order one. Thus we have already decided how we are going to scale all the components of the vector equation, and so we must do exactly the same thing here. but first we must estimate the order of magnitude of each term, exactly as before. Using our previous results, here's what we get:
|
|
Now dividing each term by , exactly as before, we obtain:
|
|
Unless you have seen this all before, the left-hand side is probably a surprise: none of the mean convection terms are of order one! In fact the only estimated term that is of order one in the whole equation is and only because we have already agreed that had to be of order one to keep a turbulence term in the -momentum equation. Of course, there cannot be an equation with only a single term equal to zero, unless we have badly over-estimated its order of magnitude. Fortunately there is the pressure term left to balance it, so to first order in , the -momentum equation reduces to simply:
| (13) |
This simple equation has an equally simple interpretation. It says thet the change in the mean pressure is only due to the radisl gradient of the transverse component of the Reynolds normal stress, .
We can integrate equation 13 across the shear layer from a given value of to infinity (or anywhere else for thet matter) to obtain:
| (14) |
If the free stream is at constant (or zero) mean velocity, then , so we can write simply:
| (15) |
Either of these can be substituted into -momentum equation to eliminate the pressure entirely, as we shall show below.
Now since typically , at least to first order in . Therefore one might be tempted to conclude that the small pressure changes across the flow are not important. And they are not, of course, if only the streamwise momentun is considered. But without this small mean pressure gradient across the flow, there would be no entrainment and no growth of the shear layer. In other words, the flow would remain parallel. clearly whether an effect is negligible or not depends on which question is being asked.
The free shear layer equations
If we assume the Reynolds number is always large enough that the viscous terms can be neglected, then as noted above, the pressure term in the terms of and . Thus, to second-order in or , the momentum equations for a free shear flow reduce to a single equation:
| (14) |
As we have seen above, the second term on the left-hand side may or may not be important, depending on whether or . Clearly this depends on the velocity deficit (or excess) relative to the free stream. The term in brackets on the right-hand side is second-order (in ) compared to the others, and so could rigthgly be neglected. It has been retained for now since in some flows it does not vanish with increasing distance from the source. Moreover, it can be quite important when considering integrals of the momentum equation, since the profiles of and do not vanish as rapidly with increasing as do the other terms.
Two-dimensional Turbulent Jets
Turbulent jets are generated by a concentrated source of momentum issuing into an ambient environment. The state of the environment and the external boundary conditions are very important in determining how the turbulent flow evolves. The simplest case to consider is that in which the environment is at rest and unbounded, so all the boundary conditions are homogeneous. The jet itself can be assumed to isssue from a line source or a slot as shown in figure ??. The jet then evolves and spreads downstream by entraining mass from the surroundings which are at most in irrotational motion induced by the vortical fluid in the jet. But no new momentum is added to the flow, and it is this fact that distinguishes the jet from all other flows. As we shall see below the rate at which momentum crosses any - plane is not quite constant at the source value, however, due to the small streamwise pressure gradient arising from the turbulence normal stresses.
The avereged continuity equation can be integrated from the centerline (where it is zero by symmetry) to obtain , i.e.,
| (17) |
It immediately follows that the velocity at is given by:
| (18) |
Twice the integral represents the total volume flow rate crossing a given - plane. Therefore it makes sense that the rate of increase of this integral is the enrainment velocity, . Obviously, cannot be zero if the jet is to spread, at least in two-dimensional flow. It is easy to see that the entrainment velocity calculated from the integral is consistent with our order of magnitude estimate above (i.e., ). Note that the integral makes no sense if the mean velocity does not go to zero with increasing .
If the free stream is assumed to have zero streamwise velocity, then the scales for the velocity and gradienits of it are the same, so the equations reduce to simply equation 16 in the next section(need to edit) and the mean continuity equation. The latter can be multiplied by to yield:
| (19) |
When this is added to equation 16 from the next section(need to edit) the terms can be combined to obtain:
| (20) |
This can be integrated across the flow for any given value of to obtain:
| (21) |
where we have assumed that , , and vanish as . (Remember these assumptions the next time someone tries to tell you that a small co-flowing stream and modest background turbulence level in the external flow are not important.)
Equation (21) can in turn be integrated from the source, say , to obtain:
| (22) |
where is the rate at which momentum per unit mass per unit length is added at the source. The first two terms of the itegral are the flux of streamwise momentum due to the mean flow and the turbulence. The last term is due to the streamwise pressure gradient obtained by integrating the -momentum equation.
This is a very powerful constraint, and together with the reduced governing equations allows determination of a similarity solution for this problem. Equation 22 also is of great value to experimentalists, since it can be used to confirm whether their measurements are valid and whether the flow indeed approximates a plane jet in an infinite environment at rest. Be fore-warned, most do not. And the explanations as to why range from hilarious to pathetic. Consider this, if you measure the velocity moments to within say 1%, then shouldn’t you be able to estimate the integral even more accurately, since the random errors from the integration should be reduced by one over the square root of the number of points you used in the integration? Be very wary of measurements which have not been tested against the momentum integral — for this and all flows. There is usually a reason the authors fail to mention it.
Similarity analysis of the plane jet
Foreword The results of this section are VERY MUCH at odds with those in all texts, including Pope 2000. There is an increasing body of literature to support the position taken, however, especially the multiplicity of solutions and the role of upstream or initial conditions. Since these ideas have a profound effect on how we think about turbulence, it is important that you consider both the old and the new, and be aware of why the conclusions are different.
The message you should take from the previous section is that regardless of how the plane jet begins, the momentum integral should be constant at . Note that Schneider (1985) has shown that the nature of the entrained flow in the neighborhood of the source can modify this value, especially for plane jets. Also the term , while important for accounting for the momentum balance in the analysis below with no loss in generality.
The search for similarity solutions to the averaged equations of motion is exactly like the approach utilized in laminar flows (eg. Batchelor 1967) except that here there are more unknows than equations, i.e. the averaged equations are not closed. Thus solutions are sought which are of the form
| (23) |
| (24) |
| (25) |
where
| (26) |
The scale functions are functions of only, i.e. and only. Note that the point of departure from the traditional turbulence analyses was passed when it was not arbitrarily decided that , etc. (c.f. Tennekes and Lumley 1972, Townsend 1976, Pope 2000).
Now let you think there is something weird about the way we have written the form of the solutions we say we are looking for, consider this. Suppose you were looking for some - dependent quantities you could normalize the profiles by to make them all collepse together - like any experimentalist almost always does. Such a search is called the search for "scales" , meaning "scales" which collapse the data. Now precisely what does "collapse the data" really mean? Exactly this: "collapse the data" means that the form of the solution must be an -independent solution to the governing equations. But this is exactly what we are looking for in equations 23 and 24. Obviously we need to plug them into the averaged equations and see if such solution is possible. We already know which terms in these equations are important. What we want to do now is to see if these govering equations admit to solutions for which these terms remain important throughout the jet development. These are called equilibrium similarity solutions. Note that for any other form of solution, one or more terms could die off, so that others dominate.
Differentiating equations 23 and 24 using the chain-rule, substituting into the averaged momentum equation, and clearing the terms yields:
| (27) |
Now, all of the -dependence is in the square bracketed terms, so equilibrium similarity solutions are possible only if all the bracketed terms have the same -dependence, i.e.,
| (28) |
The coefficients of proportionality, say and respectively, can at most depend on the details of how the flow began; i.e., the mistery argument "*".
Substitution into the momentum integral yields to first order,
| (29) |
where
| (30) |
and are parameters which also can at most depend on how the flow began (since every other dependence has been eliminated). Moreover, since is independent of , then and it follows that:
| (31) |
It's easy to see that all of the relations involving are proportional to , leaving only
| (32) |
Thus the only remaining necessary condition for similarity is
| (33) |
Note that, unlike the similarity solutions encountered in laminar flows, it's possible to have a jet which is similar without having some form of power law behavior. In fact, the -dependence of the flow may not be known at all because of the closure problem. Nonetheless, the profiles will collapse with the local length scale. In fact, it is easy to show that the profiles for all plane jets will be alike, if normalized properly, even if the growth rates are quite different. The scale velocity can be defined to be the centerline velocity by absorbing an appropriate factor into the profile function . Also, the entire factor in brackets on the right-hand side of equation 32 can be absorbed into the Reynolds stress profile, . Finally, if the length scale is chosen the same for all flows under consideration (e.g., the half-width, say , defined as the distance between the center and where ), then the similarity equation govering all jets reduces to,
| (34) |
where both has been modified so that and is just equal to with the factor absorbed into it. Thus all jets, regardless of initial conditions, will produce the same similarity velocity profile, when scaled with and . This will be not true for the Reynolds stress, however, unless the coefficient of proportionality between and is the same for all plane jets. In general, it is not; so here is where the dependence on source conditions shows up.
It is immediately obvious from equation 32 that usual arbitrary choice of considerably restricts the possible solutions to those for which , or plane jets which grow linearly. However, eben if the growth were linear under some conditions, there is nothing in the theory to this point which demands that the constant be universal and independent of how the jet began.
The idea that jets might all be described be the same asymptotic growth rate stems from the idea of a jet formed by a point source of momentum only, say . Such a source must be infinitesimal in size, since any finite size will also provide a source of mass flow. In the absence of scaling parameters other than distance from the source, , the only possibilities are , , . Obviously , and the constants of proportionality can be assumed "universal" since there is only one way to create a point source jet, and this is also contarary to the usual assumption.
The problem with the point source idea is not in the idea itself, but that it has been accepted for so long by so many without justification as the asymptotic state of finite source jets. The usual line of argument (eg. Monin and Yaglom 1972) if one is presented at all is as follows. The jet entrains mass and the rate at which this occurs increases with distance downstream, so that the rate of mass entrainment ultimately overwhelms that added at the source. It is then inferred that the latter can be neglected. And what is the experimental proof for all this? Precisely that the profiles of collapse to a single profile for all types of initial conditions. But as was demonstrated above, the mean momentum equation tells us that the mean velocity profiles will collapse to a single profile no matter what! If there is a dependence on the initial conditions it can only show up in the other turbulence moments, in profiles of the Reynolds shear stress normalized by (instead of ), and, of course, in itself. And guess what quantities no two experiments ever seem to agree on? — All of the above! Could all those experiments plus this nice theory with virtually no assumptions be wrong? Many think so. It is truly amazing how much good science people will throw away in order to hang-on to a couple of ideas that were never derived in the first place but assumed; namely, and the asymptotic independence of initial conditions. But then there are still many people who believe the world is flat.
George (1989) demonstrated for an axisymmetric jet that even simple dimensional analysis suggests that the classical theory is wrong, and the same arguments apply here. Suppose that in addition to momentum, mass is added at the source at a rate of . Now there is an additional parameter to be considered, and as a consequence, an additional length scale given by . Thus the most that can be inferred from dimensional analysis is that , and are functions of , with no insight offered into what the functions might be.
Implications of the Reynolds stress equations
George 1989 further argued that some insight into the growth rate of the jet can be obtained by considering the conditions for similarity solutions of the higher moment equations, in particular the kinetic energy equation. The choice of the kinetic energy equation for further analysis was unfortunate since it implicitly assumed that the three components of kinetic energy all scaled in the same manner. This is, in fact, true only if , which is certainly not true a priory for turbulent shear flows. Therefore, here the individual component equations of the Reynolds stress will be considered (as should have been done in George 1989).
For the plane jet the equation for can be written to first oder (Tennekes and Lumley 1972) as
| (35) |
| (35) |
where is the energy dissipation rate for .
By considering similarity forms for the new moments like
| (36) |
| (37) |
| (38) |
| (39) |
and using , it's easy to show that similarity of the -equation is possible only if
| (40) |
| (41) |
| (42) |
| (43) |
All of these are somewhat surprising: The first (even though a second moment like the Reynolds stress) because the factor of is absent; the second, third and forth because it is present.
Similar equations can be written for the ,, and -equations; i.e.
| (44) |
| (44) |
| (45) |
| (45) |
| (46) |
| (46) |
When each of the terms in these equations is expressed in similarity variables, the resulting similarity conditions are:
| (47) |
| (48) |
| (49) |
| (50) |
and the real surprise,
| (51) |
There is an additional equation which must be accounted for: namely that the sum of the pressure strain-rate terms in the component energy equations be zero (from ontinuity). Thus,
| (52) |
or in similarity variables,
| (53) |
This can be true for all only if
| (54) |
An immediate consequence is that
| (55) |
From equations 41,47 and 48 it also follows that the constraint imposed by 54 can be satisfied only if
| (51) |
But from equation 51, this can be true only if
| (57) |
The relations given by equations 56 and 55 were assumed without proof in the George 1989 analysis. The additional constraint imposed by equation 57 was not derived, however, and arises from the additional information provided by the pressure strain-rate terms.
Hence, similarity solutions of the Reynolds stress equations are possible only if
| (58) |
It is an immediate consequence of the earlier discussion on the nature of the dissipation that there are only two possibilities for this to occur:
- Either the local Reynolds number of the flow is constant so that the effect of
the dissipation on the energy containing eddies (and those producing the Reynolds stress as well) does not vary with downstream distance; or
- The local turbulence Reynolds number is high enough so that the relation is approximately valid (for a physical length ).
Unlike some flows (like the axisymmetric jet or plane wake) where the local Reynolds number is constant, for the plane jet it varies with downstream distance. Therefore the only possibility for similarity at the level of the Reynolds stresses is (ii). This can occur only when the turbulence Reynolds number is large enough, typically . Since the local Reynolds number for the plane jet continues to increase with increasing downstream distance, this state will eventually be reached. The higher the source Reynolds number, the closer to the exit plane the similarity of the moments will be realized.
Other free shear flows
The plane jet is but one of the flows which can be analyzed in the manner described above. A few of the possibilities which have already been analyzed are the axisymmetric jet, plane and axisymmetric wakes (George 1989, 1995 Hussein et al. 1994, Moser et al. 1997). Other possibilities include free shear layers, thermal plumes and the self-propelled wake to mention but a few. All of these fall into the two categories described above: Flows which evolve at constant Reynolds number, and flows which do not. The axisymmetric jet and the plane wake are of the former type, and hence when properly scaled (using the techniques described above) will yield Reynolds number and source dependent solutions. These have been already discussed in detail in the cited papers and will not be discussed further here. The second type of flows (those for which the local Reynolds number varies with streamwise distances) also fall into two types: Those for which the local Reynolds number increases downstream (like the plane jet, plume or the shear layer), and those for which it decreases (like the axisymmetric wake).
When the Reynolds number is increasing with , the flow will eventually reach the state of full similarity where all of the mean and turbulence quantities collapse. This state will be characterized by the infinite Reynolds number dissipation relation which will be manifested in the approach of to its asymptotic value. (This has been shown above to be constant for the plane jet, but will be different for wakes, for example.) Generally this approach will coincide with a turbulent Reynolds number of and the emergence in the spectra of the range. Before this, the lack of collapse will be most evident in those quantities which depend directly on , like , , etc. Other quantities like the mean flow will collapse much earlier, and as noted above will collapse to profiles independent of source conditions. The latter will not be the case for the second moment quantities since a dependence in the asymptotic value of will result in differences in the Reynolds stress equations themselves.
Perhaps the most troubling (and for that reason the most interesting) flows
are those where the local Reynolds number is decreasing — like the axisymmetric
wake. In these cases, the mean velocity profiles will collapse, assuming the
Reynolds number to be large enough that the viscous terms in the mean momentum
equation are negligible. The asymptotic growth rate (corresponding to or in the case of the axisymmetric wake, ) will only be achieved as long as the local Reynolds number is high enough (again ). As soon as it drops below this value, the growth rate and scale parameters will begin to deviate (perhaps substantially) from the asymptotic power law forms. The turbulence quantities will begin to reflect this in the lack of collapse — again first noticable in quantities with , etc. which have a direct dependence on .
The mean velocity profile, however, will continue to collapse when scaled in local variables. The same will be true for flows in which the source Reynolds number is not high enough for the flow to ever achieve the requisite turbulence Reynolds number — the mean velocity will collapse even though the -dependence will be all wrong — at least if asymptotic behavior is expected.