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

Fluid structure interaction on baloons

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 26, 2007, 18:25
Default My OF version is build all wit
  #21
Senior Member
 
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 9
jens_klostermann is on a distinguished road
My OF version is build all with cell decomposition.

Thanks to Zjelko the couplingProperties problem is solved.

But there is another strange thing, I can only use GaussSeidel as solver for U, since this is the only solver for symmetric and assymetric matrizes.
Otherwise if I use ICCG or AMG it asks for BICCG or BDCG and vice versa. How does this come?

Jens
jens_klostermann is offline   Reply With Quote

Old   January 26, 2007, 18:43
Default It's not that bad: have anothe
  #22
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,783
Rep Power: 22
hjasak will become famous soon enough
It's not that bad: have another look.

There are two fvSolution dictionaries in the case: one in system/fvSolution and another in system/solid/fvSolution so you can control the dictionaries separately. The one in "solid" is (of course) for a solid and all should be well.

If you are not picking up the proper solver, it is possible that your solid U field (requiring a symmetric solver) is hooked onto the wrong mesh, but I doubt it.

Hope this helps,

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

Old   January 26, 2007, 19:26
Default Thanks Hrv, next question:
  #23
Senior Member
 
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 9
jens_klostermann is on a distinguished road
Thanks Hrv,

next question: Can I run the icoFsiFoam in parallel?

Jens
jens_klostermann is offline   Reply With Quote

Old   January 27, 2007, 06:29
Default I recompiled everything and no
  #24
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 9
guus_kupers is on a distinguished road
I recompiled everything and now it works. Thanks Hrv!!
guus_kupers is offline   Reply With Quote

Old   January 27, 2007, 09:26
Default Well done Guus! Looking forwar
  #25
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Well done Guus! Looking forward to your first fsi results! Maybe you can post some of your results here to share your experiences:-)

Regards,
Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   February 6, 2007, 10:55
Default Hi everybody, I'm still wor
  #26
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi everybody,

I'm still working on the development of my fluid structure interaction code. At this time, I'm begining the coupling with the open source code Impact.
Thus, I have a question which might look obvious for some of you, but I don't know how to do it.
So far, I've been using a constant vectorField defined by:
vectorField fluidPatchPointsDispl = vectorField(224, vector(0,1e-02,0));
However I would like to be able to define this vectorField by reading a file which I'd put in a correct directory, like when we define the velocity field for instance, using a dictionary. Which would be the code to define such a vectorField?

Thanks in advance,
Regards,

Vincent
vinz is offline   Reply With Quote

Old   February 6, 2007, 13:04
Default Hi Vincent, Could you pleas
  #27
New Member
 
Daniel Einstein
Join Date: Mar 2009
Posts: 22
Rep Power: 9
deinstein is on a distinguished road
Hi Vincent,

Could you please give me a link for "Impact"?

Dan
deinstein is offline   Reply With Quote

Old   February 6, 2007, 13:18
Default Hi Dan, here it is: http:/
  #28
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi Dan,

here it is:
http://impact.sourceforge.net/

you'll find number of informations ther.

Vincent
vinz is offline   Reply With Quote

Old   February 7, 2007, 12:15
Default Hi, Forget the story of vec
  #29
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi,

Forget the story of vectorField. I found a way to do what I wanted.
Now, I have another problem. I would like to find a way to know the coordinates of the points which are on a boundary and their ordering inside a field.
For example, if I have a velocity field with four points, OpenFOAM can easily write the values. But I would like to know to which point of the mesh belongs the first value, to wich point belongs the second value....
Is there a way to do that? The idea would be a function mesh.boundaryfield()[boundary1] which would return the coordinates of the points on the boundary1. But I haven't been able to find it in the code. Does such a thing allready exist?

Thanks for your help,

Vincent
vinz is offline   Reply With Quote

Old   February 26, 2007, 10:16
Default Hi everybody, This is the e
  #30
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi everybody,

This is the end of my internship and I have finished my thesis report. So I give the url if some people are interested in the work that I have done on FSI on balloons and also the converter GridPro2FOAM which is largely explained in the report.
The beginning of the report won't be very interesting for many of you since it is just OpenFOAM tutorials, but the second part of the report is more on FSI and also the beginning of coupling between OF and Impact.
Unfortunately, it's only the beginning since I spent a lot of time on the discovery and testing of OpenFOAM capacities. As you can imagine, six months of internship is really short!
Anyway, if you have any questions or comments, don't hesitate.
The report can be downloaded there: Thesis report

