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

[ImmersedBoundary] Error while running IBM solver

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

Like Tree1Likes
  • 1 Post By pvpnrao

Reply
 
LinkBack Thread Tools Display Modes
Old   February 22, 2016, 22:27
Default Error while running IBM solver
  #1
Member
 
Narsimha Rao P V P
Join Date: Jan 2016
Posts: 33
Rep Power: 4
pvpnrao is on a distinguished road
I downloaded the IBM solvers and turtorials from openfoamwiki.
I am encountering the following error, while running my case

--> FOAM FATAL ERROR:
allPoints.size() < 9 : 7

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

FOAM aborting

I looked into the source code in immersedBoundaryFvPatchLeastSquaresFit.C line 91. I understand that the allpoints.size() should be grater than ncoeffs for the case to run. But I am having a hard time following the lineage to the definition of these variables.
Can somebody help.
ordinary likes this.

Last edited by pvpnrao; February 22, 2016 at 23:51.
pvpnrao is offline   Reply With Quote

Old   April 16, 2016, 12:07
Default
  #2
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
i am having the same problem. did you manage to fix it somehowe??

Last edited by urion; April 17, 2016 at 10:32.
urion is offline   Reply With Quote

Old   April 16, 2016, 13:30
Default
  #3
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick questions:
  1. Which exact installation instructions did you two follow?
  2. Which foam-extend version are you using?
__________________
wyldckat is offline   Reply With Quote

Old   April 16, 2016, 23:37
Default
  #4
Member
 
Narsimha Rao P V P
Join Date: Jan 2016
Posts: 33
Rep Power: 4
pvpnrao is on a distinguished road
Wyldckat

I used the instructions given on the following webpage

https://openfoamwiki.net/index.php/E...mersedBoundary

I am using Foam Extend 3.1. Foam 3.2 was not yet released when I started working with IBM.

Although, I did not try attempting to fix this issue yet. But, from my correspondence with Prof. Jasak, I understand, this is a mesh refinement problem, the immersed boundary is not being able to find a cell centroid for interpolation.

So, if I keep refining the mesh locally using cellSet+RefineMesh do you think I should be able to fix the problem.
pvpnrao is offline   Reply With Quote

Old   April 17, 2016, 10:48
Default
  #5
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
I don't think refining the mesh is the solution cause my case with a complex Ib runs perfectly on a coarse grid. The mentioned error only occurs after refining the mesh
urion is offline   Reply With Quote

Old   April 17, 2016, 17:54
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick answer: The best I can do on this topic is to suggest that you build from source code the "nextRelease" branch from foam-extend 3.2. It has several updates to IBM, including some related to dynamic meshing.

To get the correct branch, either clone directly:
Code:
git clone --branch nextRelease git://git.code.sf.net/p/foam-extend/foam-extend-3.2
Or after cloning, go into the folder "foam-extend-3.2" and run:
Code:
git checkout nextRelease
Use the latest instructions as a reference: http://openfoamwiki.net/index.php/In...oam-extend-3.2
wyldckat is offline   Reply With Quote

Old   April 19, 2016, 10:46
Default
  #7
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
Wyldckat
I now use foam-extend 3.2 but still get the same error. Any other suggestions?
urion is offline   Reply With Quote

Old   April 19, 2016, 18:30
Default
  #8
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quote:
Originally Posted by urion View Post
Any other suggestions?
Quick answer: I or anyone else will need more details on how to reproduce the same error that you are getting.
wyldckat is offline   Reply With Quote

Old   April 20, 2016, 04:17
Default
  #9
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
I'm running a case of turbulent channel flow where the IB consisting of a smooth plate represents the walls of the channel. this case is running well until i exchange the smooth plate IB with a complex structured IB. the new case is only running for a coarse grid. when i refine the mesh I'm getting the mentioned error. What further informations do you need?
I could send you the case then you could reproduce the error usingt the icoIbFoam solver though I'm using a modified version
urion is offline   Reply With Quote

Old   April 20, 2016, 17:47
Default
  #10
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quote:
Originally Posted by urion View Post
What further informations do you need?
Sorry, I should have pointed out that you should have followed the instructions on this thread: http://www.cfd-online.com/Forums/ope...-get-help.html

The idea is that without a test case, I can't reproduce the same error you're getting.
I could very easily spend 2h trying to try and reproduce the same error you're getting and always not get an error message .

edit: At the very least, I need detailed instructions on how to use an existing tutorial case and configure it to give the same crash you are getting.

Last edited by wyldckat; April 20, 2016 at 17:48. Reason: see "edit:"
wyldckat is offline   Reply With Quote

