Difference between revisions of "See the MRF development"
From OpenFOAMWiki
Line 26: | Line 26: | ||
Navier-Stokes equations in the inertial frame | Navier-Stokes equations in the inertial frame | ||
+ | |||
+ | <math> | ||
+ | \begin{cases} | ||
+ | \frac {D \vec u_I}{D t} = - \nabla (p/\rho) + \nu \nabla \cdot \nabla (\vec u_I) \\ | ||
+ | \nabla \cdot \vec u_I = 0 | ||
+ | \end{cases} | ||
+ | </math> '''Eqn [2]''' | ||
<math> | <math> | ||
Line 46: | Line 53: | ||
\nabla \cdot \vec u_I = 0 | \nabla \cdot \vec u_I = 0 | ||
\end{cases} | \end{cases} | ||
− | </math> | + | </math> '''Eqn [3]''' |
− | + | ||
+ | Let's look at the left-hand side of the momentum equation of Eqn [2], by taking into account Eqn [1] for the acceleration term: | ||
+ | |||
+ | <math>\frac {D \vec u_I}{D t} = \frac{D \vec u_R}{Dt} + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r</math> | ||
+ | |||
+ | <math>\frac {D \vec u_I}{D t} = \frac {\partial \vec u_R}{\partial t} + \vec u_R \cdot \nabla \vec u_R + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r</math> | ||
+ | |||
+ | <math>\frac {D \vec u_I}{D t} = \frac {\partial \vec u_R}{\partial t} + \nabla \cdot (\vec u_R \otimes \vec u_R) + \frac{d \vec \Omega}{dt} \times \vec r + 2 \vec \Omega \times \vec u_R + \vec \Omega \times \vec \Omega \times \vec r</math> since <math>\nabla \cdot \vec u_R = \nabla \cdot \vec u_I = 0</math> | ||
+ | |||
+ | <math>\nabla \cdot \vec u_I = \nabla \cdot \left [ \vec u_R + \vec \Omega \times \vec r \right ]</math> | ||
+ | |||
+ | <math>\nabla \cdot \vec u_I = \nabla \vec u_R + \underbrace{\nabla \left [ \vec \Omega \times \vec r \right ]}_{0}</math> | ||
+ | |||
+ | <math>\nabla \cdot \vec u_I = \nabla \vec u_R</math> |
Revision as of 20:34, 25 May 2009
Reynolds-Averaged Navier-Stokes formulation in the rotating frame.
Acceleration term expressed for a rotating frame around the z axis ().
Notation: I: inertial, R: rotating
General vector:
Position vector:
For the accelration, the velocity vector is:
Eqn [1]
Navier-Stokes equations in the inertial frame
Eqn [2]
Eqn [3]
Let's look at the left-hand side of the momentum equation of Eqn [2], by taking into account Eqn [1] for the acceleration term:
since