|
[Sponsors] |
[ImmersedBoundary] Is Immersed Boundary Method already available in OF 2.1.x or 2.2.x |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Member
Join Date: Jun 2011
Posts: 51
Rep Power: 14 ![]() |
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, |
|
![]() |
![]() |
![]() |
![]() |
#6 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
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 Last edited by keepfit; November 18, 2013 at 21:04. |
|
![]() |
![]() |
![]() |
![]() |
#7 |
Member
Join Date: Jun 2011
Posts: 51
Rep Power: 14 ![]() |
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, |
|
![]() |
![]() |
![]() |
![]() |
#8 |
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 13 ![]() |
Hi
I can not find IBM in OF-1.6-ext can you give me its position in application folder ? Regards |
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18 ![]() |
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? |
|
![]() |
![]() |
![]() |
![]() |
#10 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
Quote:
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. |
||
![]() |
![]() |
![]() |
![]() |
#11 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
Quote:
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. |
||
![]() |
![]() |
![]() |
![]() |
#12 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,969
Blog Entries: 45
Rep Power: 127 ![]() ![]() ![]() ![]() ![]() ![]() |
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:
![]() Best regards, Bruno
__________________
|
|
![]() |
![]() |
![]() |
![]() |
#13 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
good job, wyldckat!
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#14 |
New Member
miladrakhsha
Join Date: Aug 2012
Posts: 29
Rep Power: 12 ![]() |
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 05:38. Reason: See "Moderator note" |
|
![]() |
![]() |
![]() |
![]() |
#15 | |
Senior Member
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 20 ![]() |
Is the fictitious domain solver (mentioned below) available somewhere to test?
Thanks! Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#16 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13 ![]() |
||
![]() |
![]() |
![]() |
![]() |
#17 |
Senior Member
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 20 ![]() |
Thanks David!
|
|
![]() |
![]() |
![]() |
Tags |
ib method, vof |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Radiation in semi-transparent media with surface-to-surface model? | mpeppels | CFX | 11 | August 22, 2019 07:30 |
implementation of the Immersed Boundary Method | mi_cfd | Main CFD Forum | 19 | April 24, 2019 01:24 |
3D Windturbine simulation in SU2 | k.vimalakanthan | SU2 | 14 | February 8, 2019 14:43 |
[ImmersedBoundary] Immersed Boundary Method | DoQuocVu | OpenFOAM Community Contributions | 11 | July 13, 2018 09:36 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 15:55 |