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

Sliding Mesh or rotor stator model implimentation

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By lovecraft22

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 11, 2012, 08:58
Default Sliding Mesh or rotor stator model implimentation
  #1
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
Hi Fomers,
This a very serious question guys, and would appreciate suggestions from every one.

What is state of implementing sliding mesh in OpenFOAM ( or with OpenFOAM201)?

We require to implement VAWT ( Vertical Axis Wind Turbine), Fixing the outer domain and rotating a cylindrical Mesh inside it having blade.stl .
Now another question comes up:
What should be the approach ? ( in getting 2 meshes ie cylinder with stl blades and outer domain)



Regards
CFDkid

PS: Any other suggestions !!!!
kid is offline   Reply With Quote

Old   April 11, 2012, 10:41
Default
  #2
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
I would suggest you to have a look at the following tutorials:
Multi Reference Frames: tutorials/incompressible/MRFSimpleFoam/mixerVessel2D
Moving mesh: tutorials/incompressible/pimpleDyMFoam/propeller
lovecraft22 is offline   Reply With Quote

Old   April 12, 2012, 03:30
Default
  #3
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore

Was using openFOAM201 , @/tutorial/incompressible/pimpleDyFoam did not find propeller case. Checked out by upgrading to openFoam210 and found the propeller case, working on it now. Will come up with question soon.
Also i decided to go for sliding mesh instead of moving Frame of Reference i.e the two case you suggested, i am trying out propeller one with pimpleDyFoam.
Is my decision ok ??

regards,
CFDkid


kid is offline   Reply With Quote

Old   April 12, 2012, 05:02
Default
  #4
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Yep, the propeller tutorial came out with the latest version of OF.
About your decision, it really depend on whether you want to run a steady or unsteady simulation.
lovecraft22 is offline   Reply With Quote

Old   April 12, 2012, 07:19
Default
  #5
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore
I am just doing that propeller case, but things would get prity complicated for me if i implement a similar case.
Few questions!!!
1. What is the role of importing geometry , i do not know if at all it is required? ( At least the object which has to rotate ( suppose a turbine)).

2. How to i create the innerCylinder and innerCylinderSmall? (the sliding faces that would separate Rotor and Stator part of CFD domain)

3.How to create mesh
"surfaces=
innerCylinder
innerCylinderSmall
outerCylinder
propellerTip
propellerStem1
propellerStem2
propellerStem3
"
?

I use slome for geometry and mesh creation + simpleFoam GUI in CAELinux to for the OpenFoam case directory.
else I use discritizer.

But how should i approach for this case, really stuck man.
regards,

Last edited by kid; April 12, 2012 at 07:45. Reason: Clearity was needed
kid is offline   Reply With Quote

Old   April 12, 2012, 07:31
Default
  #6
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Have you looked already at the allRun script? That's all you need to run…
lovecraft22 is offline   Reply With Quote

Old   April 12, 2012, 07:41
Default
  #7
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore
Yes i did dig inside. Also ran the case.
Everything fine with that case.

But if i want to implement another similar case.
Like a rotating cube in a regular inlet ---> outlet domain.

Now if you see that question of mine above, you would understand what i intended to ask.

regards
CFdkid
kid is offline   Reply With Quote

Old   April 12, 2012, 08:16
Default
  #8
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
The easiest way is probably to enclose the mesh you want to rotate within a .stl cylinder then, in snappyHexMeshDict you'll need something like this:

Code:
refinementSurfaces
    {
        cylinder
        {
            level       (2 2);
            cellZone    cylinder;
            faceZone    cylinder;
            cellZoneInside  inside;
        }


    }
and then in your dynamicMeshDict:

Code:
solidBodyMotionFvMeshCoeffs
{
    cellZone        cylinder;

    solidBodyMotionFunction  rotatingMotion;
    rotatingMotionCoeffs
    {
        CofG        (0 0 0);
        radialVelocity (0 0 34); // deg/s
    }
}
HenningW likes this.
lovecraft22 is offline   Reply With Quote

Old   April 12, 2012, 08:26
Default
  #9
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore
That means , if we have a small cube to be rotated, we need to create a mesh file and enclose that mesh in cylindrical stl file. Then import that stl with mesh in OpenFoam and carry out the step you said.
kid is offline   Reply With Quote

Old   April 12, 2012, 09:51
Default
  #10
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Yes. I would you snappyHexMesh.
lovecraft22 is offline   Reply With Quote

Old   April 16, 2012, 08:01
Default
  #11
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@Lore
Man, used snappyHexMesh to generate this cube with cylindrical mesh around mvCube.stl ( should have named it rotatingCube!! forgot).

mvCube.stl was made using Salome and stl binary was placed in /constant/ triSurface/
and accordingly changes were made in snappyHexMeshDict::geometry

Now is it sufficient and should i move to edit dynamicMeshDict? Or something missing?

