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

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

Register Blogs Community New Posts Updated Threads Search

Like Tree12Likes
  • 1 Post By Phicau
  • 2 Post By cfdivan
  • 2 Post By keepfit
  • 2 Post By keepfit
  • 4 Post By wyldckat
  • 1 Post By keepfit

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 5, 2013, 20:46
Default Is Immersed Boundary Method already available in OF 2.1.x or 2.2.x
  #1
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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, 09:41
Default
  #2
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 15
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, 16:24
Default
  #3
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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, 04:58
Default
  #4
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
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, 18:06
Default
  #5
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 15
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,
Linda and miladrakhsha like this.
cfdivan is offline   Reply With Quote

Old   October 24, 2013, 13:01
Default Fictitious domain method in OF
  #6
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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, 15:29
Default
  #7
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 15
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, 02:40
Default
  #8
Senior Member
 
Join Date: Jun 2011
Posts: 163
Rep Power: 15
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, 03:35
Default
  #9
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
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, 15:37
Default
  #10
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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 21:55.
keepfit is offline   Reply With Quote

Old   October 25, 2013, 15:42
Default
  #11
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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, 15:09
Default
  #12
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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, Pirlu and 1 others like this.
__________________
wyldckat is offline   Reply With Quote

Old   April 11, 2014, 07:36
Default
  #13
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
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 4, 2014, 00:40
Default
  #14
New Member
 
miladrakhsha
Join Date: Aug 2012
Posts: 29
Rep Power: 14
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: http://www.cfd-online.com/Forums/ope...m-3-1-ext.html

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

Old   April 1, 2016, 11:31
Default Fictitious domain solver
  #15
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21
msrinath80 is on a distinguished road
Is the fictitious domain solver (mentioned below) available somewhere to test?

Thanks!

Quote:
Originally Posted by keepfit View Post
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
msrinath80 is offline   Reply With Quote

Old   April 21, 2016, 19:21
Default
  #16
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14
keepfit is on a distinguished road
Quote:
Originally Posted by msrinath80 View Post
Is the fictitious domain solver (mentioned below) available somewhere to test?

Thanks!
Hi, I will release the solver after my Ph.D thesis is finalized.
frantov likes this.
keepfit is offline   Reply With Quote

Old   April 27, 2016, 08:24
Default
  #17
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 21
msrinath80 is on a distinguished road
Thanks David!
msrinath80 is offline   Reply With Quote

Reply

Tags
ib method, vof


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
3D Windturbine simulation in SU2 k.vimalakanthan SU2 15 October 12, 2023 06:53
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
implementation of the Immersed Boundary Method mi_cfd Main CFD Forum 19 April 24, 2019 02:24
[ImmersedBoundary] Immersed Boundary Method DoQuocVu OpenFOAM Community Contributions 11 July 13, 2018 10:36
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


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