How to access to the solid body motion state from a custom interDyMFoam.C ?

jmf · June 24, 2015, 16:07

Hi every body

I'm stuck on this coding problem since weeks

My target is to expose the solid body motion caracteristics to the solver, for instance v() or a().
I tried many ways, without success. Here is the current state of my attempts.
The idea is to expose motionPtr_ via a public function.

This is the modified dynamicMotionSolverFvMesh.H :

Code:

#ifndef dynamicMotionSolverFvMesh_H
#define dynamicMotionSolverFvMesh_H
#include "dynamicFvMesh.H"

namespace Foam
{
class motionSolver;
class dynamicMotionSolverFvMesh: public dynamicFvMesh
{
    // Private data
    autoPtr<motionSolver> motionPtr_; 

    // Private Member Functions

    //- Disallow default bitwise copy construct
    dynamicMotionSolverFvMesh(const dynamicMotionSolverFvMesh&);

    //- Disallow default bitwise assignment
    void operator=(const dynamicMotionSolverFvMesh&);

public:
    //- Runtime type information
    TypeName("dynamicMotionSolverFvMesh");

    // Constructors
        //- Construct from IOobject
        dynamicMotionSolverFvMesh(const IOobject& io);

    //- Destructor
    ~dynamicMotionSolverFvMesh();

    // Member Functions

        //- Update the mesh for both mesh motion and topology change
        virtual bool update();

        // this is the additionnal public function <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
        Foam::motionSolver* motionSolverPtr() { return motionPtr_.ptr(); };  
};
#endif

This is myInterDyMFoam.C:

Code:

Application
    interDyMFoam

#include "fvCFD.H"
#include "dynamicFvMesh.H"
#include "CMULES.H"
#include "subCycle.H"
#include "immiscibleIncompressibleTwoPhaseMixture.H"
#include "turbulenceModel.H"
#include "pimpleControl.H"
#include "fvIOoptionList.H"
#include "fixedFluxPressureFvPatchScalarField.H"
#include "dynamicMotionSolverFvMesh.H" // <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    #include "createTime.H"
    #include "createDynamicFvMesh.H"
    #include "initContinuityErrs.H"

    // this is the attempt to call motionSolverPtr() <<<<<<<<<<<<<<<<<<<<<<<
    Foam::motionSolver* motionPtr; //<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
    motionPtr = Foam::dynamicMotionSolverFvMesh::motionSolverPtr(); //<<<<<<
    Info<< "motionPtr->v() = " << motionPtr->v() << nl << endl; //<<<<<<<<<<

    pimpleControl pimple(mesh);

    #include "createFields.H"
    #include "readTimeControls.H"
    #include "createPrghCorrTypes.H"