Best regards,

Vincent
vinz is offline   Reply With Quote

Old   February 28, 2007, 03:48
Default Hi everybody, As explained
  #31
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi everybody,

As explained sooner in this topic, I have to find the points which own to a specific boundary.
I found a way to do it in 2D doing:
mesh.boundaryMesh()[internalPatchID].boundaryPoints()

However, this is not working in 3D and, even if I have an idea, I don't really understand why it is not working.
So, does someone have an idea of how to fix it?
Thanks in advance,

Vincent
vinz is offline   Reply With Quote

Old   March 1, 2007, 03:34
Default Hi Frank, Thanks for your h
  #32
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi Frank,

Thanks for your hint. Finnaly I don't do exactly the same way but the localPoints() fonction helped me.
For information my code is now:
vectorList boundarypoint(n);

forAll (mesh1.boundaryMesh()[internalPatchID].localPoints(), j)
{
boundarypoint[j] = mesh1.boundaryMesh()[internalPatchID].localPoints()[j];
}

And this is working pefectly fine which is great since I can now deform my balloon shape in 3D.
By the way, do you plan to update your website soon? I would be really interested in the results that you obtain on moving meshes.

Best regards,

Vincent
vinz is offline   Reply With Quote

Old   April 23, 2007, 22:46
Default Hi Dear Foamers, where can
  #33
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hi Dear Foamers,

where can i get icoFsiFoam?
I searched everywhere and cannt find it. is it replaced with another solver now(version 1.4)?
And could anyone tell me what's the difference between icoFsiFoam and icoStructFoam?
Thanks in advance

Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   April 25, 2007, 06:00
Default Hi Daniel, I think the icoFsi
  #34
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12
dmoroian is on a distinguished road
Hi Daniel,
I think the icoFsiFoam solver comes only with the development version 1.3 from prof. Jasak. Try there!

Dragos
dmoroian is offline   Reply With Quote

Old   April 25, 2007, 06:48
Default Hi all, I think that there'
  #35
New Member
 
Guus Kupers
Join Date: Mar 2009
Location: Delft, Zuid-Holland, the Netherlands
Posts: 12
Rep Power: 9
guus_kupers is on a distinguished road
Hi all,

I think that there's not a big difference between the two solvers. There are some differences involving mesh motion and coupling.

And indeed icoFsiFoam is not (yet) available for OF-14.

Guus
guus_kupers is offline   Reply With Quote

Old   April 25, 2007, 07:26
Default Hi Dragos and Guus, Thanks
  #36
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hi Dragos and Guus,

Thanks for your reply, Much appreciated
When will it be available?(why not make it available) ...I will try to work on it for my thesis.
now, I have a copy of icoStructFoam and turbForceFoam, but i'm very new to Foam and i dont know is it impossible to combine LES (oodles) with such and such to solve FSI.
Any suggestion?

Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   June 2, 2008, 20:47
Default Hi, everybody, I am just ge
  #37
New Member
 
jin xu
Join Date: Mar 2009
Location: london, UK
Posts: 13
Rep Power: 9
jin_xu is on a distinguished road
Hi, everybody,

I am just gettine started with linux and OpenFoam. , also trying to run icoFsiFoam. After spending several days studying the solver and piecing togther the messages posted before, I think I have made some progress, but I still run into the same problems posted "By Guus Kupers on Monday, January 15, 2007 - 01:31 am".

--> FOAM FATAL ERROR : Attempt to cast type fixedValue to type fixedValue

From function
refCast<to>(From&)
in file /home/snail/OpenFOAM/OpenFOAM-1.3/src/OpenFOAM/lnInclude/typeInfo.H at line 103.

I realize that I should rebuild the executable of tetPolyPatchCellDecomp and tetPolyPatchFaceCellDecomp. then I type:

snail@snail-laptop:~/OpenFOAM/OpenFOAM-1.3/src/tetDecompositionFiniteElement/tetPolyMeshCellDecomp /tetPolyPatches/tetPolyPatch$ gcc tetPolyPatchCellDecomp.C

and I end up with errors:

