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

[ImmersedBoundary] Immersed Boundary Method in OpenFOAM-3.1-ext

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

Like Tree26Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 4, 2016, 12:21
Default
  #61
New Member
 
Join Date: Oct 2013
Posts: 12
Rep Power: 9
Nig̣ is on a distinguished road
I will be glad to do it, can you please give us some references to the cases and results that are you talking ?
I've done some search but I haven't found nothing.

Thank you in advance.
Nig̣ is offline   Reply With Quote

Old   March 6, 2016, 16:51
Default
  #62
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,957
Blog Entries: 45
Rep Power: 122
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
Quick answer: In regards to this:
Quote:
Originally Posted by hjasak View Post
We have published laminar and turbulent flow validation results at ERCOFTAC last year and did a validation study on a moving immersed boundary this December.

Why don't you start from those cases and find out where you made q mistake. All validation cases are public
The aforementioned cases are probably the ones available at one or both of the following repositories:
wyldckat is offline   Reply With Quote

Old   March 14, 2016, 18:36
Default
  #63
Member
 
crixman's Avatar
 
Christian
Join Date: Apr 2014
Posts: 74
Rep Power: 8
crixman is on a distinguished road
Hi everybody,
I am making some basic benchmarks of the IBM method.
I tried a simple laminar Hagen-Poseuille flow and the mass velocity of the flow is not exactly correct (20m/s analytical vs 17.5m/s from the simulation) and convergence after about 1000 iterations.
I cannot understand the problem, if someone is interested I can post the case.
Let me know,
Christian
crixman is offline   Reply With Quote

Old   March 20, 2016, 17:06
Default Errors runing potentialIbFoam in parallel
  #64
New Member
 
Robert Grandl
Join Date: Aug 2015
Posts: 16
Rep Power: 7
RobertG is on a distinguished road
Dear Prof. Hrvoje Jasak,
I would like to thank you, for all the effort you made by implementing the immersed boundary methode.
At the moment I'm using them more or less succesfully. I designed a 2D testcase, on which every thing worked fine.
But when I used it with a more complex meshshape in 3D, I discovered a some problems.

I managed to solve a couple of them and would like to publish their solution and discuss the problems, I did not manage to solve.

Solved problems are:
- immersedBoundaryFvPatch.C:
Running potentialIbFoam in parallel, it is possible that because of the different proceses are different fast. Which causes at some points a difference of more then 2 sec at the position of the programm on the different processes. Because of this delay, it is possible that a process freezes when it demands a value of another process, bevore it was send.
I solved it, by using MPI_Barrier( MPI_COMM_WORLD ); to synchronize the processes. It is neccessery to use this, before every up and download of a value (OPstream, IPstream and sendAndReceive(...)/immersedBoundaryFvPatchField.C)

- immersedBoundaryFvPatch.C:
cellProcCells[cellI].setSize(100) is at some cases to small and causes the potentialIbFoam to crush.
I would recommend resetting the size dynamicly if it is necessary.


Problems I was not able to solve.
- void immersedBoundaryFvPatchField<Type>::write(Ostream& os) const
It takes a lot of time to do U.write() (or so I belive. Once I waited 12 h for it, but it was not finished)
triValue() was calculated pretty fast, so I don't know what happens after it. Seriously, I have no clue what causes the extreme long run-time.
But I have one recommandation for improving the speed (if it would work regulary). If the writing-procedure is serialized, so that only one process writes at a time (It might improve the writing procedure when a hard drive is used).

When I use the potentialIbFoam in serial, it works fine. Non of the noticed problems to be seen.
But if I run decomposePar after it, it takes about one week to decompose the case.
This might be related to the last mentioned writing-problem.
But unfortunatelly, I don't finde every related class in this matter.

I hope you or somebody else can help me in this matter and give me a clue, how solve the last mentioned problem or give me a hint, in which classes are related to this problem.

Best Regards
RobertG
RobertG is offline   Reply With Quote

Old   March 21, 2016, 09:00
Default
  #65
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,863
Rep Power: 29
hjasak will become famous soon enough
Hi Robert,

Thank you for looking into this. Can you please send me the files you have changed to review the bug fixes and merge them in.

I am looking into the best way to rewrite the parallel communications in the Immersed Boundary to make everything faster, but this is a bigger project so it will require real funding. The comms we have now are fine for a smaller number of processors, but as the decomposition increases, the code will become inefficient. I will speak to Zeljko to see what we should do; however, some significant funding is required to get this fixed.

Regarding the writing of U and other fields: at this point in the code, the data from the IB points is interpolated into the points of the STL surface. The algorithm I am using is a surface walk, which is linear if the STL is clean and does not protrude outside of the background mesh. Ideally, the resolution of the STL surface should be comparable with the resolution of the background mesh.

