CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

oscillating cylinder dynamic mesh

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

Like Tree23Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 29, 2014, 13:07
Default
  #21
Member
 
Petr Furmanek
Join Date: Jan 2012
Location: Prague, Czech Republic
Posts: 63
Rep Power: 6
petr.f. is on a distinguished road
I'd say there is FSI in the standard version - take a look at $FOAM_TUTORIALS//incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam

... or do I get something wrong?
petr.f. is offline   Reply With Quote

Old   May 1, 2014, 09:15
Default
  #22
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,736
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Greetings Petr,

Quote:
Originally Posted by petr.f. View Post
I'd say there is FSI in the standard version
FSI has a more broad approach, where it allows for solid shapes to deform depending on the interaction with fluids.
That particular tutorial only relies on a rigid-solid approach (it's named "sixDoFRigidBodyMotion" in "constant/dynamicMeshDict"), which means that even though the fluid may move the solid, said solid will not change its own shape.

The difficulty with FSI itself is usually with the solid deformation using an FVM (Finite Volume Method) approach, instead of the conventional FEM (Finite Element Method), as a reflex to the interaction of forces made by the fluid on the solid.

Best regards,
Bruno
iy-a likes this.
wyldckat is offline   Reply With Quote

Old   May 1, 2014, 09:33
Default
  #23
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Dear Bruno,

many thanks for your clarification.

I have question about my case. I already finished forced oscillation of circular cylinder (case attached). Now, I have been stuck in applying the natural vibration of my case. I tried to use wingMotion as a reference, but I didn't get the result. could you help me how to apply natural (free) vibration. I knew that I have to do some changes in pointDisplacement file referring to wingMotion test case. But unfortunately, I'm still stucking on that because I lake of information of different fields in pointDisplacement file in wingMotion folder.

Many thanks in advanced. Looking forward to get your answer.

Maimouna
Attached Files
File Type: gz oscillatingCylinder.tar.gz (3.7 KB, 127 views)
Maimouna is offline   Reply With Quote

Old   May 2, 2014, 12:08
Default
  #24
Member
 
Petr Furmanek
Join Date: Jan 2012
Location: Prague, Czech Republic
Posts: 63
Rep Power: 6
petr.f. is on a distinguished road
Hi Bruno - that's what I've meant, just a simple case of FSI. Anyway - I haven't check the extend version, are they using FVM or FEM for structural deformations?
petr.f. is offline   Reply With Quote

Old   May 3, 2014, 09:46
Default
  #25
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,736
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Greetings to all!

@Maimouna:
Quote:
Originally Posted by Maimouna View Post
could you help me how to apply natural (free) vibration. I knew that I have to do some changes in pointDisplacement file referring to wingMotion test case. But unfortunately, I'm still stucking on that because I lake of information of different fields in pointDisplacement file in wingMotion folder.
You'll have to be more specific on what exactly you want to do, because I can think of several kinds of natural vibrations for such a situation:
  1. A natural vibration occurring in the fluid (ripple effect?).
  2. The fluid inducing a natural vibration on the solid cylinder (force driven?).
  3. Moving the cylinder at a previously known natural vibration frequency.
  4. Structure deformation/destruction caused by reaching a natural vibration.

@Petr:
Quote:
Originally Posted by petr.f. View Post
Hi Bruno - that's what I've meant, just a simple case of FSI. Anyway - I haven't check the extend version, are they using FVM or FEM for structural deformations?
FSI stands for Fluid-Structure Interaction. It's briefly explained here: http://www.foamcfd.org/development.html

The "simple case of FSI" you're referring to would be Fluid-Solid Interaction... which is not the specific reason why the acronym "FSI" was created . More specifically, what you're referring to is usually addressed as dynamic meshing, as the changes that are made to the geometrical boundaries will enforce changes on the mesh, or vice versa.

According to the description on that page I gave the link above, the structure side is performed using FEM in foam-extend, which currently is no yet operational in the official OpenFOAM version.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   May 7, 2014, 11:32
Default
  #26
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Quote:
You'll have to be more specific on what exactly you want to do, because I can think of several kinds of natural vibrations for such a situation:
  1. A natural vibration occurring in the fluid (ripple effect?).
  2. The fluid inducing a natural vibration on the solid cylinder (force driven?).
  3. Moving the cylinder at a previously known natural vibration frequency.
  4. Structure deformation/destruction caused by reaching a natural vibration.
Dear Bruno,

I'm looking at case number 2, but this is not force. Any idea please?

Many thanks in advanced.
Maimouna is offline   Reply With Quote

Old   May 22, 2014, 12:30
Default
  #27
Member
 
Giovanni Caramia
Join Date: Mar 2009
Location: Bari, ITALY
Posts: 37
Rep Power: 9
karamiag is on a distinguished road
Hi ok,
I downloaded your oscillating cylinder case, but I am not able to make the cylinder move. Could you please give more details about how to launch the case?
Thank you.
karamiag is offline   Reply With Quote

Old   May 22, 2014, 12:35
Default
  #28
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Quote:
Hi ok,
I downloaded your oscillating cylinder case, but I am not able to make the cylinder move. Could you please give more details about how to launch the case?
Thank you.
Hi,

the cylinder oscillates properly in my OpenFOAM versions (2.2.2 and 2.3.0).

Which OF version do you use?
Maimouna is offline   Reply With Quote

Old   May 22, 2014, 12:37
Default
  #29
Member
 
Giovanni Caramia
Join Date: Mar 2009
Location: Bari, ITALY
Posts: 37
Rep Power: 9
karamiag is on a distinguished road
I am using OpenFOAM version 2.3.0. In particular I have not found any information in your case about oscillation frequency or amplitude.
karamiag is offline   Reply With Quote

Old   May 22, 2014, 12:42
Default
  #30
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Strange
I just tried now to make sure. It's working. Did you use pimpleDyMFoam solver?
Maimouna is offline   Reply With Quote

Old   May 22, 2014, 12:48
Default
  #31
Member
 
Giovanni Caramia
Join Date: Mar 2009
Location: Bari, ITALY
Posts: 37
Rep Power: 9
karamiag is on a distinguished road
This is what I have done.
tar xf oscillatingCylinder.tar.gz
cd oscillatingCylinder
blockMesh
pimpleDyMFoam
And this is what I obtain.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0
Exec : pimpleDyMFoam
Date : May 22 2014
Time : 19:44:06
Host : "ljjjjjkllkjl"
PID : 21970
Case : /home/tesisti/OpenFOAM/tesisti-2.3.0/run/oscillatingCylinder
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 staticFvMesh


--> FOAM FATAL IO ERROR:
keyword PIMPLE is undefined in dictionary "/home/tesisti/OpenFOAM/tesisti-2.3.0/run/oscillatingCylinder/system/fvSolution"

file: /home/tesisti/OpenFOAM/tesisti-2.3.0/run/oscillatingCylinder/system/fvSolution from line 23 to line 74.

From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 643.

FOAM exiting
karamiag is offline   Reply With Quote

Old   May 22, 2014, 12:55
Default
  #32
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
This is the steps:
1. Extract tar xf oscillatingCylinder.tar.gz
2. Copy and paste in pimpleDyMFoam solver in incomressible tutorial
3. In terminal window alias 'of230'
4. Go to the case path.
5. ./Allclean
6. ./Allprepare
7. pimpleDyMFoam
8. paraFoam

Hope this help.

Regards
Maimouna is offline   Reply With Quote

Old   May 22, 2014, 13:01
Default
  #33
Member
 
Giovanni Caramia
Join Date: Mar 2009
Location: Bari, ITALY
Posts: 37
Rep Power: 9
karamiag is on a distinguished road
Extracting oscillatingCylinder.tar.gz I obtain a directory in which there is no Allclean or Allprepare file. Is it possible the targz file you uploaded is the wrong one?
karamiag is offline   Reply With Quote

Old   May 22, 2014, 15:00
Default
  #34
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Try this one. I hope it will help.

All the best
Attached Files
File Type: gz oscillatingCylinder.tar.gz (5.2 KB, 138 views)
karamiag, Bazinga and alirezahd like this.
Maimouna is offline   Reply With Quote

Old   May 23, 2014, 02:36
Default
  #35
Member
 
Giovanni Caramia
Join Date: Mar 2009
Location: Bari, ITALY
Posts: 37
Rep Power: 9
karamiag is on a distinguished road
OK, now it works!
Thank you.
karamiag is offline   Reply With Quote

Old   June 15, 2014, 12:05
Default
  #36
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,736
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Greetings to all!

@Maimouna:
Quote:
Originally Posted by Maimouna View Post
I'm looking at case number 2, but this is not force. Any idea please?
Sorry, time has been flying past me and only now did I manage to give a quick look into this. I see that you've continued to look into this and that you've found this thread: Rigid Body Mesh Motion + 6DoF

Essentially that thread does point in the right direction. My suggestion is that you also look into the tutorial case "multiphase/interDyMFoam/ras/floatingObject", comparing the versions on OpenFOAM 2.1, 2.2 and 2.3, because it has evolved over time. The idea in that case is that the floating body is moving due to impulsion of the fluid vs body density.

There is another example case... let's see if I can find it... ah, it's this thread: Questions on dynamicTopoFvMesh - that example was for a sphere that was moving inside a tube towards a hole. The detail in this thread is that it uses a specific dynamic motion library, namely this one: https://github.com/smenon/dynamicTopoFvMesh

Best regards,
Bruno
Haier and Scabbard like this.
__________________
wyldckat is offline   Reply With Quote

Old   June 16, 2014, 06:14
Default
  #37
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Dear Bruno,

many thanks for your post. But, I'm still confusing how to implement that in the cylinder case. Have I create a new solver or my case? What is the exactly steps to get the result that I'm looking.

I'm really need some help. I lack of these technical steps.

Many thanks in advanced.

Maimouna
Maimouna is offline   Reply With Quote

Old   June 16, 2014, 06:27
Default
  #38
Member
 
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 4
Scabbard is on a distinguished road
Quote:
Originally Posted by Maimouna View Post
Dear Bruno,

many thanks for your post. But, I'm still confusing how to implement that in the cylinder case. Have I create a new solver or my case? What is the exactly steps to get the result that I'm looking.

I'm really need some help. I lack of these technical steps.

Many thanks in advanced.

Maimouna
Hi Maimouna,

I have already make the case running use the sixDoFRigidBodyMotion solver, for me I looked at the WingMotions in the tutorials. For now, I just make it successfully run but have not make the validation of the result. If you need any help, may be I can help you with this case.

Scabbard
Scabbard is offline   Reply With Quote

Old   June 16, 2014, 06:31
Default
  #39
Member
 
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 4
Scabbard is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Greetings to all!

@Maimouna:

Sorry, time has been flying past me and only now did I manage to give a quick look into this. I see that you've continued to look into this and that you've found this thread: Rigid Body Mesh Motion + 6DoF

Essentially that thread does point in the right direction. My suggestion is that you also look into the tutorial case "multiphase/interDyMFoam/ras/floatingObject", comparing the versions on OpenFOAM 2.1, 2.2 and 2.3, because it has evolved over time. The idea in that case is that the floating body is moving due to impulsion of the fluid vs body density.

There is another example case... let's see if I can find it... ah, it's this thread: Questions on dynamicTopoFvMesh - that example was for a sphere that was moving inside a tube towards a hole. The detail in this thread is that it uses a specific dynamic motion library, namely this one: https://github.com/smenon/dynamicTopoFvMesh

Best regards,
Bruno
Dear Bruno,

Thanks for your keeping help. Now I want to record the trajectory of the cylinder's motion response, the centre of the cylinder displacement. Is there any function in the OpenFoam can do this? Thank you.

Best wishes,
Scabbard
Scabbard is offline   Reply With Quote

Old   June 16, 2014, 07:37
Default
  #40
Senior Member
 
ok
Join Date: Oct 2013
Posts: 344
Rep Power: 6
Maimouna is on a distinguished road
Quote:
Hi Maimouna,

I have already make the case running use the sixDoFRigidBodyMotion solver, for me I looked at the WingMotions in the tutorials. For now, I just make it successfully run but have not make the validation of the result. If you need any help, may be I can help you with this case.

Scabbard
Hi Scabbard,

of course, I need your help. I'm struggling on that for long time. If you don't mind please we can shear your idea and I would be really grateful.

My email address is: may78may@hotmail.com

Best regards
Maimouna
Maimouna 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
Dynamic Mesh "Shadow Wall" thezack FLUENT 0 June 4, 2013 22:09
dynamic mesh for drop interface IndrajitW FLUENT 0 March 30, 2013 09:03
dynamic mesh refinement and rhoCentralFoam ChrisA OpenFOAM Running, Solving & CFD 1 March 21, 2013 09:00
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10


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