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

Fluid Structure interaction

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 1, 2007, 08:53
Default Hello all, For my graduatio
  #1
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 8
guus_kupers is on a distinguished road
Hello all,

For my graduation thesis I 'm going to use a fluid structure interaction solver. So far I'm using icoStructFoam. First I want to validate this solver. To do this I have a 2D test case where water flows through a small channel in which a rigid cylinder is placed which in turn is connected to a deformable "flag". Therefore I would like to know whether some one else has run this or a similar case before as well.

Because the results I obtained differ quite a lot compared to a benchmark test executed by others. Thus did anyone else run this problem as well? And did you encounter some problems?

I'm not sure but can it be a problem that the fluid which is used is water which has a 1000 times higher density? And that because of this the communication between the fluid and the structure is not correct?

Thanks in advance!!

Guus
guus_kupers is offline   Reply With Quote

Old   March 6, 2007, 16:49
Default Hi Guus! Can it be possible
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Guus!

Can it be possible, that you ran into the "viscosity and pressure have to be divide by rho when using an incompressible solver"? Everone's been there at least once

But if it's not that there still might be a problem: If I remember it correctly (and we're refering to the same version of icoStructFoam - the one on the Wiki) the fluid pressure just gets transfered to the solid - without taking into account the density of the fluid.

Replacing the line

solidP[exchange[fI]]=fluidP[fI];

with something like

solidP[exchange[fI]]=fluidRho*fluidP[fI];

might give the correct results (the fluidRho you can hardcode or read from a dictionary).

Sorry for the mixup. But as it says on the Wiki: that thing should only be considered a demo.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   March 14, 2007, 09:46
Default I've tried this, but now my so
  #3
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 8
guus_kupers is on a distinguished road
I've tried this, but now my solution blows up almost immediately after the movement of the mesh starts.

Do you have any idea why this happens??
guus_kupers is offline   Reply With Quote

Old   March 19, 2007, 08:41
Default It looks like my problem is st
  #4
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 8
guus_kupers is on a distinguished road
It looks like my problem is stable now. However I would like to run either icoFsiFoam or icoStructFoam in parallel, has anyone tried this before?

Because I run into problems when the 2 seperate regions need to be coupled. Has anyone done this before or does anyone know how a multiple region problem can be run in parallel on 2 or more nodes??

This is the error message I get:

[0]
[0]
[0] --> FOAM FATAL ERROR : Sizes do not match
[0]
[0] From function
[0] Coupled Solver
[0] in file readCoupling.H at line 47.
[0]
FOAM parallel run exiting
[0]
[1]
[1]
[1] --> FOAM FATAL ERROR : Sizes do not match
[1]
[1] From function
[1] Coupled Solver
[1] in file readCoupling.H at line 47.
[1]
FOAM parallel run exiting
[1]

Thanks in advance!!
guus_kupers is offline   Reply With Quote

Old   March 19, 2007, 14:34
Default Another disclaimer I should pu
  #5
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Another disclaimer I should put on the icoStructFoam-Wiki-page (in fact, I just did). It is not parallelized. It is possible, I think, to parallelize it, but at first you should think about decomposing the case: if both regions are decomposed in such a way that fluid-zone on processor 42 has EXACTLY the same fluid-solid interface as the solid-zone on processor 42, it should be relatively(!) easy to parallelize the solver. On the other hand, I don't think that decomposing the case in such a way is generally applicable.

BTW: How did you decompose your case for the parallel run? I was not aware that decompesPar supports multiple mesh regions.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   March 20, 2007, 04:24
Default I don't know either whether de
  #6
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 8
guus_kupers is on a distinguished road
I don't know either whether decomposePar supports multiple mesh regions but I did it by hand. Thus I decomposed the seperate regions with a "dummy" case and after that I combined them again into one parallel case.

My case does not have a complicated mesh, so it would be possible to put the structure mesh completely inside one of the fluid meshes. In this way my meshes have the exact same interface. However it means that I don't have any mesh for my other processors. Is it possible to leave them empty in some way??

Guus
guus_kupers is offline   Reply With Quote

