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

finite area method doesn't work in parallel (foam-extend-4.0)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 6, 2017, 09:35
Default finite area method doesn't work in parallel (foam-extend-4.0)
  #1
New Member
 
Andrea
Join Date: Jan 2016
Location: Italy
Posts: 4
Rep Power: 3
andlog is on a distinguished road
Hi all,

I've noted a problem with finite area solvers in parallel using foam-extend-4.0.

I copied both liquidFilmFoam and surfactantFoam tutorials, adding in the system folder the decomposeParDict from incompressible/simpleFoam/pitzDaily tutorial.

After blockMesh and makeFaMesh, I've decompose the case using decomposePar without a hitch.

Finally, I've tried to run the case in parallel typing "runParallel surfactantFoam 4 &", getting this error:


Code:
/*---------------------------------------------------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | foam-extend: Open Source CFD                    |

|  \\    /   O peration     | Version:     4.0                                |

|   \\  /    A nd           | Web:         http://www.foam-extend.org         |

|    \\/     M anipulation  | For copyright notice see file Copyright         |

\*---------------------------------------------------------------------------*/

Build    : 4.0-4fad65ce7cac

Exec     : surfactantFoam -parallel

Date     : May 02 2017

Time     : 13:03:33

Host     : ManDybHig_clust

PID      : 5971

CtrlDict : "/stoc/andrea/OpenFOAM_EXT/andrea-4.0/tutorial/FAM/planTransport/system/controlDict"

Case     : /stoc/andrea/OpenFOAM_EXT/andrea-4.0/tutorial/FAM/planTransport

nProcs   : 4

Slaves : 

3

(

ManDybHig_clust.5972

ManDybHig_clust.5973

ManDybHig_clust.5974

)



Pstream initialized with:

    nProcsSimpleSum   : 0

    commsType         : nonBlocking

SigFpe   : Enabling floating point exception trapping (FOAM_SIGFPE).



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time



Create mesh for time = 0



[0] 

[0] 

[0] --> FOAM FATAL ERROR: 

[0] otherProcNo:2 : illegal communicator -1650115904

Communicator should be within range 0..0

[0] 

[0]     From function PstreamGlobals::checkCommunicator(const label, const label)

[0]     in file db/IOstreams/Pstreams/PstreamGlobals.C at line 72.

[0] 

FOAM parallel run aborting

[0] 

[1] 

[1] 

[1] --> FOAM FATAL ERROR: 

[1] otherProcNo:0 : illegal communicator -1104586048

Communicator should be within range 0..0

[1] 

[1]     From function PstreamGlobals::checkCommunicator(const label, const label)

[1]     in file db/IOstreams/Pstreams/PstreamGlobals.C at line 72.

[1] 

FOAM parallel run aborting

[1] 

[2] 

[2] 

[2] --> FOAM FATAL ERROR: 

[2] otherProcNo:0 : illegal communicator -886773056

Communicator should be within range 0..0

[2] 

[2]     From function PstreamGlobals::checkCommunicator(const label, const label)

[2]     in file db/IOstreams/Pstreams/PstreamGlobals.C at line 72.

[2] 

FOAM parallel run aborting

[2] 

[3] 

[3] 

[3] --> FOAM FATAL ERROR: 

[3] otherProcNo:2 : illegal communicator 1469405888

Communicator should be within range 0..0

[3] 

[3]     From function PstreamGlobals::checkCommunicator(const label, const label)

[3]     in file db/IOstreams/Pstreams/PstreamGlobals.C at line 72.

[3] 

FOAM parallel run aborting

[3] 

--------------------------------------------------------------------------

MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD 

with errorcode 1.



NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.

You may or may not see output from other processes, depending on

exactly when Open MPI kills them.

--------------------------------------------------------------------------
I get the same error with all the finite-area-method tutorials. I haven't found any suggestion for this.


The above procedure successfully run using foam-extend-3.2 instead.
Therefore, it should be a bug related to the last version (foam-extend-4.0).

Any idea what this is and how to fix it?

Best,

Andrea
andlog is offline   Reply With Quote

Old   May 6, 2017, 20:10
Default
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,802
Rep Power: 24
hjasak will become famous soon enough
Let me have a look and get back to you - this looks like some recent development.

Apologies,

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

Old   May 8, 2017, 06:49
Default
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,802
Rep Power: 24
hjasak will become famous soon enough
Hi Andrea,

My bad: simple bug in comms initialisation. I have fixed it now; if you need it urgetnylu, please drop me a line.

The bug fix will be pushed asap.

Apologies,

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

Old   May 8, 2017, 09:08
Default
  #4
New Member
 
Andrea
Join Date: Jan 2016
Location: Italy
Posts: 4
Rep Power: 3
andlog is on a distinguished road
Hi Hrvoje,

Thank you so much for your prompt reply.
I' m currently in a development stage, so I can wait for the fix.

Thank you again,

Andrea
andlog is offline   Reply With Quote

Old   July 16, 2017, 20:32
Default
  #5
New Member
 
Join Date: Jul 2017
Posts: 3
Rep Power: 2
mseka is on a distinguished road
Hi Andrea,

I am wondering whether you are working on the finite area or not? I am interested in applying finite area on a problem but there is not enough documentation on how to modify the solver for other applications in a proper way, even the original Ph.D. thesis Z. Tukoviīc, Finite volume method on domains of varying shape, is not available online.

I am wondering whether you have any documents or would like to discuss a bit on the problem?
mseka is offline   Reply With Quote

Old   July 21, 2017, 19:38
Default
  #6
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,802
Rep Power: 24
hjasak will become famous soon enough
Hi,

Sorry, I didn't realise the PhD Thesis by Tukovic is not online. We will fix this asap- it should be on

http://www.fsb.hr/cfd

under Dissemination.

There is also a paper in Coumputers and Fluids:

A moving mesh finite volume interface tracking method for surface tension dominated interfacial fluid flow
Ž Tuković, H Jasak - Computers & fluids, 2012

https://scholar.google.co.uk/citatio...J:hqOjcs7Dif8C
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   July 21, 2017, 20:02
Default
  #7
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,802
Rep Power: 24
hjasak will become famous soon enough
My bad: here's the link to prof. Zeljko Tukovic PhD Thesis

http://foam-extend.fsb.hr/wp-content...c_PhD_2005.pdf
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Reply

Tags
bug, finite area method, foam extend 4.0, parallel

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
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 83 July 22, 2017 08:26
error with reactingFoam BakedAlmonds OpenFOAM Running, Solving & CFD 4 June 22, 2016 02:21
Using starToFoam clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04
mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16
Import problem ARC Open Source Meshers: Gmsh, Netgen, CGNS, ... 0 February 27, 2010 11:56


All times are GMT -4. The time now is 08:14.