I assume you found a case where the surface walk algorithm fails, eg. if the STL is folded over or similarly broken.

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

Old   April 30, 2016, 11:00
Default
  #66
Member
 
Francis
Join Date: Jan 2014
Location: Toronto
Posts: 50
Rep Power: 8
afrotimy is on a distinguished road
Quote:
Originally Posted by Phicau View Post
Hi Hua,

this error rings a bell to me:

Code:
[pjh-Veriton-D832:8321] *** An error occurred in MPI_Recv
[pjh-Veriton-D832:8321] *** on communicator MPI_COMM_WORLD
[pjh-Veriton-D832:8321] *** MPI_ERR_TRUNCATE: message truncated
[pjh-Veriton-D832:8321] *** MPI_ERRORS_ARE_FATAL: your MPI job will now abort
It was also triggered by a parallel boundary condition I developed. The solution was found thanks to Bernhard:

Edit $WM_PROJECT_DIR/etc/controlDict and change commsType to nonBlocking .

You should try and see if this helps with your problem.

Best,

Pablo
Hi foamers,

I can see that this thread has been for a long time.

I just want to know if anyone has been able to run icoIbFoam or icoIbDyMFoam in parallel in foam-extent-3.2.

I have the decomposeParDict correctly specified and ran well for other solvers, but when running the IB solver, it gives error messages are quits as it have been reported.

I tried to change the commsType to unblocking, problem still not resolved.

I noted that in foam-extend-3.2, the etc/controlDict was kept as etc/controlDict-Example. Although, I don't know why but I copied it as etc/controlDict.

Any suggestion would be appreciated.

Last edited by afrotimy; April 30, 2016 at 11:02. Reason: Typo
afrotimy is offline   Reply With Quote

Old   June 23, 2016, 13:05
Talking immersed boundary for CHT problem
  #67
New Member
 
XuDongNa
Join Date: Feb 2012
Location: Harbin Heilongjiang Province,China
Posts: 11
Rep Power: 10
supersoldier is on a distinguished road
Send a message via ICQ to supersoldier
I use openfoam-extend-3.2 version for studying the problem of IB. I have successfully run the case cavity which is from https://github.com/wyldckat/Immersed.../cavity/cavity.

I create a new case myself which is also a lid driven cavity, but twice height than before. I create a stl file, which is shown in left picture below in blue. When I run this case, I found the immersed boundary is not work.


I study the case and immersed boundary setting again, I find that the stl surface may be need closed. So I create another stl file. This time the IB works well.


I add a temperature equation in the solver icoIbFoam, I want to simulate the case like the conjugate heat transfer problem, the upper is fluid and the bottm part is solid , the heat is transfer from the fluid to the solid.


But after run the case I get confused myself about the immersed boundary :

The first is if the immersed boundary now implemented in OpenFoam-extend3.2 can handle the the multiregion problem, such as heat transfer between two regions? I think maybe now the IB just act as a wall boundary is that correct?

The second is if the stl surface need closed. Becase I want to add cyclic boundary for the left and right side of the rectangle, if the stl surface need closed, the solid part boundary will all be IB!

I am doing some work about solid materials combustion, so if you can give me some advices , I will really appreciate it. Thank you.
Attached Files
File Type: zip cavityMultiRegion.zip (40.9 KB, 15 views)
supersoldier is offline   Reply With Quote

Old   June 27, 2016, 15:21
Default IBM for Flexible Wall Problem
  #68
New Member
 
Join Date: Aug 2011
Posts: 9
Rep Power: 11
Engineers is on a distinguished road
Dear all,

I'm tasked with the following problem, but I'm not sure if I can get it done using just the IcoIbFoam as requested. Here is the problem, we need to simulate flow around a flexible finite aspect ratio flat plate positioned normal to the incoming flow using IBM. The idea is that the plate would deform under loadings due to the vortex shedding, and we need to capture how the wake transforms due to these deformations. I've read in various places here how the FSI solver accommodates flexible elastic walls, pipes, etc. However, I am not sure if the FSI can be incorporated into the IBM solver (IcoIbFoam). What are your thoughts on this?

This is what we have already done: we have simulated the flow around a rigid stationary finite aspect ratio using IcoIbFoam. We have also simulated the flow around a rigid and pitching plate using IcoIbFoam.

Thanks a ton

AH
Engineers is offline   Reply With Quote

Old   July 19, 2016, 13:33
Default
  #69
Senior Member
 
SinaJ
Join Date: Nov 2009
Posts: 136
Rep Power: 12
sina_mech is on a distinguished road
Hi All! I am new to immersed boundary in foam-ext 3.2 and have a basic question!

