- **OpenFOAM Running, Solving & CFD**
(*http://www.cfd-online.com/Forums/openfoam-solving/*)

- - **Theoretical background of MRF library**
(*http://www.cfd-online.com/Forums/openfoam-solving/58045-theoretical-background-mrf-library.html*)

Hello,
I've derived the momHello,
I've derived the momentum equation for the relative frame of reference formulated with the absolute velocity here. The result http://www.cfd-online.com/OpenFOAM_D...your_image.gif is coincident with the literature but not with the formulation in the fluent manual. The second term on the left hand side is different, whereas the rest is equal. It can be shown that the tensor 2. rank from the dyadic product of (vec{c} \otimes vec{w}) is not symmetric, if you write it down component-by-component in cartesian coordinates, you will see that (vec{c} \otimes vec{w} \neq vec{w} \otimes vec{c}). So there seems to be a bug in the fluent manual. That's why I think, the formulation of the UEqn in the MRFSimpleFoam solver is also not valid. There is fvm::div(phi, U), in this case phi denotes to w_i (relative velocity) for each spatial direction. Am I right? |

There are two different derivations of Navier-Stokes equations for so called "Multiple Reference Frames" (MRF) or "Multiple Frames of Reference" (MFR):
- The first notation is equation 6 form here, which is also implemented in the MRF library:
- The second notation is equation 21 form here:
If anybody wants to use the actual implementation of the MRF library within a transient solver (that includes all compressible transsonic solvers), one has to be certain, that the local time derivation is calculated with the relative velocity () and not with the absolute velocity () ! As this could be a little bit tricky, I would prefer either the second notation of the MRF or using the formulation with the relative velocity in the relative frame of reference, as this is done in the SRF library, with an appropriate rotor-stator interface. |

All times are GMT -4. The time now is 17:30. |