snail@snail-laptop:~/OpenFOAM/OpenFOAM-1.3/src/tetDecompositionFiniteElement/tetPolyMeshCellDecomp /tetPolyPatches/tetPolyPatch$ gcc tetPolyPatchCellDecomp.C

In file included from tetPolyPatchCellDecomp.C:27:
tetPolyPatchCellDecomp.H:39:23: error: labelList.H: No such file or directory
tetPolyPatchCellDecomp.H:40:25: error: vectorField.H: No such file or directory
tetPolyPatchCellDecomp.H:41:25: error: triFaceList.H: No such file or directory
tetPolyPatchCellDecomp.H:76: error: expected identifier before string constant
tetPolyPatchCellDecomp.H:76: error: expected ',' or '...' before string constant
tetPolyPatchCellDecomp.H:76: error: ISO C++ forbids declaration of 'TypeName' with no type
tetPolyPatchCellDecomp.H:99: error: ISO C++ forbids declaration of 'word' with no type
tetPolyPatchCellDecomp.H:99: error: 'word' declared as a 'virtual' field
tetPolyPatchCellDecomp.H:99: error: expected ';' before '&' token
tetPolyPatchCellDecomp.H:102: error: 'label' does not name a type
tetPolyPatchCellDecomp.H:105: error: 'label' does not name a type
tetPolyPatchCellDecomp.H:108: error: 'label' does not name a type
tetPolyPatchCellDecomp.H:117: error: ISO C++ forbids declaration of 'labelList' with no type
tetPolyPatchCellDecomp.H:117: error: 'labelList' declared as a 'virtual' field
tetPolyPatchCellDecomp.H:117: error: expected ';' before '&' token
tetPolyPatchCellDecomp.H:120: error: ISO C++ forbids declaration of 'vectorField' with no type
tetPolyPatchCellDecomp.H:120: error: 'vectorField' declared as a 'virtual' field
tetPolyPatchCellDecomp.H:120: error: expected ';' before '&' token
tetPolyPatchCellDecomp.H:123: error: ISO C++ forbids declaration of 'vectorField' with no type
tetPolyPatchCellDecomp.H:123: error: 'vectorField' declared as a 'virtual' field
tetPolyPatchCellDecomp.H:123: error: expected ';' before '&' token
tetPolyPatchCellDecomp.H:127: error: 'triFaceList' does not name a type
tetPolyPatchCellDecomp.H:130: error: 'faceList' does not name a type
tetPolyPatchCellDecomp.C:36: error: expected identifier before numeric constant
tetPolyPatchCellDecomp.C:36: error: expected ',' or '...' before numeric constant
tetPolyPatchCellDecomp.C:36: error: expected constructor, destructor, or type conversion before ';' token

As I am also a beginner of linux, I am not sure my way of rebuliding the executable is correct, Could anyone give me some hints?

Thanks in advance!

jin
jin_xu is offline   Reply With Quote

Old   June 3, 2008, 10:03
Default hello again, I think I solv
  #38
New Member
 
jin xu
Join Date: Mar 2009
Location: london, UK
Posts: 13
Rep Power: 9
jin_xu is on a distinguished road
hello again,

I think I solved this problem using "./Allmake"(I should have seen it).

jin
jin_xu is offline   Reply With Quote

Old   September 14, 2016, 12:38
Default
  #39
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 35
Rep Power: 2
viraj20feb is on a distinguished road
Quote:
Originally Posted by hjasak View Post
ON, a direct question requires a direct answer. Here is a list of objects in OpenFOAM that represent a fixed value boundary:



They are all identified by the word fixedValue and you are mixing two of them.

More specifically, you are mixing:



and



Helps? (unfortunately, the forum software messes up the details of instantiation: the first one is on tetPolyPatchCellDecomp and the second on tetPolyPatchFaceDecomp.

I think I have already explained how to fix it: rebuild the executable.

Hrv

Hello,

I am trying to use icoFsiFoam for blood flow in an artery. As I give the command, the simulation stops immediately without any error message at the following step:

Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: laplace
Selecting motion diffusivity: quadratic
Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Reading field U

Reading mechanical properties

Normalising E : E/rho

Calculating Lame's coefficients

Can anyone help?
viraj20feb 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 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
Fluid-Structure Interaction Gabor BALINT Main CFD Forum 5 October 12, 1999 04:54


All times are GMT -4. The time now is 04:51.