[hsieh]

Hi,

To compute flow in rotating reference frame, centrifugal force and coriolis force are added to the momentum equation. If omega (rotation speed) is constant, then, one will compute these forces like the following:

dimensionedVector Omega
(
"Omega",
dimensionSet(0,0,-1,0,0),
vector(0, 2 *PI/60 *RPM, 0)
);
//centrifugal foce
volVectorField Fcent = (Omega ^ (Omega^mesh.C()));
// coriolis force
2*Omega ^ U)

What will happen if Omega varies with time? Is there any extra terms needed to be included?

Thanks!

Pei