[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

[hsieh]

Hi,

One more question about this.

If the computational domain locates R away from the axis of rotation, where R >> l (l is the characteristic length of the computational domain), how does it affect the centrifugal force/coriolis force calculation?

Pei

[hsieh]

OK,

I found the answer.

Another force term due to angular acceleration needs to be added as follows:

d(Omega)/dt ^ mesh.C()

pei