[dntjudme]

Dear Foamers,

I'm using OF v7 multiSolidBodyMotionSolver for simulating a wind turbine aerodynamics case.

What I want to do is to simulate the wind turbine rotating motion + surge motion simultaneously, I referred to what OpenFOAM wiki suggested in section 2.8 multiMotion: https://openfoamwiki.net/index.php/P...sh#multiMotion

Following what wiki suggested, I setup the dynamicMeshDict as follows, which has been modified slightly in OF v7 style

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolver    //solidBody;
                multiSolidBodyMotionSolver;

solidBodyMotionFunction  multiMotion;
                surgeBox
                {
                    solidBodyMotionFunction  oscillatingLinearMotion;
                    oscillatingLinearMotionCoeffs
                    {
                      amplitude     (0.125 0 0);//(0.125 0 0);
                      omega         0.785398; // (rad/s) 0.125 Hz = 0.785398 rad/s
                    }
                }

                rotor
                {
                    solidBodyMotionFunction  rotatingMotion;
                    rotatingMotionCoeffs
                    {
                      origin          (0 0 2.086);
                      axis            (1 0 0);
                      omega           25.1327;//240; //25.1327 rad/s
                    }
                    solidBodyMotionFunction  oscillatingLinearMotion;
                    oscillatingLinearMotionCoeffs
                    {
                      amplitude     (0.125 0 0);//(0.125 0 0);
                      omega         0.785398; // (rad/s) 0.125 Hz = 0.785398 rad/s
                    }
                }

The assignment of cellZone "rotor" is done in snappyHexMeshDict, includes cells inside the AMI patch; the "surgeBox" is assigned in topoSet with the following code, the line "box (-9 -7 0) (11 7 4);" covers all cells in the entire computational domain:

{
        name    surgeBoxSource;
        type    cellSet;
        action  new;
        source  boxToCell;
        sourceInfo
        {
            box (-9 -7 0) (11 7 4);
        }
    }
    {
        name    surgeBox;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        sourceInfo
        {
            set surgeBoxSource;
        }
    }

Then, I ran the moveDynamicMesh command to preview the mesh motion in ParaView, and here is what I got：1png.jpg

Here you can see that the rotation of the cellZone "rotor" is not happening while the surge motion can be properly simulated.

Does anyone know how to correctly tune the multiMotion Coeffs for this case so that the mesh motion can be properly moved with rotating motion + surge motion occurring simultaneously?

Many thanks!

Zs