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

1.6.x -> snappyHexMesh in parallel fails ...

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 31, 2009, 23:18
Default 1.6.x -> snappyHexMesh in parallel fails ...
  #1
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
Hi,

I'm experiencing some difficulties when trying to run snappyHexMesh with
mpirun and 2 cpus / I have tried the same case under 2 different machines
and have the same error. Note that the same case was working well under
1.5.x. Here is the message :

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Read mesh in = 0.11 s
[1]

[1]
[1] Cannot find file "" in directory "constant/triSurface"
[1]
[1] From function Time::findInstance(const fileName&, const word&, const IOobject::readOption)
[1] in file db/Time/findInstance.C at line 148.Overall mesh bounding box : (-22.5 -40 -0.0001) (22.5 40 18)

[1]
FOAM parallel run exiting
[1] Relative tolerance : 1e-08

Absolute matching distance : 9.3536109e-07

[0]
[0]
[0] Cannot find file "" in directory "constant/triSurface"
[0]
[0] From function Time::findInstance(const fileName&, const word&, const IOobject::readOption)
[0] in file db/Time/findInstance.C at line 148.
[0]
FOAM parallel run exiting
[0]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.


Best regards,

PO
podallaire is offline   Reply With Quote

Old   August 1, 2009, 05:07
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
As a workaround softlink or copy the constant/triSurface directory to the processor directories for now.
mattijs is offline   Reply With Quote

Old   August 1, 2009, 08:43
Default
  #3
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
ok, yes, the softlink works
podallaire is offline   Reply With Quote

Old   August 1, 2009, 09:56
Default
  #4
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
Hi,

another strange thing is that I cannot use decomposePar on a blockMesh which
was refined with refineMesh / I never had this problem before.

Regards,

PO
podallaire is offline   Reply With Quote

Old   August 1, 2009, 18:49
Default
  #5
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
sorry not to start a new thread for this / since I'm discussing meshing "challenges", I guess it's appropriate to re-use this thread ...

I believe that db/Time/timeSelector.C is causing problems whith checkMesh, it does not find the "constant" folder when no time steps are available.
podallaire is offline   Reply With Quote

Old   August 3, 2009, 03:53
Default
  #6
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Just create a dummy time directory; mkdir 0.
mattijs is offline   Reply With Quote

Old   August 3, 2009, 08:22
Default
  #7
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
Yes, that is what I tried yesterday.

Also, it appears that my problems when it comes to decompose a refined
mesh are caused by the "-overwrite" flag; I found out that not using this flag
and moving the new mesh in constant/polyMesh does not cause problems.

I also have issues with the same flag when using snappyHexMesh, that is why,
for now, I will not use it and will always move the polyMesh folder.

PO
podallaire is offline   Reply With Quote

Old   August 4, 2009, 00:42
Default Assign flow fields to decomposed mesh?
  #8
Member
 
Cem Albukrek
Join Date: Mar 2009
Posts: 50
Rep Power: 8
albcem is on a distinguished road
Relevant to "snappyhexmesh in parallel" topic stream, is there a way to assign the flow fields to the decomposed mesh that was built using snappyhexmesh in parallel?

Essentially, I am working with a 32 bit machine and I am pushing the OS memory limit on the size of mesh by using snappyhexmesh in parallel. After the decomposed mesh is built, I need to assign the flow fields to it, but obviously I am not able to first reconstruct the decomposed mesh.

Any ideas on how I can tackle the issue, without having to buy a 64bit machine?

Thanks.

Cem
albcem is offline   Reply With Quote

Old   August 4, 2009, 04:17
Default
  #9
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Regarding the -overwrite : remove the polyMesh/pointLevel, cellLevel and refinementHistory files. They should be consistent with the snappyHexMesh generated mesh so if you've regenerated a blockMesh they are no longer valid. I'll push a check into 1.6.x so at least you'll get an error message.
mattijs is offline   Reply With Quote

Old   August 4, 2009, 04:23
Default
  #10
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Quote:
Originally Posted by albcem View Post
Relevant to "snappyhexmesh in parallel" topic stream, is there a way to assign the flow fields to the decomposed mesh that was built using snappyhexmesh in parallel?

Essentially, I am working with a 32 bit machine and I am pushing the OS memory limit on the size of mesh by using snappyhexmesh in parallel. After the decomposed mesh is built, I need to assign the flow fields to it, but obviously I am not able to first reconstruct the decomposed mesh.

Any ideas on how I can tackle the issue, without having to buy a 64bit machine?

Thanks.

Cem
Please use this forum for bug reports only. Look in the documentation or on the main forum.
mattijs is offline   Reply With Quote

Old   August 4, 2009, 11:59
Default
  #11
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Quote:
Originally Posted by podallaire View Post
[1] Cannot find file "" in directory "constant/triSurface"
[1]
[1] From function Time::findInstance(const fileName&, const word&, const IOobject::readOption)

PO
I pushed a fix to 1.6.x so it now runs as before without copying the stl files.
Thanks for reporting.

Mattijs
mattijs is offline   Reply With Quote

Old   August 4, 2009, 13:57
Default
  #12
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 185
Rep Power: 8
podallaire is on a distinguished road
Thanks, I will give it a try today

cheers

PO
podallaire 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
SnappyHexMesh in Parallel bastil OpenFOAM Mesh Utilities 22 April 7, 2010 11:48
Parallel case setup boundry conditions snappyhexmesh oskar OpenFOAM Pre-Processing 5 September 11, 2009 01:12
interDyMFoam fails in parallel nikos_fb16 OpenFOAM Running, Solving & CFD 2 March 28, 2009 13:07
Ignition fails in parallel run combustion solvers msha OpenFOAM Bugs 17 January 17, 2009 04:49
Serial run OK parallel one fails r2d2 OpenFOAM Running, Solving & CFD 2 November 16, 2005 13:44


All times are GMT -4. The time now is 20:43.