How do you create STL files?!! I create them using Rhinoceros, and make a before exporting them as SLT. Is that a right way? (please check the attached stl file for an example)

My other question is that how we determine the spacing between the points of the immersed boundary? Does it even matter for the current solver?

Sorry if the questions are so trivial! a newbie here

Thank you in advance!
Attached Files
File Type: zip test.zip (1.3 KB, 6 views)
sina_mech is offline   Reply With Quote

Old   July 21, 2016, 16:03
Default
  #70
Senior Member
 
SinaJ
Join Date: Nov 2009
Posts: 136
Rep Power: 12
sina_mech is on a distinguished road
Quote:
Originally Posted by sina_mech View Post
Hi All! I am new to immersed boundary in foam-ext 3.2 and have a basic question!

How do you create STL files?!! I create them using Rhinoceros, and make a before exporting them as SLT. Is that a right way? (please check the attached stl file for an example)

My other question is that how we determine the spacing between the points of the immersed boundary? Does it even matter for the current solver?

Sorry if the questions are so trivial! a newbie here

Thank you in advance!
Dear Prof. Hrvoje Jasak,

In addition to questions above, I have another question, and I would really appreciate if you could kindly advise me about these questions.

In the presentation files, it's mentioned that we can refine the volume mesh.
do we need to refine the IB boundary as well? and if yes, how can I do that?

I could edit the existing code to refine the volume mesh, but still have problem for flow around a cylinder, and I think it might be because of IB file. This is the error I get
Code:
Starting time loop

Time = 0.01



--> FOAM FATAL ERROR:
allPoints.size() < 5 : 4

    From function void immersedBoundaryFvPatch::makeInvDirichletMatrices()
    in file immersedBoundaryFvPatch/immersedBoundaryFvPatchLeastSquaresFit.C at line 95.

FOAM aborting

Thank you very much
Sina
sina_mech is offline   Reply With Quote

Old   October 22, 2016, 18:34
Default
  #71
New Member
 
anonymous
Join Date: Jun 2016
Posts: 1
Rep Power: 0
dmsg is on a distinguished road
Hello everyone. I'm pretty new to this forum and I'm new to IBM in OpenFOAM. I am having issues installing open foam-3.1-ext with my existing OpenFOAM installs (v3.0.1). Has anyone experienced this before?

I attached the output from ./Allwmake
Attached Files
File Type: txt makelog.txt (1.6 KB, 6 views)
dmsg is offline   Reply With Quote

Old   January 10, 2017, 05:51
Default ibSamplingPoints
  #72
Member
 
Join Date: May 2014
Posts: 40
Rep Power: 8
Phil_ is on a distinguished road
Hi foamers,

does anyone have a clue why 'sampling points' are used in the IBM, especially in the context of turbulence modelling with wall functions?

Best Regards
Philip
Phil_ is offline   Reply With Quote

Old   January 10, 2017, 06:05
Default
  #73
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,863
Rep Power: 29
hjasak will become famous soon enough
Hi,

Think about how the immersed boundary works. Having the IB fit in the wall function cell would be impossible if you do not know y+ - and this data cannot be taken out of the IB fitting itself because it is unstable.

I found an original paper from Standofrd doing this on a structured mesh (badly) and I developed + implemented the idea further. There should be papers describing this, eg ECCOMAS.

Hope this helps,

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

Old   January 29, 2017, 15:14
Default LES using IBM
  #74
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 448
Rep Power: 14
syavash is on a distinguished road
Hi,

I do LES on version 2.3 and I was wondering if it is possible in extended version to perform LES using IBM?! I mean does a pimpleIbFoam-like solver exist? or is it doable to develop such a solver just like simpleIbFoam has been devlopped?!

Thanks,
Syavash

P.S: I am trying to use LES on a rough surface with a body-fitted mesh but I am not sure it is the best idea because the mesh (especially aspect ratio) is not a typical LES mesh due to surface roughness
syavash is offline   Reply With Quote

Old   May 21, 2017, 10:59
Default pressure distribution
  #75
New Member
 
Saeed
Join Date: Oct 2016
Posts: 3
Rep Power: 6
saeed_hb is on a distinguished road
Hello
I am modeling an airfoil with icoDyMIbFoam and I want to get the pressure distribution over the airfoil, in fact I want to know the pressure on ib surface ... how can I do that? is there any function like sampleDict for Ib surfaces?? The ib boundary is not real, since it has mapped in the background mesh so I have not any real boundary that can be used in sampleDict as a boundary.... this is the problem here

Thanks
saeed_hb is offline   Reply With Quote

Old   May 22, 2017, 06:16
Default
  #76
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,863
Rep Power: 29
hjasak will become famous soon enough
See foam-extend-4.0/src/immersedBoundary/immersedBoundaryForce

