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

Pivoting around a fixed point : 1-DOF free body

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 8, 2014, 17:39
Default Pivoting around a fixed point : 1-DOF free body
  #1
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Hello all,

I am trying to simulate flow around a simple hinged splitter in 2D.
Problem specifications are easy. please have a look at the below picture:



The shown splitter is fixed in the point shown in the above picture. It can freely rotate around it. The flow is entering the domain uniformly. And I want to simulate the fluid flow and its impact on the movement of the splitter.

I am using wingMotion tutorial as my base case. My dynamicMeshDict file looks like this:

Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / 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 (split);
innerDistance 1e-5;
outerDistance 2;

mass 1;
centreOfMass (1.02 0 0.125);
momentOfInertia (2.55 6.46 3.91);

orientation
(
1 0 0
0 1 0
0 0 1
);

velocity (0 0 0);
acceleration (0 0 0);
angularMomentum (0 0 0);
torque (0 0 0);
rhoName rhoInf;
rhoInf 1;
// g (0 -9.81 0);
report on;

constraints
{
fixedPoint
{
sixDoFRigidBodyMotionConstraint point;
point (1.02 0 0.125);
}

fixedAxis
{
sixDoFRigidBodyMotionConstraint axis;
axis ( 0 0 1 );
}
}
}
But here is the problem. I don't know why splitter would not rotate around the fixed point but around the center of mass ! This is strange, because when I suppress the fixed point constraint, splitter starts to move in the direction of the flow.

When I (wrongly) change the location of the center of mass to the point I want the splitter to rotate around, everything is fine.

Is there anything I am missing ?!

I would be really grateful if you can help me on this.
Thank you,
Mojtaba.
Attached Images
File Type: jpg Splitter_CFD.jpg (32.4 KB, 247 views)
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   April 10, 2014, 16:41
Default
  #2
Member
 
Join Date: Apr 2009
Posts: 36
Rep Power: 8
minger is on a distinguished road
How do you know that things are not working? It seems that this flow is stable -- I presume the steady-state solution would have the splitter in the initial position?
minger is offline   Reply With Quote

Old   April 11, 2014, 08:30
Default
  #3
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Quote:
Originally Posted by minger View Post
How do you know that things are not working? It seems that this flow is stable -- I presume the steady-state solution would have the splitter in the initial position?
Well yes, In steady state solution, splitter is in the initial position. But the in the second run when I use pimpleDyMFoam, It starts rotating around the center of mass as it can be seen through post processing.

As I mentioned I have intentionally set the point I want the splitter to rotate around, to be the center of mass, and it starts rotating just fine; but I have some serious problems in convergence. Angular velocity blows up suddenly.

I don't know if changing the center of mass to the point of rotation causes this bad behavior, thus I am just willing to solve this problem first.
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   April 11, 2014, 18:32
Default
  #4
Member
 
Join Date: Apr 2009
Posts: 36
Rep Power: 8
minger is on a distinguished road
I'm not sure I fully understand the problem. You can run the simulation steady just fine. However, you run unsteady with the dynamic mesh and it blows up.

In the dictionary you posted above, the center of mass is the same as the center of rotation. This I'm guessing blows up.

What happens when the center of rotation is "correct"? Does it blow up? Does it just not do what you expect?

For your case that's working but blowing up, try enabling maxCo -- because that particular situation is very unstable, as soon as it starts to rotate it should accelerate quickly.
minger is offline   Reply With Quote

Old   May 2, 2014, 13:00
Default
  #5
New Member
 
Nikola Pilipovic
Join Date: Nov 2013
Posts: 6
Rep Power: 3
Nik_984 is on a distinguished road
Hi Mojtaba,

Have you maybe solved the problem, I am dealing with the similar case where my object is not rotating around point that I have defined.
BTW I just NOTICED that your fixedPoint is SAME as the centre of mass(in the example you are showing in your post), maybe that's your problem.
Nik_984 is offline   Reply With Quote

Old   May 2, 2014, 13:41
Default
  #6
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Quote:
Originally Posted by minger View Post
I'm not sure I fully understand the problem. You can run the simulation steady just fine. However, you run unsteady with the dynamic mesh and it blows up.

In the dictionary you posted above, the center of mass is the same as the center of rotation. This I'm guessing blows up.

What happens when the center of rotation is "correct"? Does it blow up? Does it just not do what you expect?

For your case that's working but blowing up, try enabling maxCo -- because that particular situation is very unstable, as soon as it starts to rotate it should accelerate quickly.
Quote:
Originally Posted by Nik_984 View Post
Hi Mojtaba,

Have you maybe solved the problem, I am dealing with the similar case where my object is not rotating around point that I have defined.
BTW I just NOTICED that your fixedPoint is SAME as the centre of mass(in the example you are showing in your post), maybe that's your problem.
Dear Nikola and minger,
Unfortunately I haven't been able to solve the problem yet.
Yes I am aware of that. well even by changing the location of the fixed point or center of mass, solid just rotates around the center of mass and not the fixed point.
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   May 3, 2014, 04:05
Default
  #7
New Member
 
Nikola Pilipovic
Join Date: Nov 2013
Posts: 6
Rep Power: 3
Nik_984 is on a distinguished road
I have tried with the same aproach(center of mass coordinates = point of rotation), and I got rotation around that point, I will take a look at source code and try to find what's going on with it, if I have any progress, I will let you know.
Nik_984 is offline   Reply With Quote

Old   June 25, 2014, 09:11
Default
  #8
Member
 
Join Date: Jun 2012
Posts: 48
Rep Power: 5
maHein is on a distinguished road
It looks like this issue has been addressed in the commit 3207f8d0c7d83396cebf1958fb38de732c8b0307 :-)
maHein is offline   Reply With Quote

Old   July 10, 2014, 10:54
Default
  #9
New Member
 
Stefano
Join Date: May 2014
Location: Genova, Italy
Posts: 2
Rep Power: 0
ste_ge is on a distinguished road
Hi!

I had the same problem. A few days ago I noticed (in the latest version of OF 2.3.x, so you could have to update the sources) the existence of the keyword "centreOfRotation".

You can try to use something like:

Code:
constraints
{

   fixedPoint
   {
       sixDoFRigidBodyMotionConstraint point;
       centreOfRotation (0 0 0);
   }

}
I still have not played deeply with this feature but I think it can solve the problem.

Cheers
ste_ge is offline   Reply With Quote

Reply

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
Body Fitted Grid (BFG) vs Free Gridding (FG) CK Main CFD Forum 0 September 5, 2008 02:57
demo free flow blunt body in cfx ansys 11 jan CFX 1 July 31, 2007 19:44
gravitational force for free surface flow Jongtae Kim Main CFD Forum 1 July 2, 2000 11:57
Variable Density - Free Surface with FIDAP Vitaliy Pavlyk FLUENT 7 May 2, 2000 15:56
Any Body got a Free Surface Solution in FIDAP Dr. Ahmed Hassaneen Main CFD Forum 2 November 8, 1999 15:59


All times are GMT -4. The time now is 02:39.