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

icoFsiFoam on OpenFOAM 1.6

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 1, 2010, 15:34
Default icoFsiFoam on OpenFOAM 1.6
  #1
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi all,

Does anyone can compile icoFsiFoam on OpenFOAM 1.6 ?
Please tell me how to compile it or modifications related to the successful compilation.

I compiled successfully icoStructFoam on the 1.6 and seems to get fine results of two test cases attached the icoStructFoam, but hope to use icoFsiFoam.

Thanks in advance.

waku2005

Last edited by waku2005; June 1, 2010 at 19:03.
waku2005 is offline   Reply With Quote

Old   June 18, 2010, 11:18
Default
  #2
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Hi Waku, I've also compiled icoStructFoam on OpenFOAM-1.6, but I have some problems with the matter of running the blockMesh utility for both of the regions of the thingOnStick test case. Could you please give me some advices about it?

Thanks in advance
vkrastev is offline   Reply With Quote

Old   June 19, 2010, 01:46
Default
  #3
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi, vkratev,

I've used the icoStructFoam downloaded from
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/PhilipEvegren/IcoStructFoam_Rev561.tgz, and it works fine.

Instructionsof the icoSructFoam described in the openfoam wiki will works good
in the thingOnStick test case which was attached to the above tar-ball.

Here's my steps to run IcoStructFoam:
> tar zxvf IcoStructFoam_Rev561
> cd coStructFoam_Rev561
> cd icoStructFoam (wmake in this directory)
>
> cd ../
> pwd
> /home/ofuser/OpenFOAM/ofuser-1.6/my_orig/IcoStructFoam_Rev561
> ln -s region1/polyMesh/ thingOnStick/constant/polyMesh
> blockMesh -case thingOnStick
> rm thingOnStick/constant/polyMesh
>
> ln -s region2/polyMesh/ thingOnStick/constant/polyMesh
> blockMesh -case thingOnStick
> rm thingOnStick/constant/polyMesh
>
> cd thingOnStick
> icoStructFoam (It will run successfully without any editing dictionaries)

I've read in the E.Mogura's Blog (Japanese site) that icoStructFoam from the openfoam wiki could not be compiled on OF 1.6.

Last edited by waku2005; June 19, 2010 at 05:07.
waku2005 is offline   Reply With Quote

Old   June 24, 2010, 03:36
Default
  #4
New Member
 
Thomas Dalla Piazza
Join Date: Jun 2010
Posts: 1
Rep Power: 0
tdalla is on a distinguished road
Hi, I'm also interested to compile the icoFsiFoam for the version 1.6. I will try to compile the one from the version 1.5....
tdalla is offline   Reply With Quote

Old   June 28, 2010, 07:07
Default icoFsiFoam
  #5
New Member
 
Karl Jacob Maus
Join Date: Jun 2010
Location: Norway
Posts: 2
Rep Power: 0
karlmaus is on a distinguished road
Hi, all.

As far as I understand it (please correct me if I'm wrong, anyone ), the icoFsiFoam solver from the extend-project depends on several parts not present in the official version of OpenFOAM (OF), most notably the tetrahedral finite element routines.

The reason for this is that for a proper fluid-structure interaction (FSI) solver you would want the interface boundary to be calculated as correct as possible, hence the nodal (vertex) -based solution method of tetFEM for the mesh and thus the boundary.
To use only the cell-centered finite volume method (FVM) in OF could result in incorrect boundary calculations and wrong mesh motion (e.g., non-conserved volumes, wrong motion of boundary due to interpolated solid velocities, etc.).

Another approach would be to use a vertex-based FVM, but AFAIK this is not present in OF (yet).
karlmaus is offline   Reply With Quote

Old   October 11, 2010, 09:53
Default
  #6
Member
 
Marta's Avatar
 
Marta Lazzarin
Join Date: Jun 2009
Location: Italy
Posts: 69
Rep Power: 8
Marta is on a distinguished road
Hi All! I've just downloaded and compiled successfully IcoStructFoam for 1.6 OpenFoam version. The run didn't show any problems, but i have some difficulties in visualising the results with paraview.

What happens is that when i try to convert the format using foamToVTK command, it says it cannot find points file, which is obvious, since thi file is not inside the constant directory, but in region1 and region2 folders.
Do you have any ideas about how i can solve the problem for postprocessing?

Thanks a lot

Marta
Marta is offline   Reply With Quote

Old   December 24, 2010, 07:33
Default icostructFoam paraFoam
  #7
New Member
 
khaled
Join Date: Nov 2010
Posts: 11
Rep Power: 6
PeterII is on a distinguished road
Hi Waku, I've compiled icoStructFoam on OpenFOAM-1.7.1, but I have some problems with running the paraFoam utility for both of the regions of the thingOnStick test case. Could you please give me some advices about this?
PeterII is offline   Reply With Quote

Old   December 25, 2010, 21:28
Default
  #8
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi, khaled

I've just compiled and run the case on OF-1.7.1, successfully.
For visualization, paraFoam showed error such fv*** files not found in ./system directory, but I think this error is not so serious.

How about check two workarounds below:
1) Use formToVTK
> foamToVTK -region region1
> foamToVTK -region region2
> paraview (paraFoam)
Please use Warp filter with vector D to see region2 deformation.


2) merge fvSolution and fvSchemes files in two regions and set ./system directory
Sorry, I've not fully tested yet, but after just copying fv*** files in region1 into system directory, paraFoam works good.

Regards,
waku2005
waku2005 is offline   Reply With Quote

Old   December 27, 2010, 12:09
Default
  #9
New Member
 
