https://www.cfd-online.com/W/index.php?title=Special:Contributions/Pavitran&feed=atom&limit=50&target=Pavitran&year=&month=CFD-Wiki - User contributions [en]2016-12-10T04:46:49ZFrom CFD-WikiMediaWiki 1.16.5https://www.cfd-online.com/Wiki/User:PavitranUser:Pavitran2010-03-24T16:21:05Z<p>Pavitran: </p>
<hr />
<div>I am Dynampally Pavitran, I completed my MS in Computational and Experimental Turbulence from Chalmers University of Technology,Sweden. I have about 4 years of experience of using CFD in various domains like Aerospace, Wind Turbine and Hydroturbines. Presently I am working for a Hydroturbine company in India.</div>Pavitranhttps://www.cfd-online.com/Wiki/User:PavitranUser:Pavitran2007-01-04T08:52:18Z<p>Pavitran: </p>
<hr />
<div>I am Dynampally Pavitran, I completed my MS in Computational and Experimental Turbulence from Chalmers University of Technology,Sweden.</div>Pavitranhttps://www.cfd-online.com/Wiki/User:PavitranUser:Pavitran2006-07-24T04:09:43Z<p>Pavitran: </p>
<hr />
<div>I am Dynampally Pavitran, I completed my MS in Computational and Experimental Turbulence from Chalmers University of Technology,Sweden. Presently I am working as Design Engineer in Hindustan Aeronautics Limited(HAL).</div>Pavitranhttps://www.cfd-online.com/Wiki/User:PavitranUser:Pavitran2006-03-22T15:24:07Z<p>Pavitran: </p>
<hr />
<div>I am Dynampally Pavitran, I completed my Masters in Computational and Experimental Turbulence at Chalmers. Presently I am doing my PhD in Spain at Universidad politecnica de valencia, my PhD research is about Combustion in IC engines.</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-10-20T06:39:25Z<p>Pavitran: /* 2.4 The Origins of Turbulence */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a [[Direct numerical simulation (DNS)| direct numerical simulation]] in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this won't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems.<br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
=== The Equations Governing the Instantaneous Fluid Motions ===<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u}_i}{\partial t}+\tilde{u}_j\frac{\partial \tilde{u}_i}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T}_{ij}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u}_j\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u}_j}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T}_{ij}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u}_j\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u}_i}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T}_{ij}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T}_{kk}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T}_{ij}^{(v)}= 2\mu[\tilde{s}_{ij}-\frac{1}{3}\tilde{s}_{kk}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defined by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s}_{ij}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T}_{ij}^{(v)}= 2\mu\tilde{s}_{ij}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u}_i}{\partial t}+\tilde{u}_j\frac{\partial \tilde{u}_i}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u}_i}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u}_i}{\partial t}+\tilde{u}_j\frac{\partial \tilde{u}_i}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T}_{ij}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
=== Equations for the Average Velocity ===<br />
<br />
<br />
Although laminar solutions to the equations often exist that are consistent with the boundary conditions, perturbations to these solutions(sometimes even infinitesimal) can cause them to become turbulent. To see how this can happen, it is convenient to analyze the flow in two parts, a mean(or average) component and a fluctuating component. Thus the instantaneous velocity and stresses can be written as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\tilde {u}_i=U_i+u_i<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde p=P+p<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde T_{ij}^{(v)}=T_{ij}^{(v)}+\tau_{ij}^{(v)}<br />
&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.11)&lt;/td&gt;&lt;/tr&gt;<br />
&lt;/table&gt;<br />
<br />
<br />
Where &lt;math&gt;U_i&lt;/math&gt;, &lt;math&gt;P&lt;/math&gt; and &lt;math&gt;T_{ij}^{(v)}&lt;/math&gt; represent the mean motion, and &lt;math&gt;u_i&lt;/math&gt;, &lt;math&gt;p&lt;/math&gt; and &lt;math&gt;\tau_{ij}^{(v)}&lt;/math&gt; the fluctuating motions. This technique for decomposing the instantaneous motion is referred to as the '''''Reynolds decomposition.''''' Note that if the averages are defined as ensemble means, they are, in general, time-dependent. For the remainder of this material unless other wise stated, the density will be assumed constant so&lt;math&gt;\tilde{\rho}\equiv\rho&lt;/math&gt;,and its fluctuation is zero.<br />
<br />
Substitution of equations 2.11 into equations 2.10 yields<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial (U_i+u_i)}{\partial t}+(U_j+u_j)\frac{\partial (U_i+u_i)}{\partial x_j}\right] = -\frac{\partial (P+p)}{\partial x_i}+\frac{\partial (T_{ij}^{(v)}+\tau_{ij}^{(v)})}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.12)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This equation can now be averaged to yield an equation expressing momentum conservation for the averaged motion. Note that the operations of averaging and differentiation commute; i.e., the average of a derivative is the same as the derivative of the average. Also the average of a fluctuating quantity is zero. Thus the equation for the averaged motion reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\rho\left \langle u_j\frac{\partial u_i }{\partial x_j} \right \rangle&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.13)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the remaining fluctuating product term has been moved to the right hand side of the equation. Whether or not the last term is zero like the other fluctuating term depends on the correlation of the terms in the product. In general, these correlations are not zero.<br />
<br />
The mass conservation equation can be similarly decomposed. In incompressible form, substitution of equations 2.11 into equation 2.4 yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial (U_j+u_j)}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.14)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
of which average is <br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial U_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.15)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
It is clear from the equation 2.15 that the averaged motion satisfies the same form of the mass conservation equation as does the instantaneous motion at least for incompressible flows. How much simpler the turbulence problem would be if the same were true for the momentum! Unfortunately, as is easily seen from the equation 2.13, such is not the case.<br />
<br />
Equation 2.15 can be subtracted from equation 2.14 to yield an equation for instantaneous motion alone; i.e,<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial u_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.16)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Again, like the mean, the form of the original instantaneous equation is seen to be preserved. The reason, of course, is obvious: the continuity equation is linear. The momentum equation , on the other hand, is not; hence the difference.<br />
<br />
Equation 2.16 can be used to rewrite the last term in equation 2.13 for the mean momentum. Multiplying equation 2.16 by &lt;math&gt;u_i&lt;/math&gt; and averaging yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.17)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This can be added to:&lt;math&gt;\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle&lt;/math&gt; to obtain:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle +0=\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle+ \left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle =\frac{ \partial}{\partial x_j}{\left \langle u_iu_j\right \rangle} &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.18)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Where again the fact that arithmetic and averaging operations commute has been used.<br />
<br />
The equation for the averaged momentum, equation 2.13 can now be rewritten as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\frac{ \partial}{\partial x_j}{(\rho\left \langle u_iu_j\right \rangle)}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.19)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The last two terms on the right hand side are both divergence terms and can be combined; the result is:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial }{\partial x_j}[T_{ij}^{(v)}-{\rho\left \langle u_iu_j\right \rangle}]&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.20)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Now the terms in square brackets on the right have the dimensions of stress. The first term is, in fact , the viscous stress. The second term, on the other hand, is not a stress at all, but simply a re-worked version of the fluctuating contribution to the non-linear acceleration terms. The fact that it can be written this way, however, indicates that at least as far as the motion is concerned, it acts as though it were a stress- hence its name, the '''Reynolds stress'''. In the succeeding sections the consequences of this difference will be examined.<br />
<br />
=== The Turbulence Problem ===<br />
<br />
It is the appearance of the Reynolds stress which makes the turbulence problem so difficult - at least from the engineers perspective. Even though we can pretend it is a stress, the physics which give rise to it are very different from the viscous stress. The viscous stress can be related directly to the other flow properties by constitutive equations, which in turn depend only on the properties of the fluid( as in equation 2.5 for a Newtonian fluid). The reason this works is that when we make such closure approximations for a fluid, we are averaging over characteristic length and time scales much smaller than those of the flows we are interested in. Yet at the same time, these scales are much larger than the molecular length and time scales which characterize the molecular interactions that are actually causing the momentum transfer. (This is what the continuum approximation is all about).<br />
<br />
The '''''Reynolds stress''''', on the other hand, arises from the flow itself! Worse, the scales of the fluctuating motion which give rise to it are the scales we are interested in. This means that the closure ideas which worked so well for the viscous stress, should not be expected to work too well for the Reynolds stress. And as we shall see, they do not.<br />
<br />
This leaves us in a terrible position. Physics and engineering are all about writing equations(and boundary conditions) so we can solve them to make predictions. We don't want to have a build prototype airplanes first to see if they will they fall out of the sky. Instead we want to be able to analyze our designs before building, to save the cost in money and lives if our ideas are wrong. The same is true for dams and bridges and tunnels and automobiles. If we had confidence in our turbulence models, we could even build huge one-offs and expect them to work the first time. Unfortunately, even though turbulence models have improved to the point where we can use them in design, we still cannot trust them enough to eliminate expensive wind tunnel and model studies. And recent history is full of examples to prove this.<br />
<br />
The turbulence problem (from the engineers perspective) is then three-fold:<br />
<br />
* '''The averaged equations are not closed.''' Count the number of unknowns in equation 2.20 above. Then count the number of equations. Even with the continuity equation we have atleast six equations too few.<br />
<br />
* '''The simple ideas to provide the extra equations usually do not work.''' And even when we can fix them up for a particular class of flows (like the flow in a pipe, for example), they will most likely not be able to predict what happens in even a slightly different environment(like a bend).<br />
<br />
*'''Even the last resort of compiling engineering tables for design handbooks carries substantial risk.''' This is the last resort for the engineer who lacks equations or cannot trust them. Even when based on a wealth of experience, they require expensive model testing to see if they can be extrapolated to a particular situation. Often they cannot, so infinitely clever is Mother Nature in creating turbulence that is unique to a particular set of boundary conditions.<br />
<br />
'''Turbulent flows are indeed flows!'''. And that is the problem.<br />
<br />
==2.4 The Origins of Turbulence==<br />
<br />
Turbulent flows can often be observed to arise from laminar flows as the Reynolds number, (or someother relevant parameter) is increased. This happens because small disturbances to the flow are no longer damped by the flow, but begin to grow by taking energy from the original laminar flow. This natural process is easily visualized by watching the simple stream of water from a faucet (or even a pitcher). Turn the flow on very slow (or pour) so the stream is very smooth initially, at least near the outlet. Now slowly open the faucet (or pour faster) abd observe what happens, first far away, then closer to the spout. The surface begins to exhibit waves or ripples which appear to grow downstream . In fact, they are growing by extracting energy from the primary flow. Eventually they grow enough that the flow breaks into drops. These are capillary instabilities arisiing from surface tension, but regardless of the type of instability, the idea is the same -small (or infinitesimal ) disturbances have grown to disrupt the serenity (and simplicity) of laminar flow.<br />
<br />
The manner in which the instabilities grow naturally in a flow can be examined using the equations we have already developed above. We derived them by decomposing the motion into a mean and fluctuating part. But suppose instead we had decomposed the motion into a base flow part (the initial laminar part) and into a disturbance which represents a fluctuating part superimposed on the base flow. The result of substituting such a decomposition into the full Navier-Stokes equations and averaging is precisely that given by equations 2.13 and 2.15. But the very important difference is the additional restriction that what was previously identified as the mean (or averaged ) motion is now also the base or laminar flow.<br />
<br />
Now if the base flow is really laminar flow (which it must be by our original hypothesis), then our averaged equations governing the base flow must yield the same mean flow as the original laminar flow on which the disturbances was superimposed. But this can happen only if these new averaged equations reduce to '''exactly''' the same lamiane flow equations without any evidence of a disturbance. Clearly from equations 2.13 and 2.15, this can happen ''only if all the Reynolds stress terms vanish identically!'' Obviously this requires that the disturbances be infintesimal so the extra terms can be neglected - hence our interest in infinitesimal disturbances.<br />
<br />
So we hypothesized a base flow which was laminar and showed that it is unchanged even with the imposition of infintesimal disturbances on it - ''but only as long as the disturbances'' '''remain''' ''infinitesimal!'' What happens if the disturbance starts to grow? Obviously before we conclude that all laminar flows are laminar forever we better investigate whether or not these infinitesimal disturbances can grow to ''finite'' size. To do this we need an equation for the fluctuation itself.<br />
<br />
more to come soon.............<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/User:PavitranUser:Pavitran2005-09-15T11:54:33Z<p>Pavitran: </p>
<hr />
<div>I am Dynampally Pavitran, I am about to complete my Masters in Computational and Experimental Turbulence at Chalmers. Presently I am doing my thesis at SP Swedish National Testing and Research Institute.</div>Pavitranhttps://www.cfd-online.com/Wiki/CFD-Wiki:Community_portalCFD-Wiki:Community portal2005-09-15T11:48:16Z<p>Pavitran: /* Wiki editors - Who we are */</p>
<hr />
<div>Here you can find an overview over what is happening with this Wiki. For help with technical details and guidelines for how to contribute material to the Wiki read the [[Help:Contents]] page. We also have a [http://www.cfd-online.com/Forum/wiki.cgi discussion forum] for people that edit or add material to the Wiki. Both of these can also be reached via the links in the navigation section to the left.<br />
<br />
'''We should be focusing on two things right now:'''<br />
<br />
*Create a basic structure for the Wiki, with &quot;stubs&quot; (short articles which invite people to add more to them) for sections/pages where we still lack any significant content.<br />
<br />
*Finish a few high-quality sections/pages which other contributors can be inspired by and copy. This is important since we need to develop some form of common standard and understanding of how this Wiki should &quot;look and feel&quot;, how we write references, equations, links etc. I have started to summarize a few guidelines [[Help:Contents | here]].<br />
<br />
The intention is that we should launch the Wiki publicly by the end of October (which year depends on us ;-)<br />
<br />
== What's in the works ==<br />
<br />
You who do significant additions to the Wiki, please add some information about your work, plans and progress here so that others can see what you are working on and perhaps help, monitor, come with suggestions and most importantly, be inspired by.<br />
<br />
* I've started work on creating a first best-practise guide. I think we need some sort of first example for how one of these should look. Writing this type of guide is not easy though. I chose to start with a guide for [[Best practise guidelines for turbomachinery CFD | turbomachinery CFD]] since that is the application area that I know most about. So far I've only created a table of contents and written a few sections. What do you think about it? --[[User:Jola|Jola]] 08:09, 8 September 2005 (MDT)<br />
<br />
* I have written a fairly comprehensive description of a turbulence model - see the [[Baldwin-Lomax model]] page. I think that this is a fairly good example of how we should write this type of page about a certain model. Or what do you think? --[[User:Jola|Jola]] 08:09, 8 September 2005 (MDT)<br />
<br />
* A few people have started work on the numerics section. It now has a basic structure and I think it will get significantly more content within the next few days. --[[User:Jola|Jola]] 17:02, 12 September 2005 (MDT)<br />
<br />
* We have been allowed to base the turbulence section on an excellent intro book on turbulence written by Professor William K. George. So expect to see significant additions in this section soon. --[[User:Jola|Jola]] 17:02, 12 September 2005 (MDT)<br />
<br />
* We need to coordinate our style and structure. I've started to write some guidlines on the Help page and have also emailed a few of you trying to get things a bit more organized. --[[User:Jola|Jola]] 17:02, 12 September 2005 (MDT)<br />
<br />
* I'm trying to engage editors to write FAQ's for the largest codes (Fluent, CFX and STAR-CD) and have posted messages to the code forums here at CFD Online inviting people to contribute to these FAQ's. I'm also in direct contact with a couple of support managers at CFD companies who are wondering how they can or if they should support this FAQ initiative. We'll see what happens. --[[User:Jola|Jola]] 17:02, 12 September 2005 (MDT)<br />
<br />
== What needs to be done ==<br />
<br />
* We have many turbulence models listed in the [[turbulence modeling]] section which still lack any description. Feel free to pick a model that you are familiar with and write a description of it. --[[User:Jola|Jola]] 01:50, 13 September 2005 (MDT)<br />
<br />
* The [[FAQ's | FAQ]] section is still very thin. If you are familiar with one of the larger CFD codes please consider adding a few questions and answers to the FAQ. --[[User:Jola|Jola]] 08:28, 13 September 2005 (MDT)<br />
<br />
== Wiki editors - Who we are ==<br />
<br />
Add your name here if you make contributions to the wiki.<br />
<br />
* [[User:praveen]] - Praveen. C<br />
* [[User:Michail]] - Michail Kirichkov<br />
* [[User:jola]] - Jonas Larsson<br />
* [[User:zxaar]] - Arjun Yadav<br />
* [[User:Pavitran]] - Pavitran. D</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-15T11:44:41Z<p>Pavitran: /* 2. Equations for the Average Velocity */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a so called [[DNS simulation]] (Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''2.1 The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2.2 Equations for the Average Velocity==<br />
<br />
<br />
Although laminar solutions to the equations often exist that are consistent with the boundary conditions, perturbations to these solutions(sometimes even infinitesimal) can cause them to become turbulent. To see how this can happen, it is convenient to analyze the flow in two parts, a mean(or average) component and a fluctuating component. Thus the instantaneous velocity and stresses can be written as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\tilde {u_i}=U_i+u_i<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde p=P+p<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde T_{ij}^{(v)}=T_{ij}^{(v)}+\tau_{ij}^{(v)}<br />
&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.11)&lt;/td&gt;&lt;/tr&gt;<br />
&lt;/table&gt;<br />
<br />
<br />
Where &lt;math&gt;U_i&lt;/math&gt;, &lt;math&gt;P&lt;/math&gt; and &lt;math&gt;T_{ij}^{(v)}&lt;/math&gt; represent the mean motion, and &lt;math&gt;u_i&lt;/math&gt;, &lt;math&gt;p&lt;/math&gt; and &lt;math&gt;\tau_{ij}^{(v)}&lt;/math&gt; the fluctuating motions. This technique for decomposing the instantaneous motion is referred to as the '''''Reynolds decomposition.''''' Note that if the averages are defined as ensemble means, they are, in general, time-dependent. For the remainder of this material unless other wise stated, the density will be assumed constant so&lt;math&gt;\tilde{\rho}\equiv\rho&lt;/math&gt;,and its fluctuation is zero.<br />
<br />
Substitution of equations 2.11 into equations 2.10 yields<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial (U_i+u_i)}{\partial t}+(U_j+u_j)\frac{\partial (U_i+u_i)}{\partial x_j}\right] = -\frac{\partial (P+p)}{\partial x_i}+\frac{\partial (T_{ij}^{(v)}+\tau_{ij}^{(v)})}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.12)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This equation can now be averaged to yield an equation expressing momentum conservation for the averaged motion. Note that the operations of averaging and differentiation commute; i.e., the average of a derivative is the same as the derivative of the average. Also the average of a fluctuating quantity is zero. Thus the equation for the averaged motion reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\rho\left \langle u_j\frac{\partial u_i }{\partial x_j} \right \rangle&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.13)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the remaining fluctuating product term has been moved to the right hand side of the equation. Whether or not the last term is zero like the other fluctuating term depends on the correlation of the terms in the product. In general, these correlations are not zero.<br />
<br />
The mass conservation equation can be similarly decomposed. In incompressible form, substitution of equations 2.11 into equation 2.4 yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial (U_j+u_j)}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.14)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
of which average is <br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial U_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.15)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
It is clear from the equation 2.15 that the averaged motion satisfies the same form of the mass conservation equation as does the instantaneous motion at least for incompressible flows. How much simpler the turbulence problem would be if the same were true for the momentum! Unfortunately, as is easily seen from the equation 2.13, such is not the case.<br />
<br />
Equation 2.15 can be subtracted from equation 2.14 to yield an equation for instantaneous motion alone; i.e,<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial u_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.16)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Again, like the mean, the form of the original instantaneous equation is seen to be preserved. The reason, of course, is obvious: the continuity equation is linear. The momentum equation , on the other hand, is not; hence the difference.<br />
<br />
Equation 2.16 can be used to rewrite the last term in equation 2.13 for the mean momentum. Multiplying equation 2.16 by &lt;math&gt;u_i&lt;/math&gt; and averaging yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.17)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This can be added to:&lt;math&gt;\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle&lt;/math&gt; to obtain:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle +0=\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle+ \left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle =\frac{ \partial}{\partial x_j}{\left \langle u_iu_j\right \rangle} &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.18)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Where again the fact that arithmetic and averaging operations commute has been used.<br />
<br />
The equation for the averaged momentum, equation 2.13 can now be rewritten as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\frac{ \partial}{\partial x_j}{(\rho\left \langle u_iu_j\right \rangle)}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.19)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The last two terms on the right hand side are both divergence terms and can be combined; the result is:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial }{\partial x_j}[T_{ij}^{(v)}-{\rho\left \langle u_iu_j\right \rangle}]&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.20)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Now the terms in square brackets on the right have the dimensions of stress. The first term is, in fact , the viscous stress. The second term, on the other hand, is not a stress at all, but simply a re-worked version of the fluctuating contribution to the non-linear acceleration terms. The fact that it can be written this way, however, indicates that at least as far as the motion is concerned, it acts as though it were a stress- hence its name, the '''Reynolds stress'''. In the succeeding sections the consequences of this difference will be examined.<br />
<br />
==2.3 The Turbulence Problem==<br />
<br />
It is the appearance of the Reynolds stress which makes the turbulence problem so difficult - at least from the engineers perspective. Even though we can pretend it is a stress, the physics which give rise to it are very different from the viscous stress. The viscous stress can be related directly to the other flow properties by constitutive equations, which in turn depend only on the properties of the fluid( as in equation 2.5 for a Newtonian fluid). The reason this works is that when we make such closure approximations for a fluid, we are averaging over characteristic length and time scales much smaller than those of the flows we are interested in. Yet at the same time, these scales are much larger than the molecular length and time scales which characterize the molecular interactions that are actually causing the momentum transfer. (This is what the continuum approximation is all about).<br />
<br />
The '''''Reynolds stress''''', on the other hand, arises from the flow itself! Worse, the scales of the fluctuating motion which give rise to it are the scales we are interested in. This means that the closure ideas which worked so well for the viscous stress, should not be expected to work too well for the Reynolds stress. And as we shall see, they do not.<br />
<br />
This leaves us in a terrible position. Physics and engineering are all about writing equations(and boundary conditions) so we can solve them to make predictions. We don't want to have a build prototype airplanes first to see if they will they fall out of the sky. Instead we want to be able to analyze our designs before building, to save the cost in money and lives if our ideas are wrong. The same is true for dams and bridges and tunnels and automobiles. If we had confidence in our turbulence models, we could even build huge one-offs and expect them to work the first time. Unfortunately, even though turbulence models have improved to the point where we can use them in design, we still cannot trust them enough to eliminate expensive wind tunnel and model studies. And recent history is full of examples to prove this.<br />
<br />
The turbulence problem (from the engineers perspective) is then three-fold:<br />
<br />
* '''The averaged equations are not closed.''' Count the number of unknowns in equation 2.20 above. Then count the number of equations. Even with the continuity equation we have atleast six equations too few.<br />
<br />
* '''The simple ideas to provide the extra equations usually do not work.''' And even when we can fix them up for a particular class of flows (like the flow in a pipe, for example), they will most likely not be able to predict what happens in even a slightly different environment(like a bend).<br />
<br />
*'''Even the last resort of compiling engineering tables for design handbooks carries substantial risk.''' This is the last resort for the engineer who lacks equations or cannot trust them. Even when based on a wealth of experience, they require expensive model testing to see if they can be extrapolated to a particular situation. Often they cannot, so infinitely clever is Mother Nature in creating turbulence that is unique to a particular set of boundary conditions.<br />
<br />
'''Turbulent flows are indeed flows!'''. And that is the problem.<br />
<br />
==2.4 The Origins of Turbulence==<br />
<br />
will be updated soon.................<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-15T11:44:08Z<p>Pavitran: /* The Reynolds Averaged Equations and the Turbulence Closure Problem */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a so called [[DNS simulation]] (Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''2.1 The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2. Equations for the Average Velocity==<br />
<br />
<br />
Although laminar solutions to the equations often exist that are consistent with the boundary conditions, perturbations to these solutions(sometimes even infinitesimal) can cause them to become turbulent. To see how this can happen, it is convenient to analyze the flow in two parts, a mean(or average) component and a fluctuating component. Thus the instantaneous velocity and stresses can be written as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\tilde {u_i}=U_i+u_i<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde p=P+p<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde T_{ij}^{(v)}=T_{ij}^{(v)}+\tau_{ij}^{(v)}<br />
&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.11)&lt;/td&gt;&lt;/tr&gt;<br />
&lt;/table&gt;<br />
<br />
<br />
Where &lt;math&gt;U_i&lt;/math&gt;, &lt;math&gt;P&lt;/math&gt; and &lt;math&gt;T_{ij}^{(v)}&lt;/math&gt; represent the mean motion, and &lt;math&gt;u_i&lt;/math&gt;, &lt;math&gt;p&lt;/math&gt; and &lt;math&gt;\tau_{ij}^{(v)}&lt;/math&gt; the fluctuating motions. This technique for decomposing the instantaneous motion is referred to as the '''''Reynolds decomposition.''''' Note that if the averages are defined as ensemble means, they are, in general, time-dependent. For the remainder of this material unless other wise stated, the density will be assumed constant so&lt;math&gt;\tilde{\rho}\equiv\rho&lt;/math&gt;,and its fluctuation is zero.<br />
<br />
Substitution of equations 2.11 into equations 2.10 yields<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial (U_i+u_i)}{\partial t}+(U_j+u_j)\frac{\partial (U_i+u_i)}{\partial x_j}\right] = -\frac{\partial (P+p)}{\partial x_i}+\frac{\partial (T_{ij}^{(v)}+\tau_{ij}^{(v)})}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.12)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This equation can now be averaged to yield an equation expressing momentum conservation for the averaged motion. Note that the operations of averaging and differentiation commute; i.e., the average of a derivative is the same as the derivative of the average. Also the average of a fluctuating quantity is zero. Thus the equation for the averaged motion reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\rho\left \langle u_j\frac{\partial u_i }{\partial x_j} \right \rangle&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.13)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the remaining fluctuating product term has been moved to the right hand side of the equation. Whether or not the last term is zero like the other fluctuating term depends on the correlation of the terms in the product. In general, these correlations are not zero.<br />
<br />
The mass conservation equation can be similarly decomposed. In incompressible form, substitution of equations 2.11 into equation 2.4 yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial (U_j+u_j)}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.14)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
of which average is <br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial U_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.15)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
It is clear from the equation 2.15 that the averaged motion satisfies the same form of the mass conservation equation as does the instantaneous motion at least for incompressible flows. How much simpler the turbulence problem would be if the same were true for the momentum! Unfortunately, as is easily seen from the equation 2.13, such is not the case.<br />
<br />
Equation 2.15 can be subtracted from equation 2.14 to yield an equation for instantaneous motion alone; i.e,<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial u_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.16)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Again, like the mean, the form of the original instantaneous equation is seen to be preserved. The reason, of course, is obvious: the continuity equation is linear. The momentum equation , on the other hand, is not; hence the difference.<br />
<br />
Equation 2.16 can be used to rewrite the last term in equation 2.13 for the mean momentum. Multiplying equation 2.16 by &lt;math&gt;u_i&lt;/math&gt; and averaging yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.17)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This can be added to:&lt;math&gt;\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle&lt;/math&gt; to obtain:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle +0=\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle+ \left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle =\frac{ \partial}{\partial x_j}{\left \langle u_iu_j\right \rangle} &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.18)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Where again the fact that arithmetic and averaging operations commute has been used.<br />
<br />
The equation for the averaged momentum, equation 2.13 can now be rewritten as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\frac{ \partial}{\partial x_j}{(\rho\left \langle u_iu_j\right \rangle)}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.19)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The last two terms on the right hand side are both divergence terms and can be combined; the result is:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial }{\partial x_j}[T_{ij}^{(v)}-{\rho\left \langle u_iu_j\right \rangle}]&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.20)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Now the terms in square brackets on the right have the dimensions of stress. The first term is, in fact , the viscous stress. The second term, on the other hand, is not a stress at all, but simply a re-worked version of the fluctuating contribution to the non-linear acceleration terms. The fact that it can be written this way, however, indicates that at least as far as the motion is concerned, it acts as though it were a stress- hence its name, the '''Reynolds stress'''. In the succeeding sections the consequences of this difference will be examined.<br />
<br />
==2.3 The Turbulence Problem==<br />
<br />
It is the appearance of the Reynolds stress which makes the turbulence problem so difficult - at least from the engineers perspective. Even though we can pretend it is a stress, the physics which give rise to it are very different from the viscous stress. The viscous stress can be related directly to the other flow properties by constitutive equations, which in turn depend only on the properties of the fluid( as in equation 2.5 for a Newtonian fluid). The reason this works is that when we make such closure approximations for a fluid, we are averaging over characteristic length and time scales much smaller than those of the flows we are interested in. Yet at the same time, these scales are much larger than the molecular length and time scales which characterize the molecular interactions that are actually causing the momentum transfer. (This is what the continuum approximation is all about).<br />
<br />
The '''''Reynolds stress''''', on the other hand, arises from the flow itself! Worse, the scales of the fluctuating motion which give rise to it are the scales we are interested in. This means that the closure ideas which worked so well for the viscous stress, should not be expected to work too well for the Reynolds stress. And as we shall see, they do not.<br />
<br />
This leaves us in a terrible position. Physics and engineering are all about writing equations(and boundary conditions) so we can solve them to make predictions. We don't want to have a build prototype airplanes first to see if they will they fall out of the sky. Instead we want to be able to analyze our designs before building, to save the cost in money and lives if our ideas are wrong. The same is true for dams and bridges and tunnels and automobiles. If we had confidence in our turbulence models, we could even build huge one-offs and expect them to work the first time. Unfortunately, even though turbulence models have improved to the point where we can use them in design, we still cannot trust them enough to eliminate expensive wind tunnel and model studies. And recent history is full of examples to prove this.<br />
<br />
The turbulence problem (from the engineers perspective) is then three-fold:<br />
<br />
* '''The averaged equations are not closed.''' Count the number of unknowns in equation 2.20 above. Then count the number of equations. Even with the continuity equation we have atleast six equations too few.<br />
<br />
* '''The simple ideas to provide the extra equations usually do not work.''' And even when we can fix them up for a particular class of flows (like the flow in a pipe, for example), they will most likely not be able to predict what happens in even a slightly different environment(like a bend).<br />
<br />
*'''Even the last resort of compiling engineering tables for design handbooks carries substantial risk.''' This is the last resort for the engineer who lacks equations or cannot trust them. Even when based on a wealth of experience, they require expensive model testing to see if they can be extrapolated to a particular situation. Often they cannot, so infinitely clever is Mother Nature in creating turbulence that is unique to a particular set of boundary conditions.<br />
<br />
'''Turbulent flows are indeed flows!'''. And that is the problem.<br />
<br />
==2.4 The Origins of Turbulence==<br />
<br />
will be updated soon.................<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-15T11:43:09Z<p>Pavitran: /* 2.3 The Turbulence Problem */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a so called [[DNS simulation]] (Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''1. The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2. Equations for the Average Velocity==<br />
<br />
<br />
Although laminar solutions to the equations often exist that are consistent with the boundary conditions, perturbations to these solutions(sometimes even infinitesimal) can cause them to become turbulent. To see how this can happen, it is convenient to analyze the flow in two parts, a mean(or average) component and a fluctuating component. Thus the instantaneous velocity and stresses can be written as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\tilde {u_i}=U_i+u_i<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde p=P+p<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde T_{ij}^{(v)}=T_{ij}^{(v)}+\tau_{ij}^{(v)}<br />
&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.11)&lt;/td&gt;&lt;/tr&gt;<br />
&lt;/table&gt;<br />
<br />
<br />
Where &lt;math&gt;U_i&lt;/math&gt;, &lt;math&gt;P&lt;/math&gt; and &lt;math&gt;T_{ij}^{(v)}&lt;/math&gt; represent the mean motion, and &lt;math&gt;u_i&lt;/math&gt;, &lt;math&gt;p&lt;/math&gt; and &lt;math&gt;\tau_{ij}^{(v)}&lt;/math&gt; the fluctuating motions. This technique for decomposing the instantaneous motion is referred to as the '''''Reynolds decomposition.''''' Note that if the averages are defined as ensemble means, they are, in general, time-dependent. For the remainder of this material unless other wise stated, the density will be assumed constant so&lt;math&gt;\tilde{\rho}\equiv\rho&lt;/math&gt;,and its fluctuation is zero.<br />
<br />
Substitution of equations 2.11 into equations 2.10 yields<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial (U_i+u_i)}{\partial t}+(U_j+u_j)\frac{\partial (U_i+u_i)}{\partial x_j}\right] = -\frac{\partial (P+p)}{\partial x_i}+\frac{\partial (T_{ij}^{(v)}+\tau_{ij}^{(v)})}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.12)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This equation can now be averaged to yield an equation expressing momentum conservation for the averaged motion. Note that the operations of averaging and differentiation commute; i.e., the average of a derivative is the same as the derivative of the average. Also the average of a fluctuating quantity is zero. Thus the equation for the averaged motion reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\rho\left \langle u_j\frac{\partial u_i }{\partial x_j} \right \rangle&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.13)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the remaining fluctuating product term has been moved to the right hand side of the equation. Whether or not the last term is zero like the other fluctuating term depends on the correlation of the terms in the product. In general, these correlations are not zero.<br />
<br />
The mass conservation equation can be similarly decomposed. In incompressible form, substitution of equations 2.11 into equation 2.4 yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial (U_j+u_j)}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.14)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
of which average is <br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial U_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.15)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
It is clear from the equation 2.15 that the averaged motion satisfies the same form of the mass conservation equation as does the instantaneous motion at least for incompressible flows. How much simpler the turbulence problem would be if the same were true for the momentum! Unfortunately, as is easily seen from the equation 2.13, such is not the case.<br />
<br />
Equation 2.15 can be subtracted from equation 2.14 to yield an equation for instantaneous motion alone; i.e,<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial u_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.16)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Again, like the mean, the form of the original instantaneous equation is seen to be preserved. The reason, of course, is obvious: the continuity equation is linear. The momentum equation , on the other hand, is not; hence the difference.<br />
<br />
Equation 2.16 can be used to rewrite the last term in equation 2.13 for the mean momentum. Multiplying equation 2.16 by &lt;math&gt;u_i&lt;/math&gt; and averaging yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.17)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This can be added to:&lt;math&gt;\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle&lt;/math&gt; to obtain:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle +0=\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle+ \left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle =\frac{ \partial}{\partial x_j}{\left \langle u_iu_j\right \rangle} &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.18)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Where again the fact that arithmetic and averaging operations commute has been used.<br />
<br />
The equation for the averaged momentum, equation 2.13 can now be rewritten as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\frac{ \partial}{\partial x_j}{(\rho\left \langle u_iu_j\right \rangle)}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.19)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The last two terms on the right hand side are both divergence terms and can be combined; the result is:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial }{\partial x_j}[T_{ij}^{(v)}-{\rho\left \langle u_iu_j\right \rangle}]&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.20)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Now the terms in square brackets on the right have the dimensions of stress. The first term is, in fact , the viscous stress. The second term, on the other hand, is not a stress at all, but simply a re-worked version of the fluctuating contribution to the non-linear acceleration terms. The fact that it can be written this way, however, indicates that at least as far as the motion is concerned, it acts as though it were a stress- hence its name, the '''Reynolds stress'''. In the succeeding sections the consequences of this difference will be examined.<br />
<br />
==2.3 The Turbulence Problem==<br />
<br />
It is the appearance of the Reynolds stress which makes the turbulence problem so difficult - at least from the engineers perspective. Even though we can pretend it is a stress, the physics which give rise to it are very different from the viscous stress. The viscous stress can be related directly to the other flow properties by constitutive equations, which in turn depend only on the properties of the fluid( as in equation 2.5 for a Newtonian fluid). The reason this works is that when we make such closure approximations for a fluid, we are averaging over characteristic length and time scales much smaller than those of the flows we are interested in. Yet at the same time, these scales are much larger than the molecular length and time scales which characterize the molecular interactions that are actually causing the momentum transfer. (This is what the continuum approximation is all about).<br />
<br />
The '''''Reynolds stress''''', on the other hand, arises from the flow itself! Worse, the scales of the fluctuating motion which give rise to it are the scales we are interested in. This means that the closure ideas which worked so well for the viscous stress, should not be expected to work too well for the Reynolds stress. And as we shall see, they do not.<br />
<br />
This leaves us in a terrible position. Physics and engineering are all about writing equations(and boundary conditions) so we can solve them to make predictions. We don't want to have a build prototype airplanes first to see if they will they fall out of the sky. Instead we want to be able to analyze our designs before building, to save the cost in money and lives if our ideas are wrong. The same is true for dams and bridges and tunnels and automobiles. If we had confidence in our turbulence models, we could even build huge one-offs and expect them to work the first time. Unfortunately, even though turbulence models have improved to the point where we can use them in design, we still cannot trust them enough to eliminate expensive wind tunnel and model studies. And recent history is full of examples to prove this.<br />
<br />
The turbulence problem (from the engineers perspective) is then three-fold:<br />
<br />
* '''The averaged equations are not closed.''' Count the number of unknowns in equation 2.20 above. Then count the number of equations. Even with the continuity equation we have atleast six equations too few.<br />
<br />
* '''The simple ideas to provide the extra equations usually do not work.''' And even when we can fix them up for a particular class of flows (like the flow in a pipe, for example), they will most likely not be able to predict what happens in even a slightly different environment(like a bend).<br />
<br />
*'''Even the last resort of compiling engineering tables for design handbooks carries substantial risk.''' This is the last resort for the engineer who lacks equations or cannot trust them. Even when based on a wealth of experience, they require expensive model testing to see if they can be extrapolated to a particular situation. Often they cannot, so infinitely clever is Mother Nature in creating turbulence that is unique to a particular set of boundary conditions.<br />
<br />
'''Turbulent flows are indeed flows!'''. And that is the problem.<br />
<br />
==2.4 The Origins of Turbulence==<br />
<br />
will be updated soon.................<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-15T09:04:47Z<p>Pavitran: /* 2. Equations for the Average Velocity */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a so called [[DNS simulation]] (Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''1. The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2. Equations for the Average Velocity==<br />
<br />
<br />
Although laminar solutions to the equations often exist that are consistent with the boundary conditions, perturbations to these solutions(sometimes even infinitesimal) can cause them to become turbulent. To see how this can happen, it is convenient to analyze the flow in two parts, a mean(or average) component and a fluctuating component. Thus the instantaneous velocity and stresses can be written as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\tilde {u_i}=U_i+u_i<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde p=P+p<br />
&lt;/math&gt;<br />
:&lt;math&gt;<br />
\tilde T_{ij}^{(v)}=T_{ij}^{(v)}+\tau_{ij}^{(v)}<br />
&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.11)&lt;/td&gt;&lt;/tr&gt;<br />
&lt;/table&gt;<br />
<br />
<br />
Where &lt;math&gt;U_i&lt;/math&gt;, &lt;math&gt;P&lt;/math&gt; and &lt;math&gt;T_{ij}^{(v)}&lt;/math&gt; represent the mean motion, and &lt;math&gt;u_i&lt;/math&gt;, &lt;math&gt;p&lt;/math&gt; and &lt;math&gt;\tau_{ij}^{(v)}&lt;/math&gt; the fluctuating motions. This technique for decomposing the instantaneous motion is referred to as the '''''Reynolds decomposition.''''' Note that if the averages are defined as ensemble means, they are, in general, time-dependent. For the remainder of this material unless other wise stated, the density will be assumed constant so&lt;math&gt;\tilde{\rho}\equiv\rho&lt;/math&gt;,and its fluctuation is zero.<br />
<br />
Substitution of equations 2.11 into equations 2.10 yields<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial (U_i+u_i)}{\partial t}+(U_j+u_j)\frac{\partial (U_i+u_i)}{\partial x_j}\right] = -\frac{\partial (P+p)}{\partial x_i}+\frac{\partial (T_{ij}^{(v)}+\tau_{ij}^{(v)})}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.12)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This equation can now be averaged to yield an equation expressing momentum conservation for the averaged motion. Note that the operations of averaging and differentiation commute; i.e., the average of a derivative is the same as the derivative of the average. Also the average of a fluctuating quantity is zero. Thus the equation for the averaged motion reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\rho\left \langle u_j\frac{\partial u_i }{\partial x_j} \right \rangle&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.13)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the remaining fluctuating product term has been moved to the right hand side of the equation. Whether or not the last term is zero like the other fluctuating term depends on the correlation of the terms in the product. In general, these correlations are not zero.<br />
<br />
The mass conservation equation can be similarly decomposed. In incompressible form, substitution of equations 2.11 into equation 2.4 yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial (U_j+u_j)}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.14)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
of which average is <br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial U_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.15)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
It is clear from the equation 2.15 that the averaged motion satisfies the same form of the mass conservation equation as does the instantaneous motion at least for incompressible flows. How much simpler the turbulence problem would be if the same were true for the momentum! Unfortunately, as is easily seen from the equation 2.13, such is not the case.<br />
<br />
Equation 2.15 can be subtracted from equation 2.14 to yield an equation for instantaneous motion alone; i.e,<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial u_j}{\partial x_j}=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.16)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Again, like the mean, the form of the original instantaneous equation is seen to be preserved. The reason, of course, is obvious: the continuity equation is linear. The momentum equation , on the other hand, is not; hence the difference.<br />
<br />
Equation 2.16 can be used to rewrite the last term in equation 2.13 for the mean momentum. Multiplying equation 2.16 by &lt;math&gt;u_i&lt;/math&gt; and averaging yields:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle=0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.17)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This can be added to:&lt;math&gt;\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle&lt;/math&gt; to obtain:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle +0=\left \langle u_j\frac{\partial u_i}{\partial x_j}\right \rangle+ \left \langle u_i\frac{\partial u_j}{\partial x_j}\right \rangle =\frac{ \partial}{\partial x_j}{\left \langle u_iu_j\right \rangle} &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.18)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Where again the fact that arithmetic and averaging operations commute has been used.<br />
<br />
The equation for the averaged momentum, equation 2.13 can now be rewritten as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial T_{ij}^{(v)}}{\partial x_j}-\frac{ \partial}{\partial x_j}{(\rho\left \langle u_iu_j\right \rangle)}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.19)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The last two terms on the right hand side are both divergence terms and can be combined; the result is:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial U_i}{\partial t}+U_j\frac{\partial U_i}{\partial x_j}\right] = -\frac{\partial P}{\partial x_i}+\frac{\partial }{\partial x_j}[T_{ij}^{(v)}-{\rho\left \langle u_iu_j\right \rangle}]&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.20)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Now the terms in square brackets on the right have the dimensions of stress. The first term is, in fact , the viscous stress. The second term, on the other hand, is not a stress at all, but simply a re-worked version of the fluctuating contribution to the non-linear acceleration terms. The fact that it can be written this way, however, indicates that at least as far as the motion is concerned, it acts as though it were a stress- hence its name, the '''Reynolds stress'''. In the succeeding sections the consequences of this difference will be examined.<br />
<br />
==2.3 The Turbulence Problem==<br />
will be updated soon........<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/Talk:Introduction_to_turbulenceTalk:Introduction to turbulence2005-09-14T09:13:53Z<p>Pavitran: /* 2. Equations for the Average Velocity */</p>
<hr />
<div></div>Pavitranhttps://www.cfd-online.com/Wiki/Talk:Introduction_to_turbulenceTalk:Introduction to turbulence2005-09-14T09:13:43Z<p>Pavitran: /* The Reynolds Averaged Equations and the Turbulence Closure Problem */</p>
<hr />
<div><br />
<br />
==2. Equations for the Average Velocity==</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-14T09:13:02Z<p>Pavitran: /* What is Turbulence? */</p>
<hr />
<div>==What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional [[vorticity]]. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
===Why Study Turbulence?===<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous [[turbulence closure problem]].<br />
<br />
Of course, closure is not a problem when performing a so called [[DNS simulation]] (Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
===The cost of our ignorance===<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
===What do we really know for sure?===<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone (Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''1. The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2. Equations for the Average Velocity==<br />
<br />
==Credits==<br />
'''This text was based on &quot;Introduction to Turbulence&quot; by Professor William K.George, Chalmers University of Technology, Sweden.'''</div>Pavitranhttps://www.cfd-online.com/Wiki/Talk:Introduction_to_turbulenceTalk:Introduction to turbulence2005-09-14T09:09:32Z<p>Pavitran: </p>
<hr />
<div>==The Reynolds Averaged Equations and the Turbulence Closure Problem==<br />
<br />
'''1. The Equations Governing the Instantaneous Fluid Motions'''<br />
<br />
All fluid motions, whether turbulent or not, are governed by the dynamical equations for a fluid. These can be written using Cartesian tensor notation as:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.1)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}\right]+ \tilde{\rho}\frac{\partial \tilde{u_j}}{\partial x_j}= 0 &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.2)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
<br />
where &lt;math&gt;\tilde{u_i}(\vec{x},t)&lt;/math&gt; represents the i-the component of the fluid velocity at a point in space,&lt;math&gt;[\vec{x}]_i=x_i&lt;/math&gt;, and time,t. Also <br />
&lt;math&gt;\tilde{p}(\vec{x},t)&lt;/math&gt; represents the static pressure, &lt;math&gt;\tilde{T_{ij}}^{(v)}(\vec{x},t)&lt;/math&gt;, the viscous(or deviatoric) stresses, and &lt;math&gt;\tilde\rho&lt;/math&gt; the fluid density. The tilde over the symbol indicates that an instantaneous quantity is being considered. Also the Einstein summation convention has been employed[1].<br />
<br />
In equation 2.1, the subscript i is a free index which can take on the values 1,2 and 3. Thus equation 2.1 is in reality three separate equations. These three equations are just Newton's second law written for a continuum in a spatial(or Eulerian) reference frame. Together they relate the rate of change of momentum per unit mass &lt;math&gt;(\rho{u_i})&lt;/math&gt;,a vector quantity, to the contact and body forces.<br />
<br />
Equation 2.2 is the equation for mass conservation in the absence of sources(or sinks) of mass. Almost all flows considered in this material will be incompressible, which implies that derivative of the density following the fluid material[the term in brackets] is zero. Thus for incompressible flows, the mass conservation equation reduces to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{D \tilde{\rho}}{Dt}=\frac{\partial \tilde{\rho}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{\rho}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.3)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From equation 2.2 it follows that for incompressible flows,<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\frac{\partial \tilde{u_i}}{\partial x_j}= 0&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.4)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscous stresses(the stress minus the mean normal stress) are represented by the tensor&lt;math&gt;\tilde{T_{ij}}^{(v)}&lt;/math&gt;. From its definition,&lt;math&gt;\tilde{T_{kk}}^{(v)}&lt;/math&gt;=0. In many flows of interest, the fluid behaves as a Newtonian fluid in which the viscous stress can be related to the fluid motion by a constitutive relation of the form.<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu[\tilde{s_{ij}}-\frac{1}{3}\tilde{s_{kk}}\delta_{ij}] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.5)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
The viscosity, &lt;math&gt;\mu&lt;/math&gt;, is a property of the fluid that can be measured in an independent experiment. &lt;math&gt;\tilde s_{ij}&lt;/math&gt; is the instantaneous strain rate tensor defoned by<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{s_{ij}}= \frac{1}{2}\left[\frac{\partial \tilde u_i}{\partial x_j}+\frac{\partial \tilde u_j}{\partial x_i}\right] &lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.6)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
From its definition, &lt;math&gt;\tilde s_{kk}=\frac{\partial \tilde u_k}{\partial x_k}&lt;/math&gt;. If the flow is incompressible, &lt;math&gt;\tilde s_{kk}=0&lt;/math&gt; and the Newtonian constitutive equation reduces to<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\tilde{T_{ij}}^{(v)}= 2\mu\tilde{s_{ij}}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.7)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Throughout this material, unless explicitly stated otherwise, the density &lt;math&gt;\tilde\rho=\rho&lt;/math&gt; and the viscosity &lt;math&gt;\mu&lt;/math&gt; will be assumed constant. With these assumptions, the instantaneous momentum equations for a Newtonian Fluid reduce to:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac {1}{\tilde\rho}\frac{\partial \tilde{p}}{\partial x_i}+\nu\frac{\partial^2 {\tilde{u_i}}}{\partial x_j^2}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.8)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
where the kinematic viscosity, &lt;math&gt;\nu&lt;/math&gt;, has been defined as:<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
&lt;math&gt;\nu\equiv\frac{\mu}{\rho}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.9)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
Note that since the density is assumed conastant, the tilde is no longer necessary.<br />
<br />
Sometimes it will be more instructive and convenient to not explicitly include incompressibilty in the stress term, but to refer to the incompressible momentum equation in the following form:<br />
<br />
&lt;table width=&quot;100%&quot;&gt;<br />
&lt;tr&gt;&lt;td&gt;<br />
:&lt;math&gt;<br />
\rho\left[\frac{\partial \tilde{u_i}}{\partial t}+\tilde{u_j}\frac{\partial \tilde{u_i}}{\partial x_j}\right] = -\frac{\partial \tilde{p}}{\partial x_i}+\frac{\partial \tilde{T_{ij}}^{(v)}}{\partial x_j}&lt;/math&gt;<br />
&lt;/td&gt;&lt;td width=&quot;5%&quot;&gt;(2.10)&lt;/td&gt;&lt;/tr&gt;&lt;/table&gt;<br />
<br />
This form has the advantage that it is easier to keep track of the exact role of the viscous stresses.<br />
<br />
==2. Equations for the Average Velocity==</div>Pavitranhttps://www.cfd-online.com/Wiki/Introduction_to_turbulenceIntroduction to turbulence2005-09-07T07:57:10Z<p>Pavitran: </p>
<hr />
<div>==1 What is Turbulence?==<br />
<br />
Turbulence is that state of fluid motion which is characterized by apparently random and chaotic three-dimensional vorticity. When turbulence is present, it usually dominates all other flow phenomena and results in increased energy dissipation, mixing, heat transfer, and drag.<br />
<br />
For a long time scientists were not really sure in which sense turbulence is 'random', but they were pretty sure it was. Like any one who is trained in physics, we believe the flows we see around us must be the solution to some set of equations which govern. (This is after all what mechanics is about- writing equations to describe and predict the world around us) But because of the nature of the turbulence, it wasn't clear whether the equations themselves had some hidden randomness, or just the solutions. And if the latter, was it something the equations did to them, or a consequence of the intial conditions<br />
<br />
==1.2 Why Study Turbulence?==<br />
<br />
There really are the two reasons for studying turbulence- engineering and physics! And they are not necessarily complementary, atleast in the short run.<br />
<br />
Certainly a case can be made that we don't know enough about the turbulence to even start to consider engineering problems. To begin with, we always have fewer equations that unknowns in any attempt to predict anything other than the instantaneous motions. This is the famous '''turbulence closure problem'''.<br />
<br />
Of course, closure is not a problem with the so-called DNS(Direct Numerical Simulations) in which we numerically produce the instantaneous motions in a computer using the exact equations governing the fluid. Unfortunately we won't be able to perform such simulations for real engineering problems until atleast a few hundred generations of computers have come and gone. And this wo't really help us too much, since even when we now perform a DNS simulation of a really simple flow, we are already overwhelmed by the amount of data and its apparent random behaviour. This is because without some kind of theory, we have no criteria for selecting from it in a single lifetime what is important.<br />
<br />
The engineer's counter argument to the scientist's lament above is:<br />
<br />
* airplanes must fly,<br />
* weather must be forecast,<br />
* sewage and water management systems must be built<br />
* society needs ever more energy-efficient hardware and gadgets.<br />
<br />
Thus the engineer argues, no matter the inadequate state of our knowledge, '''we have the responsibilty as engineers to do the best we can with what we have'''. Who, considering the needs, could seriously argue with this? Almost incredibly - some physicists do!<br />
<br />
It seems evindent then that there must be at least two levels of assault on turbulence. At one level, the very nature of turbulence must be explored. At the other level, our current state of knowledge- however inadequate it might be- must be stretched to provide engineering solutions to real problems. <br />
<br />
==1.3 The cost of our ignorance==<br />
<br />
It is difficult to place a price tag on the cost of our limited understanding of turbulence, but it requires no imagination at all to realize that it must be enormous. Try to estimate, for example, the aggregate cost to society of our limited turbulence prediction abilities which result in inadequate weather-forecasts alone. Or try to place a value on the increased cost to the consumer need of the designer of virtually every fluid-thermal system-from heat exchangers to hypersonic planes- to depend on empiricism and experimentation, with the resulting need for abundant safety factors and non-optimal performance by all but the crudest measures.Or consider the frustration to engineers and cost to management of the never-ending need for 'code-validation' experiments every time a new class of flows is encounteredor major design change is contemplated. The whole idea of 'codes' in the first place was to be able to evaluate designs wihtout having to do experiments or build prototypes.<br />
<br />
==1.4 What do we really know for sure?==<br />
<br />
Turbulence is a subject on which still studies are going on. We really don't know a whole lot for sure about turbulence. And worse, we even disagree about what we think we know! There are indeed some things some researchers think we understand pretty well - like for example the kolmogorov similarity theory for the dissipative scales and the Law of the Wall for wall-bounded flows. These are based on assumptions and logical constructions about how we believe turbulence behaves in the limit of infinite Reynolds number. But even these ideas have never been tested in controlled laboratory ecperiments in the limits of high Reynolds number, because no one has ever had the large scale facilities required to do so.<br />
<br />
It seems to be a characteristic of humans(and contrary to popular beleif, scientists and engineers are indeed human) that we tend to accept ideas which have been around a while as fact, instead of just working hypotheses that are still waiting to be tested. One can reasonably argue that the acceptance of most ideas in turbulence is perhaps more due to the time lapsed since they were proposed and found to be in resonable agreement with limited data base, than that they have been subjected to experimental tests over the range of their assumed validity. Thus it might be wise to view most 'established' laws and theories of turbulence as more like religious creeds than matters of fact.<br />
<br />
The whole situation is a bit analogous to the old idea that the sun and stars revolved around the earth - it was a fine idea, and even good today for navigational purposes. The only problem was that one day someone(Copernicus, Brahe and Galileo among them) looked up and realized it wasn't true. So it may be with a lot of what we believe today to be true about turbulence - some day you may be the one to look at evidence in a new way and decide that things we thought to be true are wrong.<br />
<br />
== ==<br />
'''Source:Introduction to Turbulence by Professor William K.George,Chalmers Univ, Sweden.'''</div>Pavitran