CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   AMR and parallel run (http://www.cfd-online.com/Forums/openfoam/98280-amr-parallel-run.html)

lichmaster March 7, 2012 06:36

[Fixed] AMR and parallel run
 
Hi all,

I'm trying to run the damBreakWithObstacle tutorial, which solves the classical dam break problem with automatic mesh refinement. No problems with the serial run.
It seems also to work in parallel mode (decomposed with scotch).
However, when i try to run reconstructPar, it gives

Quote:

--> FOAM FATAL IO ERROR:
cannot find file

file: /home/lorenzo/Scrivania/damBreakWithObstacle/processor0/0.02/polyMesh/pointProcAddressing at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting
Any ideas?

nimasam March 7, 2012 08:49

first you should use "reconstructParMesh" then use "reconstructPar"
i saw the following command somewhere in forum which help me:
Quote:

timeList=$(ls processor0/ | awk '$1!="constant"');

for line in $timeList
do
echo "Recontructing mesh for t = "$line" s";
reconstructParMesh -time $line > logTmp;

echo "Reconstructing fields for t = "$line" s";
reconstructPar -time $line > logTmp;

echo "";
done

rm logTmp;

lichmaster March 8, 2012 04:10

Thanks a lot, had to do some minor changes to the code you posted but works great :)

kumar March 9, 2012 04:28

information regarding AMR
 
Hello guys,
I would like to know some basic things about the AMR in Openfoam.
In which version of OF is the AMR available.
And what kind of meshes can be used along with AMR. Any documentation or tutorial about the use .
Thanks for your answers.

bye
regards
K.Suresh kumar

lichmaster March 9, 2012 04:45

Hi,

I don't know in which version of OpenFOAM AMR first became operative, in this (2.1.0) surely it is. Maybe even in 2.0.0.

For now I have just tried to follow the tutorial in tutorial/multiphase/interDymFoam/ras/damBreakWithObstacle, and works flawlessly.

Given OF's capabilities to operate with generic polyhedral meshes, AMR should work with any kind of mesh.

Hrushi April 15, 2012 06:20

Quote:

Originally Posted by nimasam (Post 348132)
first you should use "reconstructParMesh" then use "reconstructPar"
i saw the following command somewhere in forum which help me:


Code:

timeList=$(ls processor0/ | awk '$1!="constant"');

for line in $timeList
do
echo "Recontructing mesh for t = "$line" s";
reconstructParMesh -time $line > logTmp;

echo "Reconstructing fields for t = "$line" s";
reconstructPar -time $line > logTmp;

echo "";
done

rm logTmp;

Hi Nima,

I ran the above script. It reconstructs the mesh but while reconstructing the fields it gives error. I am posting the output result of the script.

Recontructing mesh for t = 0 s


--> FOAM FATAL ERROR:
Your time was specified as 0 but there is no polyMesh/points in that time.
(there is a points file in "constant")
Please rerun with the correct time specified (through the -constant, -time or -latestTime (at your option).


From function reconstructParMesh
in file reconstructParMesh.C at line 449.

FOAM exiting

Reconstructing fields for t = 0 s


--> FOAM FATAL ERROR:
No times selected

From function reconstructPar
in file reconstructPar.C at line 139.

FOAM exiting


Recontructing mesh for t = 0.005 s
Reconstructing fields for t = 0.005 s
#0 Foam::error::printStack(Foam::Ostream&) in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::objectRegistry::checkOut(Foam::regIOobject&) const in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#4 Foam::regIOobject::~regIOobject() in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#5 Foam::IOList<int>::~IOList() in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#6 Foam::processorMeshes::read() in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libreconstruct.so"
#7 Foam::processorMeshes::readUpdate() in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libreconstruct.so"
#8 main in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/reconstructPar"
#9 __libc_start_main in "/lib64/libc.so.6"
#10 _start in "/Storage1/cfd/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/reconstructPar"
./recd: line 4: 5337 Segmentation fault (core dumped) reconstructPar -time $line >logTmp

Recontructing mesh for t = 0.01 s


After running the script, I have polyMesh directory in all time dumps in my case directory. But there are no fields.

What can be the reason for this?

Hrushikesh

MadsR May 3, 2012 07:23

Actually I begin to see the same error all of a sudden. Normally the reconstructParMesh hasn't been necessary and it doesn't even help here.

Mads


All times are GMT -4. The time now is 15:35.