khaled
Join Date: Nov 2010
Posts: 11
Rep Power: 6
PeterII is on a distinguished road
Hi Waku,
thank you so much for this information. I compiled with success the case thingOnStick but I cannot compile the second example icoStructFoamTest and I receive the following error:

--> FOAM FATAL IO ERROR:
patch type 'wall' not constraint type 'symmetryPlane'
for patch top of field D in file "/home/khaled/OpenFOAM/khaled-1.7.1/applications/solvers/IcoStructFoam_Rev561/icoStructFoamTest/0/region2/D"

file: /home/khaled/OpenFOAM/khaled-1.7.1/applications/solvers/IcoStructFoam_Rev561 /icoStructFoamTest/0/region2/D::boundaryField::top from line 43 to line 43.
From function symmetryFvPatchField<Type>::symmetryFvPatchField
(
const fvPatch& p,
const Field<Type>& field,
const dictionary& dict
)
in file fields/fvPatchFields/constraint/symmetry/symmetryFvPatchField.C at line 99.

FOAM exiting



please anny help

thank you in advance

best Regards
khaled
PeterII is offline   Reply With Quote

Old   December 27, 2010, 20:57
Default
  #10
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi khaled,

Error message show you that there is somthing wrong in "top" patch of region 2.
Edit constant/region2/polyMesh/boundary file like below and please run again.

> vi constant/region2/polyMesh/boundary

top
{
// type wall;
type symmetryPlane;
nFaces 60;
startFace 4760;
}

Hope this helps you.
best regards,
waku.
waku2005 is offline   Reply With Quote

Old   December 31, 2010, 19:01
Default icofsiopenfoam
  #11
New Member
 
khaled
Join Date: Nov 2010
Posts: 11
Rep Power: 6
PeterII is on a distinguished road
hello waku,

thanks for your very useful suggestions. I would like to do simulation of blood flow including fluid-structure interaction. I'm trying to install icoFsiFoam on OpenFOAM 1.7.1.
Thanks in advance for any suggestions
khaled
PeterII is offline   Reply With Quote

Old   October 17, 2012, 04:28
Default
  #12
New Member
 
fatemeh sadeghpour
Join Date: May 2012
Posts: 13
Rep Power: 5
fatemeh.sadeghpour is on a distinguished road
Hi khaled
has your problem solved?i cant run icoStructFoamTest in oF 2.1.1 , it has a lot off errors like:

when i use:
ln -s region1/polyMesh/ icoStructFoamTest/constant/polyMesh

it reply:


create time



--> FOAM FATAL ERROR:
Cannot open mesh description file
"/home/sadeghpour/OpenFOAM/sadeghpour-2.1.1/run/tutorials/icoStructFoamTest/constant/polyMesh/blockMeshDict"


From function blockMesh
in file blockMeshApp.C at line 148.

FOAM exiting


and when i use the command :
foamToVTK . icoStructFoamTest -mesh region1

it reply:

Create time

Create mesh region1 for time = 0



--> FOAM FATAL IO ERROR:
cannot open file

file: /home/sadeghpour/OpenFOAM/sadeghpour-2.1.1/run/tutorials/icoStructFoamTest/system/region1/fvSchemes at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 87.

FOAM exiting

please someone help me what should i do?
fatemeh.sadeghpour is offline   Reply With Quote

Old   February 14, 2013, 15:14
Default icoStructFoam, where is mesh.update() ??
  #13
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 129
Rep Power: 8
ziemowitzima is on a distinguished road
Hello,
I am working currently on a solver with moving mesh. Actually one boundary is changing according to solution of other equation.
I implemented mesh changes same as in icoSructFoam solver. But to update my mesh I need to add following line:
HTML Code:
 mesh.update();
which is not included in icoSructFoam...

I am confused because my mesh is not get updated without this command, but mesh in icoSructFoam is updated...

Moreover my solution, after mesh is updated, behaves somehow strange.

I hope anyone will have any suggestions.

ZMM
ziemowitzima is offline   Reply With Quote

Old   February 14, 2013, 18:27
Default
  #14
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 129
Rep Power: 8
ziemowitzima is on a distinguished road
Hi,

HTML Code:
mesh.update();
is not a problem anymore, I found the differences between my solver and icoStructFoam. Now it works both with and without "mesh.update();".

But I still have the problem with energy being added to the system after mesh is changed.
It happens for simple diffusion equation:
HTML Code:
fvScalarMatrix omEqn
        (
        dimDt*fvm::ddt(om)
       - fvm::laplacian(dimensionedScalar("1",dimensionSet(0, 2, 0, 0, 0),1), om)
         );        
        omEqn.solve();
if mesh does not change, or change only "a little" then it behaves as diff equation, namely field "om" diffuses out, and getting weaker.
But if mesh is changed "little more than a little" then "om" gets stronger and stronger with each time step...

Do you have any clue why it is like that ??

thanks
ZMM
ziemowitzima is offline   Reply With Quote

Reply

Tags
icofsifoam, openfoam 1.6

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
Installing OpenFOAM 1.6 on Ubuntu 10.4 sda OpenFOAM Installation 5 May 27, 2010 20:19
kOmegaSST in openfoam 1.6 Gearb0x OpenFOAM 2 March 3, 2010 07:02
Cross-compiling OpenFOAM 1.6 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 7 January 19, 2010 16:39
OpenFOAM 1.6 package - CentOS 5.3 x86_64 linnemann OpenFOAM Installation 36 September 2, 2009 12:46
OpenFOAM Version 1.6 Released opencfd OpenFOAM Announcements from ESI-OpenCFD 0 July 27, 2009 17:55


All times are GMT -4. The time now is 00:58.