Old   April 21, 2016, 03:06
Default icoIbFoam channel flow with komplex IB, --> FOAM FATAL ERROR: allPoints.size() < 9
  #11
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
I modified the channelFoam tutorial to run the case with icoIbFoam-solver and a complex IB as channel walls. A sketch of the geometry can be found under the link:https://www.dropbox.com/s/euwbfw0b05...leyIB.jpg?dl=0
and a screendump of the mesh under the link:https://www.dropbox.com/s/nvbpnj04n7t29li/mesh.jpg?dl=0
The output of checkmesh can be found under the link: https://www.dropbox.com/s/smhvofhqnz...kMesh.log?dl=0
Link to fvSchemes: https://www.dropbox.com/s/ciaumy6j2h...fvSchemes?dl=0
Link to fvSolution: https://www.dropbox.com/s/a7z956vzm0...vSolution?dl=0
The output of the terminal is: https://www.dropbox.com/s/3nneljz2j5...bFoam.log?dl=0
And a link to the full case: https://www.dropbox.com/sh/d5bw0xi8l...eVZJILh1a?dl=0
note that if I reduce the mesh resulotion to the half and again add the IB-statement in the boundary file the case runs without the mentioned error
urion is offline   Reply With Quote

Old   May 1, 2016, 19:34
Default
  #12
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
@urion - Quick answer: I only managed to take a quick look into this the other day and I haven't managed to do any tests yet.

What comes to mind is that the problem might be related to the fact that your STL files are on the borders of the mesh, which might be the reason for this problem.
At least if I understood you correctly, the problem also occurs when the whole mesh is refined, which could lead to numerical issues.

Which reminds me: please try increasing the refinement in the STL files as well.
For example, in ParaView you can use the filter "Subdivide", which will divide each triangle on the STL into several smaller triangles. This might help with numerical issues, because each point on the mesh will be able to better "see" the nearby vertices of the triangles, therefore leading to smaller numerical errors... at least in theory.
wyldckat is offline   Reply With Quote

Old   May 2, 2016, 02:25
Default prevent pressure correction in IB
  #13
New Member
 
K.H
Join Date: Feb 2016
Posts: 15
Rep Power: 4
urion is on a distinguished road
hi wyldckat thanks for your answer. i will try refining the stl-file as well. I have another question regarding the solver I am using. I expanded the icoIbFoam solver in foam extend 3.1 by the prssure-correction made in the solver channelFoam so i can use the new solver as DNS-solver for turbulent channel flow with IB. The problem is that the pressure correction seems to be performed in the whole channel but i dont want to perform the pressure correction in the IB. what changes do i have to make to the solver especially the pressure correction to prevent this.?
Link to used solver:
https://www.dropbox.com/sh/wfgf194ul...u3_FkDEAa?dl=0
urion is offline   Reply With Quote

Old   December 6, 2016, 05:35
Default
  #14
New Member
 
Dong Zhang
Join Date: Dec 2016
Posts: 2
Rep Power: 0
what is on a distinguished road
you can try to exame your surface normals in your stl file, the surface normals of all the triangles should point outward.
what is offline   Reply With Quote

Old   June 25, 2017, 19:57
Default
  #15
New Member
 
Chung Shen
Join Date: Nov 2016
Posts: 1
Rep Power: 0
chung1shen is on a distinguished road
Dear All,

I installed foam-extend-4.0 on Ubuntu 16.04 and tried to run the tutorial movingCynlinderInChannelIco in parallel with 8 processors.

The only thing I modified is, in system/decomposeParDict:
numberOfSubdomains 8;
method simple;
simpleCoeffs
{
n (4 2 1);
delta 0.001;
}

It failed with the following messages:
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] allPoints.size() < 5 : 3
[1]
[1] From function void immersedBoundaryFvPatch::makeInvDirichletMatrices( )
[1] in file immersedBoundaryFvPatch/immersedBoundaryFvPatchLeastSquaresFit.C at line 95.
[1]
FOAM parallel run aborting
[1]

[5]
[5]
[5] --> FOAM FATAL ERROR:
[5] allPoints.size() < 5 : 4
[5]
[5] From function void immersedBoundaryFvPatch::makeInvDirichletMatrices( )
[5] in file immersedBoundaryFvPatch/immersedBoundaryFvPatchLeastSquaresFit.C at line 95.
[5]
FOAM parallel run aborting
[5]

How should I proceed from here? I tried subdividing the immersed boundary mesh and/or refining the domain mesh but still getting the same error.

Thanks in advance for your help.

Chung
chung1shen 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
Running external application in solver dariodario32 OpenFOAM Programming & Development 2 July 10, 2017 08:52
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
No running solver after successful installation of OF-2.1.x on SGI Altix UV 1000 sebastian OpenFOAM Installation 16 November 13, 2012 06:10
Statically Compiling OpenFOAM Issues herzfeldd OpenFOAM Installation 21 January 6, 2009 10:38
Compiling and Running UDFs on an IBM Power3 S. Hollman FLUENT 0 January 6, 2005 12:23


All times are GMT -4. The time now is 09:00.