There are function objects for this.
saeed_hb likes this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 25, 2017, 11:20
Default FSI with Ib method
  #77
New Member
 
Saeed
Join Date: Oct 2016
Posts: 3
Rep Power: 6
saeed_hb is on a distinguished road
Is there any way that immersed boundary method can be used for fluid structure interaction problems in foam-extend-4.0?
saeed_hb is offline   Reply With Quote

Old   November 19, 2017, 12:59
Default IB in OpenFoam
  #78
Member
 
Ben 017
Join Date: Nov 2017
Posts: 70
Rep Power: 4
Ben UWIHANGANYE is on a distinguished road
Hello Dear Foamers,

Please help to do tutorials of IB in OpenFoam
- Cylinder in a Cavity
1. How to Define volume mesh in constant/polyMesh/blockMeshDict dictionary?
2. How to Create polyMesh using blockMesh?

3. how to Include immersedBoundaryPolyPatch into polyMesh boundary.

4. How to Include immersedBoundaryFvPatchField into boundary field of pressure and velocity fields.

Through which way, i can make this?
and I have tried to run IcoIbFoam but said that such command doesn't exist

Thank you!
Ben UWIHANGANYE is offline   Reply With Quote

Old   November 19, 2017, 13:49
Default
  #79
New Member
 
Robert Grandl
Join Date: Aug 2015
Posts: 16
Rep Power: 7
RobertG is on a distinguished road
Quote:
Originally Posted by Ben UWIHANGANYE View Post
Hello Dear Foamers,

Please help to do tutorials of IB in OpenFoam
- Cylinder in a Cavity
1. How to Define volume mesh in constant/polyMesh/blockMeshDict dictionary?
2. How to Create polyMesh using blockMesh?

3. how to Include immersedBoundaryPolyPatch into polyMesh boundary.

4. How to Include immersedBoundaryFvPatchField into boundary field of pressure and velocity fields.

Through which way, i can make this?
and I have tried to run IcoIbFoam but said that such command doesn't exist

Thank you!

Hello Ben,
thous are very basic question.
Tharefore, I whould suggest, you should have a loock at,
https://www.openfoam.com/documentati.../blockMesh.php
This will provide a answer for question 1.+2.

regarding 3.+4.
may be you should habe a look at the tutorials you can find in:
foam-extend-4.0/tutorials/immersedBoundary

If want to you use immersedBoundary, you should use foam-extend-4.0

There are 3 reasons why icoIbFoam whount start.
1. Immersed Boundaries are officialy available since foam-extend-3.2
2. You did not compile foram-extend-3.2 or 4.0 correctly.
3. You did not start the script foram-extend-X.X/etc/bashrc
(4. You are not using foam-extend at all)

Best Regards
RobertG
RobertG is offline   Reply With Quote

Old   February 2, 2018, 03:39
Thumbs up Put my own algorithm in OpenFOAM-IBM
  #80
New Member
 
Tongwei-Zhang
Join Date: Jan 2018
Posts: 5
Rep Power: 4
fireztw is on a distinguished road
Quote:
Originally Posted by RobertG View Post
Hello Ben,
thous are very basic question.
Tharefore, I whould suggest, you should have a loock at,
https://www.openfoam.com/documentati.../blockMesh.php
This will provide a answer for question 1.+2.

regarding 3.+4.
may be you should habe a look at the tutorials you can find in:
foam-extend-4.0/tutorials/immersedBoundary

If want to you use immersedBoundary, you should use foam-extend-4.0

There are 3 reasons why icoIbFoam whount start.
1. Immersed Boundaries are officialy available since foam-extend-3.2
2. You did not compile foram-extend-3.2 or 4.0 correctly.
3. You did not start the script foram-extend-X.X/etc/bashrc
(4. You are not using foam-extend at all)

Best Regards
RobertG
Hello RobertG,
As a beginner of IBM in OF like Ben, I have run some tutorials in OF-ex3.2 and changed the mesh/IB(.stl) by my needs. But i can't understand the specific processing methods(eg: solver etc.),following questions:
1. What changes are there in the Ib solver compared with original (eg: icoFoam to icoIbFoam)
2. I want to use my own algorithm to solve relevant problems, what files are need to understanding and make changes.

Best Regards
fireztw
fireztw is offline   Reply With Quote

Reply

Tags
immersed boundary method, openfoam-1.6-ext, openfoam-3.1-ext

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
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
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 11:53
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
[ImmersedBoundary] who has the paper about immersed boundary method in openfoam blueshit OpenFOAM Community Contributions 1 November 18, 2013 08:16


All times are GMT -4. The time now is 18:40.