    volScalarField rAU
    (
        IOobject
        (
            "rAU",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedScalar("rAUf", dimTime/rho.dimensions(), 1.0)
    );

    #include "correctPhi.H"
    #include "createUf.H"
    #include "CourantNo.H"
    #include "setInitialDeltaT.H"

    Info<< "\nStarting time loop\n" << endl;

    while (runTime.run())
    {
        #include "readControls.H"
        #include "alphaCourantNo.H"
        #include "CourantNo.H"

        #include "setDeltaT.H"

        runTime++;

        Info<< "Time = " << runTime.timeName() << nl << endl;

        // --- Pressure-velocity PIMPLE corrector loop
        while (pimple.loop())
        {
            if (pimple.firstIter() || moveMeshOuterCorrectors)
            {
                scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();

                mesh.update();

                if (mesh.changing())
                {
                    Info<< "Execution time for mesh.update() = "
                        << runTime.elapsedCpuTime() - timeBeforeMeshUpdate
                        << " s" << endl;

                    gh = g & mesh.C();
                    ghf = g & mesh.Cf();
                }

                if (mesh.changing() && correctPhi)
                {
                    // Calculate absolute flux from the mapped surface velocity
                    phi = mesh.Sf() & Uf;

                    #include "correctPhi.H"

                    // Make the flux relative to the mesh motion
                    fvc::makeRelative(phi, U);

                    mixture.correct();
                }

                if (mesh.changing() && checkMeshCourantNo)
                {
                    #include "meshCourantNo.H"
                }
            }

            #include "alphaControls.H"
            #include "alphaEqnSubCycle.H"

            mixture.correct();

            #include "UEqn.H"

            // --- Pressure corrector loop
            while (pimple.correct())
            {
                #include "pEqn.H"
            }

            if (pimple.turbCorr())
            {
                turbulence->correct();
            }
        }

        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

    Info<< "End\n" << endl;

    return 0;
}

Result : this does not compile. The compilation error are :

Code:

interDyMFoam.C: In function ‘int main(int, char**)’:
interDyMFoam.C:57:17: error: ‘motionSolverPtr’ is not a member of ‘Foam::dynamicMotionSolverFvMesh’
interDyMFoam.C:58:54: error: invalid use of incomplete type ‘class Foam::motionSolver’
In file included from interDyMFoam.C:44:0:
/home/nous/OpenFOAM/OpenFOAM-2.3.x/src/dynamicFvMesh/lnInclude/dynamicMotionSolverFvMesh.H:45:7: error: forward declaration of ‘class Foam::motionSolver’

I can't find a solution, I don't understand this.

And furthermore the next question is : would motionPtr->v() give something ?

I would appreciate any help

Best regards

Jean-Michel

jmf · November 19, 2016, 07:09

Hi

The solution for me was to use dictionaries in the object registry. Some example here after.

Create the dictionary and write something:

Code:

   /* interDyMFoam.C
    */

    // Creates the exchange dictionary in the registry
    dictionary dataDict;
    mesh.thisDb().store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );

   // Open and write
    const dictionary& motionDict = 
            mesh.thisDb().lookupObject<IOdictionary>("motion");

    // Update for the motion solver
    dictionary updateDbDictionary = &motionDict;
    updateDbDictionary.set("centreOfRotation", initalCentreOfMass_);
    const_cast<dictionary& > (motionDict) = updateDbDictionary ;

Read something:

Code:

/* sixDoFRigidBodyMotionSolver.C
 */
    
    if ( mesh().thisDb().foundObject<IOdictionary >("motion"))
    {
        const dictionary& dbDictionary = 
                  mesh().thisDb().lookupObject <IOdictionary >("motion");
        point centreOfRotation 
            = dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0));
    }

This works everywhere one has access to mesh or time.

The const_cast trick allows overwriting everywhere. It come from the post #14 of this thread:
http://www.cfd-online.com/Forums/ope...rfield-bc.html

Without an access to time or mesh, for instance in sixDoFRigidBodyMotion.C, I did not found better solution than adding setter and getter functions, which is hassle. It a pity that not all the classes have access to registry.

Hope this will be useful

Jean-Michel

yuhan1991 · February 23, 2017, 07:33

Hi! Thank you for your contribution,I also need to get the solid body motion state in the main Foam-solver to calculate the quaternion . But I didn't make it work when I did what you say.Could you please help me?

I added the following code into myinterDyMFoam.C:

Code:

#include "sixDoFRigidBodyMotion.H"
#include "dynamicMotionSolverFvMesh.H"

...
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Info<< "\nStarting time loop\n" << endl;

    while (runTime.run())
    {

...
// Creates the exchange dictionary in the registry
    dictionary dataDict;
    mesh.thisDb().store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );

   // Open and write
    const dictionary& motionDict = 
            mesh.thisDb().lookupObject<IOdictionary>("motion");

    // Update for the motion solver
    dictionary updateDbDictionary = &motionDict;
    updateDbDictionary.set("centreOfRotation", initalCentreOfMass_);
    const_cast<dictionary& > (motionDict) = updateDbDictionary ;

But the compiler said that:

Code:

error: ‘initialCentreOfMass_’ was not declared in this scope
     updateDbDictionary.set("centreOfRotation", initialCentreOfMass_);
                                                ^
