Introduction to turbulence/Free turbulent shear flows

From CFD-Wiki

Revision as of 21:45, 14 November 2010

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,

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:

y-component:

In addition, we have the two-dimensional mean continuity equation which reduces to:

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.

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 lenght scale, say , to mean:

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.,

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:

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:

From our considerations above, we know that:

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,

It follows immediately that the order of magnitude of the cross-stream velocity is:

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