Bug in reconstructPar

david · January 9, 2009, 08:14

Hi everybody

I found a bug in the reconstruction of decomposed cases with moving (not deforming) meshes. The new positions of the mesh points are reconstructed correctly if reconstructPar is used without any optional parameters. When the parameters -time or -latestTime are used, the points are not being reconstructed. The problem occurs with OpenFOAM-1.4.1 and 1.5.

Regards
David

see also:
http://www.cfd-online.com/OpenFOAM_D...es/1/9644.html

mattijs · January 12, 2009, 07:40

Thanks. Attached a fix. Could you let us know whether it works in general and I'll push it into 1.5.x.

reconstructPar.C

david · January 12, 2009, 11:52

Hi Mattijs

The fix is working. Thanks a lot.

Regards
David

feymark · April 22, 2009, 05:51

Hi everyone!

I'm having problem using reconstructPar -latestTime on a moving mesh. The problem is that the polyMesh-directory is not reconstructed. reconstructPar -time TIME works aslong as I don't use the latest time as TIME. I have tried the latest version of reconstructPar from the svn and also the code Mattijs posted here, it did not work.

Regards,
Feymark

mattijs · April 23, 2009, 14:35

Hi Feymark,

I've pushed a change to 1.5.x which worked on my testcase. Could you let me know if there are any problems with it?

Thanks,

Mattijs

feymark · April 27, 2009, 08:33

It works, thanks a lot!

/Feymark

feymark · April 27, 2009, 10:50

Hi,

I have another problem regarding but decomposePar this time. When I try to run decomposePar the polyMesh-dir isn't decomposed.

Best regards,
Feymark

feymark · April 27, 2009, 11:02

I should add that I have, before doing reconstructPar and removing the processor*-folders, 2 different polyMesh-folders under each processor*/-folder. After performing decomposePar the polyMesh in processor*/constant/ is there but not the one in processor*/TIME.

Regards,
Feymark

mattijs · April 28, 2009, 05:44

Does it work with -time? (so directly point at the time-with-polyMesh)

feymark · April 28, 2009, 06:50

reconstructPar works perfectly fine now. This is a problem with decomposePar and -time is not a valid option there. The polyMesh-folder that is suppose to be in the processor*/TIME-folder is not there after decomposing, only the polyMesh in processor*/constant.

/Feymark

mattijs · May 26, 2009, 13:11

Hi Feymark,

I pushed a few small changes into 1.5.x which should fix the problem for decomposing cases with meshes in time directories. If the original mesh (at the specified time) originates from anything other than 'constant' it will write the mesh to the time directory. E.g. original case has
- full mesh in constant
- polyMesh/points in 0.3
and current time is 0.5 it will create a full mesh in processor*/0.5/polyMesh and not create a constant folder.

Thanks for reporting.

January 9, 2009, 08:14	Hi everybody I found a bug	#1
david Member David Hora Join Date: Mar 2009 Location: Zürich, Switzerland Posts: 63 Rep Power: 17	Hi everybody I found a bug in the reconstruction of decomposed cases with moving (not deforming) meshes. The new positions of the mesh points are reconstructed correctly if reconstructPar is used without any optional parameters. When the parameters -time or -latestTime are used, the points are not being reconstructed. The problem occurs with OpenFOAM-1.4.1 and 1.5. Regards David see also: http://www.cfd-online.com/OpenFOAM_D...es/1/9644.html

January 12, 2009, 07:40	Thanks. Attached a fix. Could	#2
mattijs Senior Member Mattijs Janssens Join Date: Mar 2009 Posts: 1,419 Rep Power: 26	Thanks. Attached a fix. Could you let us know whether it works in general and I'll push it into 1.5.x. reconstructPar.C

January 12, 2009, 11:52	Hi Mattijs The fix is worki	#3
david Member David Hora Join Date: Mar 2009 Location: Zürich, Switzerland Posts: 63 Rep Power: 17	Hi Mattijs The fix is working. Thanks a lot. Regards David

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
ReconstructPar	kupiainen	OpenFOAM Post-Processing	12	April 21, 2009 02:47
[mesh manipulation] Error with reconstructPar	skabilan	OpenFOAM Meshing & Mesh Conversion	3	June 10, 2008 19:34
ReconstructPar	maka	OpenFOAM Bugs	6	August 22, 2007 05:23
Problem with reconstructPar	fabianpk	OpenFOAM	5	August 14, 2007 10:17
ReconstructPar error	r2d2	OpenFOAM Pre-Processing	7	September 12, 2006 10:37

April 22, 2009, 05:51		#4
feymark New Member Andreas Feymark Join Date: Mar 2009 Location: Gothenburg, Sweden Posts: 18 Rep Power: 17	Hi everyone! I'm having problem using reconstructPar -latestTime on a moving mesh. The problem is that the polyMesh-directory is not reconstructed. reconstructPar -time TIME works aslong as I don't use the latest time as TIME. I have tried the latest version of reconstructPar from the svn and also the code Mattijs posted here, it did not work. Regards, Feymark

April 23, 2009, 14:35		#5
mattijs Senior Member Mattijs Janssens Join Date: Mar 2009 Posts: 1,419 Rep Power: 26	Hi Feymark, I've pushed a change to 1.5.x which worked on my testcase. Could you let me know if there are any problems with it? Thanks, Mattijs

April 27, 2009, 08:33		#6
feymark New Member Andreas Feymark Join Date: Mar 2009 Location: Gothenburg, Sweden Posts: 18 Rep Power: 17	It works, thanks a lot! /Feymark

April 27, 2009, 10:50		#7
feymark New Member Andreas Feymark Join Date: Mar 2009 Location: Gothenburg, Sweden Posts: 18 Rep Power: 17	Hi, I have another problem regarding but decomposePar this time. When I try to run decomposePar the polyMesh-dir isn't decomposed. Best regards, Feymark

April 27, 2009, 11:02		#8
feymark New Member Andreas Feymark Join Date: Mar 2009 Location: Gothenburg, Sweden Posts: 18 Rep Power: 17	I should add that I have, before doing reconstructPar and removing the processor-folders, 2 different polyMesh-folders under each processor/-folder. After performing decomposePar the polyMesh in processor/constant/ is there but not the one in processor/TIME. Regards, Feymark

April 28, 2009, 05:44		#9
mattijs Senior Member Mattijs Janssens Join Date: Mar 2009 Posts: 1,419 Rep Power: 26	Does it work with -time? (so directly point at the time-with-polyMesh)

April 28, 2009, 06:50		#10
feymark New Member Andreas Feymark Join Date: Mar 2009 Location: Gothenburg, Sweden Posts: 18 Rep Power: 17	reconstructPar works perfectly fine now. This is a problem with decomposePar and -time is not a valid option there. The polyMesh-folder that is suppose to be in the processor/TIME-folder is not there after decomposing, only the polyMesh in processor/constant. /Feymark

May 26, 2009, 13:11		#11
mattijs Senior Member Mattijs Janssens Join Date: Mar 2009 Posts: 1,419 Rep Power: 26	Hi Feymark, I pushed a few small changes into 1.5.x which should fix the problem for decomposing cases with meshes in time directories. If the original mesh (at the specified time) originates from anything other than 'constant' it will write the mesh to the time directory. E.g. original case has - full mesh in constant - polyMesh/points in 0.3 and current time is 0.5 it will create a full mesh in processor*/0.5/polyMesh and not create a constant folder. Thanks for reporting.