Favre averaged NavierStokes equations
From CFDWiki
Line 95:  Line 95:  
</table>  </table>  
  Note that the corresponding expression  +  Note that the corresponding expression (15) for Favre averaged turbulent flows contains an extra term related to the turbulent energy. 
Equations (1)(9), supplemented with gas data for <math>\gamma</math>, <math>Pr</math>, <math>\mu</math> and perhaps <math>R</math>, form a closed set of partial differential equations, and need only be complemented with boundary conditions.  Equations (1)(9), supplemented with gas data for <math>\gamma</math>, <math>Pr</math>, <math>\mu</math> and perhaps <math>R</math>, form a closed set of partial differential equations, and need only be complemented with boundary conditions.  
Line 130:  Line 130:  
Note that with the above definitions <math>\overline{\Phi'} = 0</math>, but <math>\overline{\Phi''} \neq 0</math>.  Note that with the above definitions <math>\overline{\Phi'} = 0</math>, but <math>\overline{\Phi''} \neq 0</math>.  
+  
+  == Open Turbulent Equations ==  
+  
+  In order to obtain an averaged form of the governing equations, the instantaneous continuity equation (1), momentum  
+  equation (2) and energy equation (3) are timeaveraged. Introducing a density weighted time average decomposition (10) of <math>u_i</math> and <math>e_0</math>, and a standard time average decomposition (11) of <math>\rho</math> and <math>p</math> gives the following exact open equations:  
+  
+  <table width="100%"><tr><td>  
+  :<math>  
+  \frac{\partial \overline{\rho}}{\partial t} +  
+  \frac{\partial}{\partial x_i}\left[ \overline{\rho} \widetilde{u_i} \right] = 0  
+  </math>  
+  </td><td width="5%">(12)</td></tr></table>  
+  
+  <table width="100%"><tr><td>  
+  :<math>  
+  \frac{\partial}{\partial t}\left( \overline{\rho} \widetilde{u_i} \right) +  
+  \frac{\partial}{\partial x_j}  
+  \left[  
+  \overline{\rho} \widetilde{u_i} \widetilde{u_j} + \overline{p} \delta_{ij} +  
+  \overline{\rho u''_i u''_j}  \overline{\tau_{ji}}  
+  \right]  
+  = 0  
+  </math>  
+  </td><td width="5%">(13)</td></tr></table>  
+  
+  <table width="100%"><tr><td>  
+  :<math>  
+  \frac{\partial}{\partial t}\left( \overline{\rho} \widetilde{e_0} \right) +  
+  \frac{\partial}{\partial x_j}  
+  \left[  
+  \overline{\rho} \widetilde{u_j} \widetilde{e_0} +  
+  \widetilde{u_j} \overline{p} + \overline{u''_j p} +  
+  \overline{\rho u''_j e''_0} + \overline{q_j}  \overline{u_i \tau_{ij}}  
+  \right]  
+  = 0  
+  </math>  
+  </td><td width="5%">(14)</td></tr></table>  
+  
+  The density averaged total energy <math>\widetilde{e_0}</math> is given by:  
+  
+  <table width="100%"><tr><td>  
+  :<math>\widetilde{e_0} \equiv \widetilde{e} + \frac{\widetilde{u_k} \widetilde{u_k}}{2} + k</math>  
+  </td><td width="5%">(15)</td></tr></table>  
+  
+  Where the turbulent energy, $k$, is defined by~(\ref{eq:turb_energy}):  
+  
+  <table width="100%"><tr><td>  
+  :<math>k \equiv \fav{\frac{u''_k u''_k}{2}}</math>  
+  </td><td width="5%">(16)</td></tr></table>  
+  
+  Equation (12), (13) and (14) are referred to as the Favre averaged NavierStokes equations. <math>\overline{\rho}</math>, <math>\widetilde{u_i}</math> and <math>\widetilde{e_0}</math> are the primary solution variables.  
+  
+  == Approximations ==  
<math>  <math> 
Revision as of 09:22, 5 September 2005
Contents 
Instantaneuos Equations
The instantaneous continuity equation (1), momentum equation (2) and energy equation (3) for a compressible fluid can be written as:
 (1) 
 (2) 
 (3) 
For a Newtonian fluid, assuming Stokes Law for monoatomic gases, the viscous stress is given by:
 (4) 
Where the traceless viscous strainrate is defined by:
 (5) 
The heatflux, , is given by Fourier's law:
 (6) 
Where the laminar Prandtl number is defined by:
 (7) 
To close these equations it is also necessary to specify an equation of state. Assuming a calorically perfect gas the following relations are valid:
 (8) 
Where , , and are constant.
The total energy is defined by:
 (9) 
Note that the corresponding expression (15) for Favre averaged turbulent flows contains an extra term related to the turbulent energy.
Equations (1)(9), supplemented with gas data for , , and perhaps , form a closed set of partial differential equations, and need only be complemented with boundary conditions.
Favre Averaged Equations
It is not possible to solve the instantaneous equations directly for the applications of interest here. At the Reynolds numbers typically present in a turbine these equations have very chaotic turbulent solutions, and it is necessary to model the influence of the smallest scales. All turbulence models used in this work are based on onepoint averaging of the instantaneous equations. The averaging procedure used is described in the next sections.
Averaging
Let be any dependent variable. It is convenient to define two different types of averaging of :
 Classical time average (Reynolds average):

(10) 

 Density weighted time average (Favre average):
 (11) 

Note that with the above definitions , but .
Open Turbulent Equations
In order to obtain an averaged form of the governing equations, the instantaneous continuity equation (1), momentum equation (2) and energy equation (3) are timeaveraged. Introducing a density weighted time average decomposition (10) of and , and a standard time average decomposition (11) of and gives the following exact open equations:
 (12) 
 (13) 
 (14) 
The density averaged total energy is given by:
 (15) 
Where the turbulent energy, $k$, is defined by~(\ref{eq:turb_energy}):
 (16) 
Equation (12), (13) and (14) are referred to as the Favre averaged NavierStokes equations. , and are the primary solution variables.
Approximations