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

Problem with icoFsiFoam in the latest OFdev

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 17, 2009, 22:24
Default Hi people ! I used to run 2
  #1
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 76
Rep Power: 8
mathieu is on a distinguished road
Hi people !

I used to run 2D fluid-structure interaction simulation with "closed" fluid-solid interface (see case A in the attachment). My current installation is OpenFOAM 1.4.1-dev revision 784.

Now, if I update my OpenFOAM installation to the latest 1.4.1-dev or even to the latest 1.5-dev version, those kind of simulation don't work anymore. However, if the fluid-solid interface is "open" (see case B in the attachment), then there is no problem at all (in all version).

The following package contains both cases A and B that can be run with the standard icoFsiFoam solver. More details are also provided in a pdf file.

Any help would be appreciated ! Kind regards,

Mathieu


mathieu is offline   Reply With Quote

Old   February 17, 2009, 22:29
Default http://www.cfd-online.com/Ope
  #2
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 76
Rep Power: 8
mathieu is on a distinguished road
package.tgz
mathieu is offline   Reply With Quote

Old   February 18, 2009, 07:47
Default Sorry, I do not get it. I hav
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,762
Rep Power: 21
hjasak will become famous soon enough
Sorry, I do not get it. I have just run caseA and I can see forces being transferred, a flow field, a structural deformation field and a moving mesh for both the fluid and the solid.

When you say "it does not work any more", what precisely does not work:
- the case does not run
- you get a floating point exception
- you do not get the flow field
- you do not get structural deformation
- you do not get mesh motion
- force transfer is not correct

I am missing something, but don't know what...

More info please,

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

Old   February 18, 2009, 08:36
Default Got it: this is to do with a s
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,762
Rep Power: 21
hjasak will become famous soon enough
Got it: this is to do with a search algorithm for a multiply projected boundary (thank Gavin and his FSI artery simulation).

Whta happens is that for a closed surface with escapes, you cannot guarantee the projection will be correct using the fast algorithm. I have changed the switch in

~/OpenFOAM/OpenFOAM-1.5-dev/etc/controlDict

to use

nSquaredProjection 1;

in the OptimisationSwitches section and all works well. The deformation on the solut surface looks a bit stupid (you may want to use mesh motion to deform that side as well), but at least we know what the problem (and solution) are.

Enjoy,

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

Old   February 18, 2009, 13:12
Default Dear Hrv, Thanks for the re
  #5
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 76
Rep Power: 8
mathieu is on a distinguished road
Dear Hrv,

Thanks for the reply ! I changed nSquaredProjection to 1 and the problem seems to be gone. Thanks again,

Mathieu
mathieu is offline   Reply With Quote

Old   February 19, 2009, 14:16
Default Dear Hrv, I investigated th
  #6
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 76
Rep Power: 8
mathieu is on a distinguished road
Dear Hrv,

I investigated this issue further and I found that there is still a problem with the mesh when using OF-1.5-dev with nSquaredProjection set to 1 (no more problems in the latest revision of 1.4.1-dev though!). See the right tip of the plate in the following pictures.

1.4.1-dev:


1.5-dev:


Is there something else I am missing or is this a bug ?

By the way, what did you mean by "solut surface" in your last message ?

Regards,

Mathieu
mathieu is offline   Reply With Quote

Old   February 19, 2009, 14:37
Default This is the kind of work I do
  #7
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,762
Rep Power: 21
hjasak will become famous soon enough
This is the kind of work I do under a support contract. Alternatively, you have a dig through the code yourself.

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

Old   September 23, 2010, 08:22
Default
  #8
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 80
Rep Power: 8
maruthamuthu_venkatraman is on a distinguished road
Quote:
Originally Posted by mathieu View Post
Dear Hrv,

I investigated this issue further and I found that there is still a problem with the mesh when using OF-1.5-dev with nSquaredProjection set to 1 (no more problems in the latest revision of 1.4.1-dev though!). See the right tip of the plate in the following pictures.

1.4.1-dev:


1.5-dev:


Is there something else I am missing or is this a bug ?

By the way, what did you mean by "solut surface" in your last message ?

Regards,

Mathieu

Hello Mathiew,
Have you fixed the error in 1.5-dev version for icoFSIFoam solver. As u stated 1.4-dev vrsion has no problems, Do I lack something if i use that version.

Regards
Muthu
maruthamuthu_venkatraman 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
Icofsifoam varun OpenFOAM Running, Solving & CFD 8 April 27, 2011 06:10
IcoFsiFoam vinz OpenFOAM Running, Solving & CFD 14 November 3, 2010 07:20
Errors in running a icoFsiFoam case jin_xu OpenFOAM Pre-Processing 0 June 9, 2008 06:48
Errors in compling icoFsiFoam jin_xu OpenFOAM Running, Solving & CFD 5 June 4, 2008 20:15
IcoFsiFoam tutorial pbo OpenFOAM Running, Solving & CFD 0 March 6, 2008 10:02


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