make: *** [Make/linux64GccDPInt32Opt/myinterDyMFoam.o] Error 1

I also put the second part of codes into sixDoFRigidBodyMotionStateIO.C(in void Foam::sixDoFRigidBodyMotionState::write(Ostream& os) const ,because I thought that the motion state would be exported with the sixDoFRigidBodyMotionState )

Code:

void Foam::sixDoFRigidBodyMotionState::write(Ostream& os) const
{
    os.writeKeyword("centreOfRotation")
        << centreOfRotation_ << token::END_STATEMENT << nl;
    os.writeKeyword("orientation")
        << Q_ << token::END_STATEMENT << nl;
    os.writeKeyword("velocity")
        << v_ << token::END_STATEMENT << nl;
    os.writeKeyword("acceleration")
        << a_ << token::END_STATEMENT << nl;
    os.writeKeyword("angularMomentum")
        << pi_ << token::END_STATEMENT << nl;
    os.writeKeyword("torque")
        << tau_ << token::END_STATEMENT << nl;


    
    if ( mesh().thisDb().foundObject<IOdictionary >("motion"))
    {
        const dictionary& dbDictionary = 
                  mesh().thisDb().lookupObject <IOdictionary >("motion");
        point centreOfRotation 
            = dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0));
    }
}

But the compile didn’t pass either.The compiler gave another information:

Code:

sixDoFRigidBodyMotion/sixDoFRigidBodyMotionStateIO.C:59:15: error: ‘mesh’ was not declared in this scope
     if ( mesh().thisDb().foundObject<IOdictionary >("motion"))
               ^
sixDoFRigidBodyMotion/sixDoFRigidBodyMotionStateIO.C:59:38: error: ‘IOdictionary’ was not declared in this scope
     if ( mesh().thisDb().foundObject<IOdictionary >("motion"))

The version of OpenFOAM is 3.0.1. Could you please give me some advice? I would appreciate any help. Thank you very much.

jmf · February 25, 2017, 12:23

Hi

First you have to declare and initialize initialCentreOfMass_, possibly from dynamicMeshDict, like this :

Code:

IOdictionary dynamicMeshDict (
        IOobject (
            "dynamicMeshDict",
            runTime.constant(),
            mesh,
            IOobject::MUST_READ,
            IOobject::NO_WRITE
        )
    );
    dictionary subDict = dynamicMeshDict.subDict("tightBodyMotionCoeffs");
    vector initalCentreOfMass_ = subDict.lookup("centreOfMass");

I never tried to use sixDoFRigidBodyMotionState.

But I'm afraid that like sixDoFRigidBodyMotion, sixDoFRigidBodyMotionStateIO doesnt have any access to mesh(), mesh, neither to time or time(), which carry the registry.

If you really need to exchange data with them, a inelegant way is to write setter and getter methods in sixDoFRigidBodyMotionStateIO.C. Here is an example you will have to adapt to your needs :

Code:

void Foam::sixDoFRigidBodyMotion::set_v0(Foam::vector v) { v0_ = v;} 
Foam::vector Foam::sixDoFRigidBodyMotion::get_v0() { return v0_;}

Good luck
Jean-Michel

yuhan1991 · March 6, 2017, 09:07

Hi Fontaine! Thank you for your help again! Recently I have already made this work successfully . Now I can get the 6dof motion of the body from sixDofRigidBodyDisplacementPointPatchVectorField:: updateCoeffs() and write them to a dictionary in the main solver . So you can use them anywhere while the program is running(but the user must be able to access the object registry).This is the code snippet in the main solver:

