FSI on a rotating reference frame
 

I'm trying to study the start up performance of turbines. I want the instantaneous force balance (of the turbine) to assign an angular velocity to the rotating reference frame which encompasses the turbine patch. Turbine is rigid and does not deform.

Here's what I'm thinking to do..

1. Use pimpleDymFoam with the dynamicMotionSolverFvMesh
-this should produce the acceleration of the turbine/RRF
-use acceleration as in input to the solidBodyMotionFvMesh with a rotatingMotion as a motionFunction.
-aka update the radialVelocity in rotatingMotionCoeffs each time step.


I want to perform these actions in the sourcecode and develop a new solver. Any suggestions as to the procedure to "merge these two motion solvers (which are coupled)" in a new solver?
- which headers file to include from each?
- how to pass the acceleration/RPM data between solvers in the int main()
- I have good C++ OOP experience but not used to the built in types, classes and structures.

Any suggestions will be greatly appreciated... for now I'll just look at the "guts" of dynamicMotionSolverFvMesh & solidBodyMotionFvMesh and how their called.