CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

OpenFoam MRF Torque calculations and assumptions

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   September 24, 2019, 12:03
Default OpenFoam MRF Torque calculations and assumptions
Join Date: Apr 2011
Posts: 55
Rep Power: 14
Ruonin is on a distinguished road
Hi folks,

I know this topic has been thoroughly discussed in several posts but I have not seem a definite answer but rather conflicting information.

I would like to calculate the Power Number for a multiple level impeller.
I use simpleFoam with two MRF zones one per impeller level.
This is straightforward as seen in several threads.

I use the following (typical) code in controlDict in duplicated

functions (
forces 1
{ type forces;
functionObjectLibs ("");
outputControl timeStep;
outputInterval 1;
patches (impeller bottom);
rhoName rhoInf;
rhoInf 1000;
CofR (-0.35 -0.55 2.07);
} );

1. My first assumption is the following: my understanding is that CofR are the coordinates of the center of the impeller (i.e., the intersection of the impeller with the shaft)?

After convergence we get the following data (dummy data for example purposes):

forces output:
forces(pressure,viscous)((-0.1 -0.2 -0.3),(-0.1 -0.2 -0.3),(-0.1 -0.2 -0.3))
moment(pressure,viscous)((-0.1 -0.2 -0.3),(-0.1 -0.2 -0.3),(-0.1 -0.2 -0.3))

My second assumption is this: to get the torque (assuming Z axis as center of rotation) we sum the red values in the above result, correct?

Assuming my assumptions are correct, I simulated my case with 1M+ nodes and got a convergence without issues.

However, when I calculate the power number for both levels individually the bottom impeller value differs from the value supplied by the vendors

The vendor states ~0.5 Power Number for the bottom impeller and I'm getting around 0.27.

I double checked everything but no luck on understanding the difference.

My BC for the simulation are:
- symmetry for the fluid surface
- wall for reactor walls and baffles
- movingWall with -6 rad/s for the shaft and impeller levels (individually)
- MRF with - 6 rad/s

Can anyone provide some insight?


Ruonin is offline   Reply With Quote

Old   September 26, 2019, 14:19
Join Date: Apr 2011
Posts: 55
Rep Power: 14
Ruonin is on a distinguished road
Dare to bump!
Ruonin is offline   Reply With Quote


mrf simplefoam, mrf. moving wall, openfoam 4.1, power number

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

All times are GMT -4. The time now is 21:05.