|
|
| Line 1: |
Line 1: |
| {{Refimprove|date=November 2010}}
| | Hi there! :) My name is Bell, I'm a student studying Nutritional Sciences from Schoonhoven, Netherlands.<br><br>my site: [http://www.party24-7.com/board_Qyic63/675426 Www.Bookbyte.Com Coupon] |
| | |
| The '''Reynolds-averaged Navier–Stokes equations ''' (or '''RANS''' equations) are time-averaged<ref>
| |
| The true time average (<math> \bar{X} </math>) of a variable (<math> x </math>) is defined by
| |
| : <math> \bar{X} = \lim_{T \to \infty}\frac{1}{T}\int_{t_0}^{t_0+T} x\, dt.</math>
| |
| For this to be a well-defined term, the limit (<math> \bar{X} </math>) must be independent of the initial condition at <math>t_0</math>. In the case of a [[chaotic dynamical system]], which the equations under turbulent conditions are thought to be, this means that the system can have only one strange attractor, a result that has yet to be proved for the Navier-Stokes equations. However, assuming the limit exists (which it does for any bounded system, which fluid velocities certainly are), there exists some <math>T</math> such that integration from <math>t_0</math> to <math>T</math> is arbitrarily close to the average. This means that given transient data over a sufficiently large time, the average can be numerically computed within some small error. However, there is no analytical way to obtain an upper bound on <math>T</math>.</ref>
| |
| equations of motion for [[fluid flow]]. The idea behind the equations is [[Reynolds decomposition]], whereby an instantaneous quantity is decomposed into its time-averaged and fluctuating quantities, an idea first proposed by [[Osborne Reynolds]].<ref>
| |
| Reynolds, Osborne, 1895: "On the Dynamical Theory of Incompressible Viscous Fluids and the Determination of the Criterion." ''Philosophical Transactions of the Royal Society of London. A'', v. 186, pp. 123-164. Available online from [http://www.jstor.org/stable/90643 JSTOR].</ref> The RANS equations are primarily used to describe [[turbulent flow]]s. These equations can be used with approximations based on knowledge of the properties of flow [[turbulence]] to give approximate time-averaged solutions to the [[Navier–Stokes equations]].
| |
| For a [[Stationary process|stationary]], incompressible [[Newtonian fluid]], these equations can be written in [[Einstein notation]] as:
| |
| | |
| : <math>\rho\bar{u}_j \frac{\partial \bar{u}_i }{\partial x_j}
| |
| = \rho \bar{f}_i
| |
| + \frac{\partial}{\partial x_j}
| |
| \left[ - \bar{p}\delta_{ij}
| |
| + \mu \left( \frac{\partial \bar{u}_i}{\partial x_j} + \frac{\partial \bar{u}_j}{\partial x_i} \right)
| |
| - \rho \overline{u_i^\prime u_j^\prime} \right ].
| |
| </math>
| |
| | |
| The left hand side of this equation represents the change in mean momentum of fluid element owing to the unsteadiness in the [[mean flow]] and the convection by the mean flow. This change is balanced by the mean body force, the isotropic stress owing to the mean pressure field, the viscous stresses, and apparent stress <math> \left( - \rho \overline{u_i^\prime u_j^\prime} \right)</math> owing to the fluctuating velocity field, generally referred to as the [[Reynolds stresses|Reynolds stress]]. This nonlinear Reynolds stress term requires additional modeling to close the RANS equation for solving, and has led to the creation of many different [[turbulence models]]. The time-average operator <math>\overline{.}</math> is a [[Reynolds operator]].
| |
| | |
| == Derivation of RANS equations ==
| |
| The basic tool required for the derivation of the RANS equations from the instantaneous [[Navier–Stokes]] equations is the [[Reynolds decomposition]]. Reynolds decomposition refers to separation of the flow variable (like velocity <math>u</math>) into the mean (time-averaged) component (<math>\overline{u}</math>) and the fluctuating component (<math>u^{\prime}</math>). Because the mean operator is a [[Reynolds operator]], it has a set of properties. One of these properties is that the mean of the fluctuating quantity being equal to zero(<math>\bar{u^\prime} = 0</math>). Thus,
| |
| | |
| :<math> u(\boldsymbol{x},t) = \bar{u}(\boldsymbol{x}) + u^\prime(\boldsymbol{x},t) \,</math>, where <math> \boldsymbol{x} = (x,y,z) </math> is the position vector. Some authors<ref>{{cite book|last=Tennekes|first=H.|title=A first course in turbulence|year=1992|publisher=MIT Press|location=Cambridge, Mass. [u.a.]|isbn=978-0-262-20019-6|edition=14. print.|coauthors=Lumley, J. L.}}</ref> prefer using <math>U</math> instead of <math> \bar{u} </math> for the mean term (since an overbar is sometimes used to represent a vector). In this case, the fluctuating term <math>u^\prime</math> is represented instead by <math>u</math>. This is possible because the two terms do not appear simultaneously in the same equation. To avoid confusion, the notation <math> u, \bar{u}, \mbox{ and } u^{\prime} </math> will be used to represent the instantaneous, mean, and fluctuating terms, respectively.
| |
| | |
| The properties of [[Reynolds operator]]s are useful in the derivation of the RANS equations. Using these properties, the Navier–Stokes equations of motion, expressed in tensor notation, are (for an incompressible Newtonian fluid):
| |
| | |
| : <math> \frac{\partial u_i}{\partial x_i} = 0 </math>
| |
| : <math> \frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j}
| |
| = f_i
| |
| - \frac{1}{\rho} \frac{\partial p}{\partial x_i}
| |
| + \nu \frac{\partial^2 u_i}{\partial x_j \partial x_j}
| |
| </math>
| |
| | |
| where <math>f_i</math> is a vector representing external forces.
| |
| | |
| Next, each instantaneous quantity can be split into time-averaged and fluctuating components, and the resulting equation time-averaged,
| |
| | |
| <ref>
| |
| Splitting each instantaneous quantity into its averaged and fluctuating components yields,
| |
| | |
| : <math> \frac{\partial \left( \bar{u_i} + u_i^\prime \right)}{\partial x_i} = 0 </math>
| |
| | |
| : <math> \frac{\partial \left( \bar{u_i} + u_i^\prime\right)}{\partial t}
| |
| + \left( \bar{u_j} + u_j^\prime\right) \frac{\partial \left( \bar{u_i} + u_i^\prime\right)}{\partial x_j}
| |
| = \left( \bar{f_i} + f_i^\prime\right)
| |
| - \frac{1}{\rho} \frac{\partial \left(\bar{p} + p^\prime\right)}{\partial x_i}
| |
| + \nu \frac{\partial^2 \left( \bar{u_i} + u_i^\prime\right)}{\partial x_j \partial x_j}.
| |
| </math>
| |
| | |
| Time-averaging these equations yields,
| |
| : <math> \overline{\frac{\partial \left( \bar{u_i} + u_i^\prime \right)}{\partial x_i}} = 0 </math>
| |
| : <math> \overline{\frac{\partial \left( \bar{u_i} + u_i^\prime\right)}{\partial t}}
| |
| + \overline{\left( \bar{u_j} + u_j^\prime\right) \frac{\partial \left( \bar{u_i} + u_i^\prime\right)}{\partial x_j}}
| |
| = \overline{\left( \bar{f_i} + f_i^\prime\right)}
| |
| - \frac{1}{\rho} \overline{\frac{\partial \left(\bar{p} + p^\prime\right)}{\partial x_i}}
| |
| + \nu \overline{\frac{\partial^2 \left( \bar{u_i} + u_i^\prime\right)}{\partial x_j \partial x_j}}.
| |
| </math>
| |
| | |
| Note that the nonlinear terms (like <math> \overline{u_i u_i} </math>) can be simplified to,
| |
| | |
| <math> \overline{u_i u_i}
| |
| = \overline{\left( \bar{u_i} + u_i^\prime \right)\left( \bar{u_i} + u_i^\prime \right) }
| |
| = \overline{\bar{u_i}\bar{u_i} + \bar{u_i}u_i^\prime + u_i^\prime\bar{u_i} + u_i^\prime u_i^\prime}
| |
| = \bar{u_i}\bar{u_i} + \overline{u_i^\prime u_i^\prime} </math>
| |
| </ref> to yield:
| |
| | |
| : <math> \frac{\partial \bar{u_i}}{\partial x_i} = 0</math>
| |
| : <math> \frac{\partial \bar{u_i}}{\partial t}
| |
| + \bar{u_j}\frac{\partial \bar{u_i} }{\partial x_j}
| |
| + \overline{u_j^\prime \frac{\partial u_i^\prime }{\partial x_j}}
| |
| = \bar{f_i}
| |
| - \frac{1}{\rho}\frac{\partial \bar{p}}{\partial x_i}
| |
| + \nu \frac{\partial^2 \bar{u_i}}{\partial x_j \partial x_j}. </math>
| |
| | |
| The momentum equation can also be written as,<ref>
| |
| This follows from the mass conservation equation which gives,
| |
| : <math> \frac{\partial u_i}{\partial x_i} = \frac{\partial \bar{u_i}}{\partial x_i} + \frac{\partial u_i^\prime}{\partial x_i} = 0</math>
| |
| </ref>
| |
| | |
| : <math> \frac{\partial \bar{u_i}}{\partial t}
| |
| + \bar{u_j}\frac{\partial \bar{u_i} }{\partial x_j}
| |
| = \bar{f_i}
| |
| - \frac{1}{\rho}\frac{\partial \bar{p}}{\partial x_i}
| |
| + \nu \frac{\partial^2 \bar{u_i}}{\partial x_j \partial x_j}
| |
| - \frac{\partial \overline{u_i^\prime u_j^\prime }}{\partial x_j}.
| |
| </math>
| |
| On further manipulations this yields,
| |
| : <math>\rho \frac{\partial \bar{u_i}}{\partial t}
| |
| + \rho \bar{u_j} \frac{\partial \bar{u_i} }{\partial x_j}
| |
| = \rho \bar{f_i}
| |
| + \frac{\partial}{\partial x_j}
| |
| \left[ - \bar{p}\delta_{ij}
| |
| + 2\mu \bar{S_{ij}}
| |
| - \rho \overline{u_i^\prime u_j^\prime} \right ]
| |
| </math>
| |
| | |
| where,
| |
| <math>
| |
| \bar{S_{ij}} = \frac{1}{2}\left( \frac{\partial \bar{u_i}}{\partial x_j} + \frac{\partial \bar{u_j}}{\partial x_i} \right)
| |
| </math>
| |
| is the mean rate of strain tensor.
| |
| | |
| Finally, since integration in time removes the time dependence of the resultant terms, the time derivative must be eliminated, leaving:
| |
| : <math>\rho \bar{u_j}\frac{\partial \bar{u_i} }{\partial x_j}
| |
| = \rho \bar{f_i}
| |
| + \frac{\partial}{\partial x_j}
| |
| \left[ - \bar{p}\delta_{ij}
| |
| + 2\mu \bar{S_{ij}}
| |
| - \rho \overline{u_i^\prime u_j^\prime} \right ].
| |
| </math>
| |
| | |
| ==Notes==
| |
| {{Reflist}}
| |
| | |
| {{DEFAULTSORT:Reynolds-averaged Navier-Stokes equations}}
| |
| [[Category:Fluid dynamics]]
| |
| [[Category:Turbulence]]
| |
| [[Category:Turbulence models]]
| |