Regards
CFDkid

@Lore
Thanks a ton man , this helpped me to learn snaapyHexMesh , though rotating a cube is still pending.
Attached Images
File Type: png snap_mvCube.png (5.2 KB, 145 views)
kid is offline   Reply With Quote

Old   April 16, 2012, 09:11
Default
  #12
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Are you sure you enclosed your cube within a cylindrical surface? Cannot see it in the image you attached…
lovecraft22 is offline   Reply With Quote

Old   April 16, 2012, 09:31
Default
  #13
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
Quote:
Originally Posted by lovecraft22 View Post
Are you sure you enclosed your cube within a cylindrical surface? Cannot see it in the image you attached…
SnappyHexMeshDict::
************************************************** *************
geometry
{
mvCube.stl
{
type triSurfaceMesh;
name mvCube;
}

cylinder1 //This cylinder is used.
{
type searchableCylinder;
point1 (100 50 0);
point2 (100 50 10);
radius 30;
}
}
castellatedMeshControls::
refinementSurfaces
{ mvCube
{
level (2 4);
}

}
refinementRegions
{

cylinder1
{
mode inside;
levels ((1e+15 2));
}
}

***************************************
This is what was used.
When you say covering the cube with cylinder. I am interpreting the one implimented.
But if you differ from this implementation. How to go about it ? (enclosed your cube within a cylindrical surface)

Regards
kid is offline   Reply With Quote

Old   April 16, 2012, 09:40
Default
  #14
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@Lore
Were you talking about this cylinder in the image attached. Now i could see it.
The patch you mentioned in your previous reply, stuck me now after your reply.
Is it correct now?

regards
Attached Images
File Type: png snap_mvCube-2.png (9.7 KB, 70 views)
kid is offline   Reply With Quote

Old   April 16, 2012, 09:48
Default
  #15
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Yes, should be ok now. I'm running a similar case, I'll enclose it here in case you wanted to have a look. Unfortunately I'm struggling with it because it keeps blowing up… anyway the setup should be ok.
My cube is called cubo.stl and my cylinder is called cylinder.stl.

Refer to the allRun file for the commands I ran.
Attached Files
File Type: zip 2_cubo.zip (35.6 KB, 64 views)
lovecraft22 is offline   Reply With Quote

Old   April 17, 2012, 04:30
Default
  #16
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore
Correct man, it is blowing up.

Help man. But do not know why?

Regards

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.0-0bc225064152
Exec : pimpleDyMFoam
Date : Apr 17 2012
Time : 12:34:13
Host :
PID : 26986
Case pen-foam/mvCube
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone cylinder1
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::solidBodyMotionFvMesh::solidBodyMotionFvMesh (Foam::IOobject const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#4 Foam::dynamicFvMesh::addIOobjectConstructorToTable <Foam::solidBodyMotionFvMesh>::New(Foam::IOobjec t const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#5 Foam::dynamicFvMesh::New(Foam::IOobject const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#6
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#7 __libc_start_main in "/lib/libc.so.6"
#8
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
Segmentation fault
kid is offline   Reply With Quote

Old   April 17, 2012, 05:41
Default
  #17
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Don't know… is your setup the same as mine?
lovecraft22 is offline   Reply With Quote

Old   April 17, 2012, 06:16
Default
  #18
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@lore
My setup was not same as your will try again.
Why pimpleDyFoam is so challenging to implement
Also looking at the thesis mentioned in this thread. (http://www.cfd-online.com/Forums/ope...ng-motion.html)

@lore
Should we also explore AMI feature , which i suppose should be implimented, as GGi is no longer present in
OpenFoam-210.
Do not know if atall these two things are related to this case, we are solving.
In propeller tutorial saw some commands in Allrun.pre which we have not implemented.
Could be because of that solver is crashing.

regards

Last edited by kid; April 17, 2012 at 07:56. Reason: missing thoughts
kid is offline   Reply With Quote

Old   April 17, 2012, 18:28
Default
  #19
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 460
Rep Power: 18
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Hi, seems like my setup was wrong.

Have a look here:
http://www.cfd-online.com/Forums/ope...lowing-up.html

and especially here:
http://www.cfd-online.com/Forums/ope...tml#post355232

Sorry about that…
lovecraft22 is offline   Reply With Quote

Old   April 18, 2012, 01:41
Default
  #20
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
@Lore
Seeing that discussion in a flash. This is cool man, need to put lot of effort
from my side.

Thanks a lot and great guiding line , now it feels this is right track.
Would work now and let you know the status soon.

regards
kid is offline   Reply With Quote

Reply


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 07:42
Sliding mesh error alice FLUENT 3 February 21, 2012 06:53
[snappyHexMesh] external flow with snappyHexMesh chelvistero OpenFOAM Meshing & Mesh Conversion 11 January 15, 2010 20:43
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09


All times are GMT -4. The time now is 01:15.