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

Is Immersed Boundary Method already available in OF 2.1.x or 2.2.x

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

Like Tree9Likes
  • 1 Post By Phicau
  • 1 Post By cfdivan
  • 2 Post By keepfit
  • 2 Post By keepfit
  • 3 Post By wyldckat

Reply
 
LinkBack Thread Tools Display Modes
Old   October 5, 2013, 19:46
Default Is Immersed Boundary Method already available in OF 2.1.x or 2.2.x
  #1
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
Hi Foamers,

I just wonder, is IB (Immersed Boundary) method already integrated in OF 2.2.x (or OF 2.1.x), or just available with 1.6-ext version?

Also I am interested in VOF + IB method, can someone recommend related journal papers on the multiphase flow (VOF or similar) with IB method?

Thanks a lot.


Best regards,

David
keepfit is offline   Reply With Quote

Old   October 6, 2013, 08:41
Default
  #2
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi David,

AFAIK the IBM method is not available on "standard" OF versions.

This post could help u out - meshless

Regards,
cfdivan is offline   Reply With Quote

Old   October 6, 2013, 15:24
Default
  #3
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
Quote:
Originally Posted by cfdivan View Post
Hi David,

AFAIK the IBM method is not available on "standard" OF versions.

This post could help u out - meshless

Regards,
Can we use the source code of IB method, and compile within OF2.x.x. Since I dont want to install extra 1.6-ext.
keepfit is offline   Reply With Quote

Old   October 7, 2013, 03:58
Default
  #4
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: London, UK
Posts: 242
Rep Power: 7
Phicau is on a distinguished road
Hi,

of course it is possible. It may take longer than implementing it from scratch, though, as there are thousands of changes to be made. OpenFOAM has changed a lot from version 1.6 to 2.x : some methods have disappeared and have been substituted by others, some of them are now called using other syntax, and so on.

Good luck!

Pablo
miladrakhsha likes this.
Phicau is offline   Reply With Quote

Old   October 11, 2013, 17:06
Default
  #5
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi keepfit,

It will be nice if you succeed compiling an working version of IBM for 2.x. I already tried, but then I just realized that I was not able to sort out all the problems.

I am looking forward for your success story !

Regards,
miladrakhsha likes this.
cfdivan is offline   Reply With Quote

Old   October 24, 2013, 12:01
Default Fictitious domain method in OF
  #6
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
Hi Foamers,

Instead of implementing IB method in OF 2.1.x/2.2.x, a so called fictitious domain (FD) method is used to present solid phase in fluid domain.

Compared to the original IB method in OF-1.6-ext that inside the immersed boundary the p/U/alpha1 field is void and the BC is fixed, the fictitious domain (FD) method treats the immersed body as fluid, plus a rigid constraint (D(U)=0) to force the "fluid" inside IB boundary behaviour as rigid body.

This new solver currently works fine in general (needs improvement), here is some rough results:

A solid ball drops from air into water

1. Initial status



2. ball dropping into water



3. air-water interface after ball dropped into water



4. Fictitious domain U field (ball is slowed down in water due to the drag force)




Cheers,

David
Mojtaba.a and Pirlu like this.

Last edited by keepfit; November 18, 2013 at 22:04.
keepfit is offline   Reply With Quote

Old   October 24, 2013, 14:29
Default
  #7
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi keepfit,

Pretty interesting approach! Have you tried to model external flows with this methodology...flow around a profile or sphere. It would be interesting to have a close view of the mesh at the wall/fluid interface.

Regards,
cfdivan is offline   Reply With Quote

Old   October 25, 2013, 01:40
Default
  #8
Senior Member
 
Join Date: Jun 2011
Posts: 141
Rep Power: 6
mechy is on a distinguished road
Hi
I can not find IBM in OF-1.6-ext
can you give me its position in application folder ?

Regards
mechy is offline   Reply With Quote

Old   October 25, 2013, 02:35
Default
  #9
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: London, UK
Posts: 242
Rep Power: 7
Phicau is on a distinguished road
Hey David,

