[efirvida]

Hi i´m trying to convert the propeller example into a 6dof case, to see how the fluid move the propeller. I just change the dynamicMeshDict as show below, but when I run the simulation always throws --> FOAM FATAL ERROR: Could not find rho. An I don´t know where to put the rho value. because I don´t see it on any example of 6dof.
I´m using Openfoam 4.1 here is my case file modified:

dynamicMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh      dynamicMotionSolverFvMesh;

motionSolverLibs ("libsixDoFRigidBodyMotion.so");

solver            sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
    patches         ("propeller.*");
    innerDistance   0.3;
    outerDistance   1;

    mass            1.6;
    centreOfMass    (0 0 0);
    momentOfInertia (0.5 0.5 0.5);
    orientation
    (
        1 0 0
        0 1 0
        0 0 1
    );
    angularMomentum (0 0 0);
    g               (0 0 -9.81);
    rhoName         rhoInf;
    rhoInf          1;
    report          on;

    solver
    {
        type symplectic;
    }

    constraints
    {
        yAxis
        {
            sixDoFRigidBodyMotionConstraint axis;
            axis                (0 1 0);
        }
        center
        {
            sixDoFRigidBodyMotionConstraint point;
            axis                (0 0 0);
        }
    }

    restraints
    {

    }
}


// ************************************************************************* //

pointDisplacement

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       pointVectorField;
    location    "0.01";
    object      pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    "propeller.*"
    {
        type            calculated;
        value           uniform (0 0 0);
    }

    ".*"
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


// ************************************************************************* //

controlDict

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     pimpleDyMFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         5.0;

deltaT          1e-5;
writeControl    adjustableRunTime;
writeInterval   0.001;

////- For testing with moveDynamicMesh
//deltaT          0.01;
//writeControl    timeStep;
//writeInterval   1;

purgeWrite      0;

writeFormat     binary;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           2;
/*
functions
{
    #includeFunc Q
    #include "surfaces"
    #include "forces"
}
*/
// ************************************************************************* //

[Tobi]

I had to smile about your thread because two threads below yours there is the same question and I solved it yesterday already. I refer you to the solution -> http://www.cfd-online.com/Forums/ope...tml#post626241

Just as a suggestion:

• Read the error message (cannot find rho -> so FOAM already tells you what is missing)
• Use the search engine of this forum
• Check out the tutorial cases in FOAM-4.x
• If you are familiar with the source code, use Doxygen

Good luck.

[efirvida]

Originally Posted by Tobi

I had to smile about your thread because two threads below yours there is the same question and I solved it yesterday already. I refer you to the solution -> http://www.cfd-online.com/Forums/ope...tml#post626241

What happens is that yesterday did not exist that question when I wrote mine, but it was only until today that I could put it. And after seeing the answer, I did not see how to delete my question. I already have the case running to see what happens.

I'm really not familiar with the doxygen use it once but I generate too much complicated documentation. I wrote my own script to extract what I needed from the headers of the .H code, but I have not had much time to review it. Here is my generate doc but only for the $FOAM_SRC folder I have to run it on the applications folder

many thanks anyway

[bye bye my blue]

Originally Posted by efirvida

Hi i´m trying to convert the propeller example into a 6dof case, to see how the fluid move the propeller. I just change the dynamicMeshDict as show below, but when I run the simulation always throws --> FOAM FATAL ERROR: Could not find rho. An I don´t know where to put the rho value. because I don´t see it on any example of 6dof.
I´m using Openfoam 4.1 here is my case file modified:

dynamicMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh      dynamicMotionSolverFvMesh;

motionSolverLibs ("libsixDoFRigidBodyMotion.so");

solver            sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
    patches         ("propeller.*");
    innerDistance   0.3;
    outerDistance   1;

    mass            1.6;
    centreOfMass    (0 0 0);
    momentOfInertia (0.5 0.5 0.5);
    orientation
    (
        1 0 0
        0 1 0
        0 0 1
    );
    angularMomentum (0 0 0);
    g               (0 0 -9.81);
    rhoName         rhoInf;
    rhoInf          1;
    report          on;

    solver
    {
        type symplectic;
    }

    constraints
    {
        yAxis
        {
            sixDoFRigidBodyMotionConstraint axis;
            axis                (0 1 0);
        }
        center
        {
            sixDoFRigidBodyMotionConstraint point;
            axis                (0 0 0);
        }
    }

    restraints
    {

    }
}


// ************************************************************************* //

pointDisplacement

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       pointVectorField;
    location    "0.01";
    object      pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    "propeller.*"
    {
        type            calculated;
        value           uniform (0 0 0);
    }

    ".*"
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


// ************************************************************************* //

controlDict

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     pimpleDyMFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         5.0;

deltaT          1e-5;
writeControl    adjustableRunTime;
writeInterval   0.001;

////- For testing with moveDynamicMesh
//deltaT          0.01;
//writeControl    timeStep;
//writeInterval   1;

purgeWrite      0;

writeFormat     binary;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           2;
/*
functions
{
    #includeFunc Q
    #include "surfaces"
    #include "forces"
}
*/
// ************************************************************************* //

Instead of your "functions" in ContorlDict, use my script

functions

{

forces

{

type forces;

libs ( "libforces.so" );

writeControl timeStep;

writeInterval 10;

patches (walls);

rho rhoInf;

log true;

rhoInf 1;

CofR (0 0 0);

}

}

[Manideep Reddy]

Originally Posted by bye bye my blue

Instead of your "functions" in ContorlDict, use my script

functions

{

forces

{

type forces;

libs ( "libforces.so" );

writeControl timeStep;

writeInterval 10;

patches (walls);

rho rhoInf;

log true;

rhoInf 1;

CofR (0 0 0);

}

}

Thank you very much