Code:

    dictionary dataDict;
    
        runTime.store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                runTime,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );
     
    const dictionary& dbDictionary = runTime.lookupObject <IOdictionary >("motion");
    Info <<"Get centreOfRotation from database:"<< dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get initialCentreOfRotation from database:"<< dbDictionary.lookupOrDefault("initialCentreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get orientation from database:"<< dbDictionary.lookupOrDefault("orientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;
    Info <<"Get initialOrientation from database:"<< dbDictionary.lookupOrDefault("initialOrientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;

This is the modified sixDoFRigidBodyDisplacementPointPatchVectorField.C :

Code:

void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
{
........
.......
 if ( t.foundObject<IOdictionary >("motion"))
    {
       const dictionary& motionDict = 
       t.lookupObject<IOdictionary>("motion");
       dictionary updateDbDictionary = &motionDict;
       updateDbDictionary.set("centreOfRotation", motion_.state().centreOfRotation());
       updateDbDictionary.set("initialCentreOfRotation", motion_.getinitialCentreOfRotation());
       updateDbDictionary.set("orientation", motion_.orientation());
       updateDbDictionary.set("initialOrientation", motion_.getinitialQ());   
       const_cast<dictionary& > (motionDict) = updateDbDictionary ;  //without this code the result will still be ZERO!!!
       

    }

I hope this is useful.
Thank you for you help again Fontaine !

wdx_cfd · February 26, 2019, 20:45

Quote:

Originally Posted by jmf

Hi

First you have to declare and initialize initialCentreOfMass_, possibly from dynamicMeshDict, like this :

Code:

IOdictionary dynamicMeshDict (
        IOobject (
            "dynamicMeshDict",
            runTime.constant(),
            mesh,
            IOobject::MUST_READ,
            IOobject::NO_WRITE
        )
    );
    dictionary subDict = dynamicMeshDict.subDict("tightBodyMotionCoeffs");
    vector initalCentreOfMass_ = subDict.lookup("centreOfMass");

I never tried to use sixDoFRigidBodyMotionState.

But I'm afraid that like sixDoFRigidBodyMotion, sixDoFRigidBodyMotionStateIO doesnt have any access to mesh(), mesh, neither to time or time(), which carry the registry.

If you really need to exchange data with them, a inelegant way is to write setter and getter methods in sixDoFRigidBodyMotionStateIO.C. Here is an example you will have to adapt to your needs :

Code:

void Foam::sixDoFRigidBodyMotion::set_v0(Foam::vector v) { v0_ = v;} 
Foam::vector Foam::sixDoFRigidBodyMotion::get_v0() { return v0_;}

Good luck
Jean-Michel

Hi Jean-Michel,

I am facing a very similar problem recently. I have accessed the motion state by your method and I want to use it as a judgement of the stop of the PIMPLE outer loop.

For example:

I change the code:

Code:

while(pimple.loop())
{
/*  pimpleloop */
}

to

Code:

while(true)
{
/*  pimpleloop */

if(condition)         // condition is the parameter of rigid body that 
{                        // calculated in sixDOFSolver and obtained from 
    break;            // your method. In other words, the pimple outer 
}                        // loop will stop if the condition becomes true.   

}

This code has no problem when running in the serial model. However, an error occurs once the program ran parallel. The erroe is just like this:

Code:

[1] --> FOAM FATAL IO ERROR: 
[1] error in IOstream "IOstream" for operation operator>>(Istream&, List<T>&) : reading first token
[1] 
[1] file: IOstream at line 0.
[1] 
[1]     From function void Foam::IOstream::fatalCheck(const char*) const
[1]     in file db/IOstreams/IOstreams/IOstream.C at line 109.
[1] 
FOAM parallel run exiting

It seems that there is a problem of the IOstream which maybe due to the dictionary we define in the interFoam. This error appears when the first time step is finished(i.e. break by the "condition" obtained from sixDOFSolver). Furthermore, it happens at the moment of accessing the primitiveField when calculating the courant number:

Code:

fvc::surfaceSum(mag(phi))().primitiveField()

There is a problem of reading the internal field. I am using openfoam v6 and this problem makes me confused. Would you please give me some suggestions? Thank you very much!

wdx

jmf · March 2, 2019, 09:04

Hi

In parallel each processor do its own job, ignoring the existence of the other procs. You have to ensure the consistency of the data that are used for your stop condition, and after that, to publish the condition to all the processors at the same time.

So you may be interested in functions like this, to build your condition and publish it :

Code:

if (Pstream::master())
{
    x = 1;
}
Pout<< x << endl;
Pstream::scatter(x);
Pout<< x << endl;

See also in forces.C ways to gather the data from all procs:

Code:

Pstream::listCombineGather(force_, plusEqOp<vectorField>());

Hope it helps

J-Michel

wdx_cfd · March 3, 2019, 21:17

Quote:

Originally Posted by jmf

Hi

In parallel each processor do its own job, ignoring the existence of the other procs. You have to ensure the consistency of the data that are used for your stop condition, and after that, to publish the condition to all the processors at the same time.

So you may be interested in functions like this, to build your condition and publish it :

Code:

if (Pstream::master())
{
    x = 1;
}
Pout<< x << endl;
Pstream::scatter(x);
Pout<< x << endl;

See also in forces.C ways to gather the data from all procs:

Code:

Pstream::listCombineGather(force_, plusEqOp<vectorField>());

Hope it helps

J-Michel

Hi Michel,

Thank you for your advice!

I will try it recently and I will publish the result here.

Thanks a lot.

wdx.

Zhi Cheng · September 19, 2019, 20:59

Hi,yuhan! I don't know if you could get this message after long time. I am using your method, but the output was always zero. I found maybe my "motion" in my dyfoamsolver is not connected to the "motion" in sixDoFRigidBodyDisplacementPointPatchVectorField.C . Because if I change "motion" to "motion1" in my solver, it also compile. Do anyone knows why?
Best regards,
ZHI CHENG

Quote:

Originally Posted by yuhan1991

Hi Fontaine! Thank you for your help again! Recently I have already made this work successfully . Now I can get the 6dof motion of the body from sixDofRigidBodyDisplacementPointPatchVectorField:: updateCoeffs() and write them to a dictionary in the main solver . So you can use them anywhere while the program is running(but the user must be able to access the object registry).This is the code snippet in the main solver:

Code:

    dictionary dataDict;
    
        runTime.store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                runTime,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );
     
    const dictionary& dbDictionary = runTime.lookupObject <IOdictionary >("motion");
    Info <<"Get centreOfRotation from database:"<< dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get initialCentreOfRotation from database:"<< dbDictionary.lookupOrDefault("initialCentreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get orientation from database:"<< dbDictionary.lookupOrDefault("orientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;
    Info <<"Get initialOrientation from database:"<< dbDictionary.lookupOrDefault("initialOrientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;

This is the modified sixDoFRigidBodyDisplacementPointPatchVectorField.C :

Code:

void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
{
........
.......
 if ( t.foundObject<IOdictionary >("motion"))
    {
       const dictionary& motionDict = 
       t.lookupObject<IOdictionary>("motion");
       dictionary updateDbDictionary = &motionDict;
       updateDbDictionary.set("centreOfRotation", motion_.state().centreOfRotation());
       updateDbDictionary.set("initialCentreOfRotation", motion_.getinitialCentreOfRotation());
       updateDbDictionary.set("orientation", motion_.orientation());
       updateDbDictionary.set("initialOrientation", motion_.getinitialQ());   
       const_cast<dictionary& > (motionDict) = updateDbDictionary ;  //without this code the result will still be ZERO!!!
       

    }

I hope this is useful.
Thank you for you help again Fontaine !

quarkz · January 29, 2022, 19:31

Quote:

Originally Posted by yuhan1991

Hi Fontaine! Thank you for your help again! Recently I have already made this work successfully . Now I can get the 6dof motion of the body from sixDofRigidBodyDisplacementPointPatchVectorField:: updateCoeffs() and write them to a dictionary in the main solver . So you can use them anywhere while the program is running(but the user must be able to access the object registry).This is the code snippet in the main solver:

Code:

    dictionary dataDict;
    
        runTime.store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                runTime,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );
     
    const dictionary& dbDictionary = runTime.lookupObject <IOdictionary >("motion");
    Info <<"Get centreOfRotation from database:"<< dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get initialCentreOfRotation from database:"<< dbDictionary.lookupOrDefault("initialCentreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get orientation from database:"<< dbDictionary.lookupOrDefault("orientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;
    Info <<"Get initialOrientation from database:"<< dbDictionary.lookupOrDefault("initialOrientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;

This is the modified sixDoFRigidBodyDisplacementPointPatchVectorField.C :

Code:

void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
{
........
.......
 if ( t.foundObject<IOdictionary >("motion"))
    {
       const dictionary& motionDict = 
       t.lookupObject<IOdictionary>("motion");
       dictionary updateDbDictionary = &motionDict;
       updateDbDictionary.set("centreOfRotation", motion_.state().centreOfRotation());
       updateDbDictionary.set("initialCentreOfRotation", motion_.getinitialCentreOfRotation());
       updateDbDictionary.set("orientation", motion_.orientation());
       updateDbDictionary.set("initialOrientation", motion_.getinitialQ());   
       const_cast<dictionary& > (motionDict) = updateDbDictionary ;  //without this code the result will still be ZERO!!!
       

    }

I hope this is useful.
Thank you for you help again Fontaine !

Hi, I need to get the orientation (pitch) of the body in my 6dof problem at each time step. I then want this info in my modified oscillatingDisplacement mesh routine. Is it possible? If so, I need some help understanding this. I'm using overPimpleFoam.

So do I add this:

Code:

    dictionary dataDict;
    
        runTime.store
    (   
        new IOdictionary
        (   
            IOobject
            (   
                "motion",
                runTime.timeName(),
                runTime,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            dataDict
        )
    );
     
    const dictionary& dbDictionary = runTime.lookupObject <IOdictionary >("motion");
    Info <<"Get centreOfRotation from database:"<< dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get initialCentreOfRotation from database:"<< dbDictionary.lookupOrDefault("initialCentreOfRotation", point(0, 0, 0)) <<endl;
    Info <<"Get orientation from database:"<< dbDictionary.lookupOrDefault("orientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;
    Info <<"Get initialOrientation from database:"<< dbDictionary.lookupOrDefault("initialOrientation", tensor(0, 0, 0, 0, 0, 0, 0, 0, 0)) <<endl;

into the end of overPimpleFoam.c ? if not, where do i add them?

And also add:

Code:

void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
{
........
.......
 if ( t.foundObject<IOdictionary >("motion"))
    {
       const dictionary& motionDict = 
       t.lookupObject<IOdictionary>("motion");
       dictionary updateDbDictionary = &motionDict;
       updateDbDictionary.set("centreOfRotation", motion_.state().centreOfRotation());
       updateDbDictionary.set("initialCentreOfRotation", motion_.getinitialCentreOfRotation());
       updateDbDictionary.set("orientation", motion_.orientation());
       updateDbDictionary.set("initialOrientation", motion_.getinitialQ());   
       const_cast<dictionary& > (motionDict) = updateDbDictionary ;  //without this code the result will still be ZERO!!!
       

    }

into sixDoFRigidBodyDisplacementPointPatchVectorField.C ?

Lastly, how can i access the orientatio tensor from my modified oscillatingDisplacement mesh routine?

Sorry, I am a Fortran programmer so my knowledge of c++ is limited. Thanks for the help.

November 19, 2016, 07:09	[Solved] How to access to the solid body motion state from interDyMFoam.C ?	#2
jmf Member Jean-Michel FONTAINE Join Date: Aug 2009 Location: Orleans - France Posts: 55 Rep Power: 16	Hi The solution for me was to use dictionaries in the object registry. Some example here after. Create the dictionary and write something: Code: /* interDyMFoam.C / // Creates the exchange dictionary in the registry dictionary dataDict; mesh.thisDb().store ( new IOdictionary ( IOobject ( "motion", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), dataDict ) ); // Open and write const dictionary& motionDict = mesh.thisDb().lookupObject<IOdictionary>("motion"); // Update for the motion solver dictionary updateDbDictionary = &motionDict; updateDbDictionary.set("centreOfRotation", initalCentreOfMass_); const_cast<dictionary& > (motionDict) = updateDbDictionary ; Read something: Code: / sixDoFRigidBodyMotionSolver.C / if ( mesh().thisDb().foundObject<IOdictionary >("motion")) { const dictionary& dbDictionary = mesh().thisDb().lookupObject <IOdictionary >("motion"); point centreOfRotation = dbDictionary.lookupOrDefault("centreOfRotation", point(0, 0, 0)); } This works everywhere one has access to mesh or time. The const_cast trick allows overwriting everywhere. It come from the post #14 of this thread: http://www.cfd-online.com/Forums/ope...rfield-bc.html Without an access to time or mesh, for instance in sixDoFRigidBodyMotion.C, I did not found better solution than adding setter and getter functions, which is hassle. It a pity that not all the classes have access to registry. Hope this will be useful Jean-Michel zkdkeen, yuhan1991 and wdx_cfd like this. Last edited by jmf; November 19, 2016 at 07:52. Reason: Added [Solved]*

February 25, 2017, 12:23	How to access to the solid body motion state from interDyMFoam.C	#4
jmf Member Jean-Michel FONTAINE Join Date: Aug 2009 Location: Orleans - France Posts: 55 Rep Power: 16	Hi First you have to declare and initialize initialCentreOfMass_, possibly from dynamicMeshDict, like this : Code: IOdictionary dynamicMeshDict ( IOobject ( "dynamicMeshDict", runTime.constant(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE ) ); dictionary subDict = dynamicMeshDict.subDict("tightBodyMotionCoeffs"); vector initalCentreOfMass_ = subDict.lookup("centreOfMass"); I never tried to use sixDoFRigidBodyMotionState. But I'm afraid that like sixDoFRigidBodyMotion, sixDoFRigidBodyMotionStateIO doesnt have any access to mesh(), mesh, neither to time or time(), which carry the registry. If you really need to exchange data with them, a inelegant way is to write setter and getter methods in sixDoFRigidBodyMotionStateIO.C. Here is an example you will have to adapt to your needs : Code: void Foam::sixDoFRigidBodyMotion::set_v0(Foam::vector v) { v0_ = v;} Foam::vector Foam::sixDoFRigidBodyMotion::get_v0() { return v0_;} Good luck Jean-Michel Last edited by jmf; February 25, 2017 at 12:30. Reason: Was unclear

March 2, 2019, 09:04	Works in serial mode, and not in paralle	#7
jmf Member Jean-Michel FONTAINE Join Date: Aug 2009 Location: Orleans - France Posts: 55 Rep Power: 16	Hi In parallel each processor do its own job, ignoring the existence of the other procs. You have to ensure the consistency of the data that are used for your stop condition, and after that, to publish the condition to all the processors at the same time. So you may be interested in functions like this, to build your condition and publish it : Code: if (Pstream::master()) { x = 1; } Pout<< x << endl; Pstream::scatter(x); Pout<< x << endl; See also in forces.C ways to gather the data from all procs: Code: Pstream::listCombineGather(force_, plusEqOp<vectorField>()); Hope it helps J-Michel

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Calculation of the Governing Equations	Mihail	CFX	7	September 7, 2014 06:27
[DesignModeler] Question about converting surface body to solid body in DesignModeler	lnk	ANSYS Meshing & Geometry	2	November 24, 2013 00:02
[DesignModeler] how to subtract a surface body from a solid body?	yahya	ANSYS Meshing & Geometry	4	January 7, 2013 11:00
cooling a solid body	Ralf Schmidt	FLUENT	1	April 16, 2007 10:54
Two-Phase Buoyant Flow Issue	Miguel Baritto	CFX	4	August 31, 2006 12:02