Old   March 20, 2007, 13:15
Default Hi Guus! Yes it would be po
  #7
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Guus!

Yes it would be possible (I think, I haven't tried, and I won't for reasons given below): You can decompose your fluid mesh. Put all stuff in processorX/constant/polyMesh into a processorX/region1/constant/polyMesh (or is it the other way round?). If you want to calculate the solid phase on only one processor put the mesh you have from constant/region2/polymesh into processor42/region2/constant/polyMesh and for all processors x!=42 create an "empty" (cell lists etc have length 0, but must be present) mesh in processorX/region2/constant/polymesh. Still I would be pleasantly surprised if that works.

Still the problem remains (as discussed elsewhere on this board) that if the solid phase is calculated on only one processor, you're bound to get very inefficient parallelization.

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   March 20, 2007, 13:46
Default Hi Berhard, In fact, the de
  #8
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Hi Berhard,

In fact, the decomposition is only a part of the problem. For proper you still need parallelism, you still need to make sure that the bits of interpolation get properly reduced to the right processors + you need to take into account the cases where the processor boundary goes right across the sluid-structure interface.

There's an easy way around the problem, but I'm toofussy to do it that way. What you do is to make a single mesh out of both fluid and structure domain and using fvMeshSubset make a separate fluid and structure domain out of it. (I did something similar for complex cases where you've got lots of bits of over-lapping and non-overlapping physics: a fuel cell). You then decompose the underlying mesh into parallel domains and use the global mesh for data mapping - this will parallelise naturally.

On the negative side, I don't feel like doing this right now. Maybe we can talk about it on the FSI session on the Workshop in Zagreb (if there's sufficient interest)...

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 21, 2007, 08:39
Default Hi Hrv! My solution was ref
  #9
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Hrv!

My solution was refering to the very special topology Guus was refering to: only one processor has the complete fluid-struct interface. And for that case icoStructFoam MIGHT work (although more by accident than by design)

But that fvMeshSubset looks interesting. So there is always a collateral benefit with that kind of discussions
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   March 21, 2007, 09:30
Default I've tried to do what you expl
  #10
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 8
guus_kupers is on a distinguished road
I've tried to do what you explained. I emptied out all the existing mesh files and also I created an empty cells file for region2. However when icoStructFoam tries to create the "empty" mesh for region2, the solid region, I get the following error message:

--> FOAM FATAL ERROR : points deallocated

From function
const pointField& polyMesh::allPoints() const
in file meshes/polyMesh/polyMesh.C at line 642.
FOAM parallel run aborting

I think this has to do with the fact that the mesh is empty. Is there a way to avoid this or have I done something wrong??

Guus
guus_kupers is offline   Reply With Quote

Old   March 21, 2007, 09:32
Default Beware: the original version o
  #11
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Beware: the original version of fvMeshSubset I wrote a while ago really annoyed me because of the class hierarchy it implies (I could not run off the subset)... so I've changed it. As an additional benefit, you get all the mapping done for free, but there are issues with multiple boundaries I am still thinking about.

Give me a shout if you want to take it for a spin.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   June 17, 2011, 11:53
Default my solver
  #12
New Member
 
Rasoul
Join Date: Nov 2010
Posts: 3
Rep Power: 6
rasoulgharaei is on a distinguished road
Hi everyone,
I made new solver (interStructFoam) by combining the interFoam and solidDisplacementFaom According to the icoStructFoam Instruction, but I could not compile my solver in OP-1.7.0.
Can anyone help me?
I Attach my solver with this message.


regards,
rasoul gharaei
Attached Files
File Type: zip my solver.zip (48.1 KB, 11 views)
rasoulgharaei 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
fluid structure interaction taru agrawal FLUENT 4 September 10, 2007 03:12
FLUID STRUCTURE INTERACTION Javier Main CFD Forum 5 October 11, 2004 04:22
Fluid & Structure Interaction... nesusa FLUENT 1 March 31, 2002 21:30
Fluid-Structure Interaction for VOF Zi-Wei Chiou Main CFD Forum 4 November 27, 2001 10:26
Fluid Structure Interaction Delly Main CFD Forum 2 November 26, 2001 13:01


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