CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

"centerOfRotation" handling in sixDoFRigidBodyMotion

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2015, 14:34
Default "centerOfRotation" handling in sixDoFRigidBodyMotion
  #1
New Member
 
Stefano
Join Date: May 2014
Location: Genova, Italy
Posts: 2
Rep Power: 0
ste_ge is on a distinguished road
Hi everyone,

I am a bit confused having an issue with the sixDoFRigidBodyMotion solver.

First of all, what is exactly, from a physical point of view, the object called "centreOfRotation"?
I presume it represents the point for which the equation that rules the rotation (i.e. Newton's second law for angular momentum) is written. In fact, all the moment arms are calculated as it was.
If I'm right, the theory says that one can choose any arbitrary point to be the center of rotation, having the same result (the moments and the moment of inertia will vary accordingly).
Of course, this point moves with the body, i.e. the distance between centreOfRotation and centreOfMass is a constant vector in the body local reference, while in the fixed reference only the magnitude of the vector is constant.

In the problem that I am studying, I need to choose "smartly" this point, so that the contribution of certain forces is zero, and the point does not coincide with the center of mass.

Now, I'm not being able to set the case successfully. I tried with setting in the constant/dynamicMeshDict an additional centreOfRotation entry, for example:
Code:
...
centreOfRotation (-5 0 0);
centreOfMass (0 0 0);
...
that I think is somehow read by the code, but I see in the output that "centreOfMass" is overwritten by "centreOfRotation".

Looking at the source code at src/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion.C (lines 103-111), I found strange the following:
Code:
...
    initialCentreOfMass_
    (
        dict.lookupOrDefault
        (
            "initialCentreOfMass",
            vector(dict.lookup("centreOfMass"))
        )
    ),
    initialCentreOfRotation_(initialCentreOfMass_),
...
Why is the centreOfRotation assumed to be coincident with the centreOfMass?

Did anyone have a similar problem before?

Any help is appreciated!
Stefano

P.S.: the solver application I use is pimpleDyMFoam.

Last edited by ste_ge; January 30, 2015 at 04:40.
ste_ge is offline   Reply With Quote

Reply

Tags
centerofrotation, pimpledyfoam, sixdofrigidbodymotion

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implicit handling of sources which couple components using Vector-Coupled Solver fumiya OpenFOAM Programming & Development 1 September 22, 2014 03:42
Tips for Handling Complex Geometries, Very Large Grids tcarrigan Pointwise & Gridgen 5 August 6, 2014 14:44
Handling bad cells at the solver level Akshay OpenFOAM Running, Solving & CFD 4 June 18, 2014 01:35
Material Handling Equipment aayushinfo56 Main CFD Forum 0 August 29, 2010 00:22
macros handling vectors sivakumar FLUENT 0 January 18, 2002 17:16


All times are GMT -4. The time now is 03:55.