it certainly looks very good, well done! Just a question, since you are dealing with a 2-phase flow, can you confirm that the mass (VOF) is conserved with this implementation?
Phicau is offline   Reply With Quote

Old   October 25, 2013, 14:37
Default
  #10
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
Quote:
Originally Posted by Phicau View Post
Hey David,

it certainly looks very good, well done! Just a question, since you are dealing with a 2-phase flow, can you confirm that the mass (VOF) is conserved with this implementation?
well, since the solid phase is treated as liquid (alpha1 is forced to 1), the alphaEqn is solved on the whole domain, the mass conservation can be satisfied in general.

However, a little "liquid" is observed leaking from within the IB boundary (only happens when solid phase in gas phase) . The reason is that the constraint is applied onto the solid phase after UEqn and pEqn is solved. It means that during the sub time-stepping (solving p and U), there is no CONSTRAINT on the solid phase. To avoid this problem the constraint might be applied in every sub time-steps, but it's computational expensive.

If someone has good advice, please feel free contact me.

Last edited by keepfit; October 25, 2013 at 20:55.
keepfit is offline   Reply With Quote

Old   October 25, 2013, 14:42
Default
  #11
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
Quote:
Originally Posted by mechy View Post
Hi
I can not find IBM in OF-1.6-ext
can you give me its position in application folder ?

Regards
there you are:

slides:

Source codes of IB method in OF-1.6-ext

Note that the IB method in OF-1.6-ext uses IBMask (0 - inside IB boundary, 1 - fluid domain) to calculate various fields (e.g. p. U), hence physical behaviour inside the immersed body is not solved.
wyldckat and chery1986 like this.
keepfit is offline   Reply With Quote

Old   April 5, 2014, 14:09
Default
  #12
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,530
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

Thank you David for posting the direct links to those files!

I've finally come around to git'ifying this source code and creating a wiki page for the Immersed Boundary toolkit:
Anyone and everyone, please feel free to edit the wiki page and fork the repository!

Best regards,
Bruno
Toorop, johnkh and Pirlu like this.
wyldckat is online now   Reply With Quote

Old   April 11, 2014, 06:36
Default
  #13
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 55
Rep Power: 5
keepfit is on a distinguished road
good job, wyldckat!

Quote:
Originally Posted by wyldckat View Post
Greetings to all!

Thank you David for posting the direct links to those files!

I've finally come around to git'ifying this source code and creating a wiki page for the Immersed Boundary toolkit:
Anyone and everyone, please feel free to edit the wiki page and fork the repository!

Best regards,
Bruno
keepfit is offline   Reply With Quote

Old   September 3, 2014, 23:40
Default
  #14
New Member
 
miladrakhsha
Join Date: Aug 2012
Posts: 27
Rep Power: 5
miladrakhsha is on a distinguished road
Hi all,
I just wanted to ask if there has been any update in regard to implementation of IBM
in official openFoam or at least in OpenFOAM 3.0-ext.
I am working with the one that has been implemented in OpemFOAM 1.6-ext but I am not satisfied with the solver and stability of the solver at all.
The solver terminates when bounding k, or omega warning appears. And, I realized that bounded schemes are not available in of1.6-ext.
Any comments would be greatly appreciated.

-----------
Moderator note: This topic is discussed further here: Immersed Boundary Method in OpenFOAM-3.1-ext

Last edited by wyldckat; September 14, 2014 at 05:38. Reason: See "Moderator note"
miladrakhsha is offline   Reply With Quote

Reply

Tags
ib method, vof

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
implementation of the Immersed Boundary Method mi_cfd Main CFD Forum 14 June 1, 2015 22:59
Immersed boundary method (IBM) vs Dynamic meshing tomyangbath Main CFD Forum 5 April 17, 2009 03:02
Immersed boundary method with momentum forcing Fuka Main CFD Forum 1 February 6, 2006 04:46
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55
New topic on same subject - Flow around race car Tudor Miron CFX 15 April 2, 2004 06:18


All times are GMT -4. The time now is 13:07.