[Tobi]

Hi all,

i run snappy with 6 cores. After it, i will reconstruct it again. But how?!

tobi@tobi:~/OpenFOAM/tobi-1.7.x/run/teg$ reconstructPar -constant
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.7.x-2154baf2ac24
Exec   : reconstructPar -constant
Date   : Nov 10 2010
Time   : 14:10:51
Host   : tobi
PID    : 2687
Case   : /home/tobi/OpenFOAM/tobi-1.7.x/run/teg
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0



--> FOAM FATAL IO ERROR: 
cannot open file

file: /home/tobi/OpenFOAM/tobi-1.7.x/run/teg/processor0/0/polyMesh/pointProcAddressing at line 0.

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

FOAM exiting

the file "pointProcAddressing" is missing.
Any ideas.
Tobi

[stevenvanharen]

Why do you use -constant?

Have you tried it without?

[Tobi]

Hi,

yes i ve tried it with all tags.

but i realized, that i need "reconstructParMesh" instead of reconstructPar after mesh generation. I am so stupid

tobi

[sandy]

Originally Posted by Tobi

Hi,

yes i ve tried it with all tags.

but i realized, that i need "reconstructParMesh" instead of reconstructPar after mesh generation. I am so stupid

tobi

Hi Tobi, what is the difference between reconstrucParMesh and reconstructPar?

[Tobi]

Originally Posted by sandy

Hi Tobi, what is the difference between reconstrucParMesh and reconstructPar?

hi Sandy,...

if you split your mesh to solve sHM with more cores you have to reconstruct your mesh with the command

reconstructMeshPar

if you run your equation with more cores you have to use after finishing

reconstructPar

Couse you have to reconstruct the mesh and all other files like (p / U / T / etc.)

Hope it 's helpful.

-> reconstructPar for solver
-> reconstructParMesh for sHM

I hope its correct - but i am certain sure about it.

Tobi

[sandy]

Hi Tobi, thank you for your reply. However, now I use lines step by step as follows:

1. blockMesh
2. snappyHexMesh
3. decomposePar
4. mpirun -np 8 interFoam -parallel
.....

Could you tell me which one I should use, reconstructPar or reconstructParMesh? Waiting for your help again.

Sandy

[Tobi]

Originally Posted by sandy

Hi Tobi, thank you for your reply. However, now I use lines step by step as follows:

1. blockMesh
2. snappyHexMesh
3. decomposePar
4. mpirun -np 8 interFoam -parallel
.....

Could you tell me which one I should use, reconstructPar or reconstructParMesh? Waiting for your help again.

Sandy

hi

5. reconstructPar

I added a bash file you can execute it with ./solve and change it - you should use your solver

http://www.file-upload.net/download-...solve.zip.html

Tobi

[stephie]

Hello Tobi,

maybe you can help me. I got the same error message:

--> FOAM FATAL IO ERROR:
cannot find file

file: /home/stephanie/Schreibtisch/FFW_E2/processor0/0.001/polyMesh/pointProcAddressing at line 0.

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

FOAM exiting

I am using a prepared SHM. Fore my case I am using interDyMFoam with an adaptiveMesh at the phase boundary.

I run these commands:

setFields >> log.setFields
decomposePar >> log.decomposePar
mpirun -np 2 interDyMFoam -case /home/stephanie/Schreibtisch/FFW_E2/ -parallel >> log.interDyMFoam
reconstructPar >> log.reconstructPar

There seemed to be no problem until the command reconstructPar. There are both processor-folders with the timesteps inside.

Might you have any idea, where I did a mistake?
I would be very grateful for your help.

Thank you and best regards,
Stephie

[TobM]

Hi Stephie,

with adaptive mesh refinement you have each time step a different mesh.
Use
1. reconstructParMesh for the mesh for each time step
2. reconstructPar for the fields

Best regards

Tobias

[Coke Rivas Ordenes]

Hi Tobi and Stephie,

Hope you can help me out with something that has been keeping my mind busy, related with reconstructing the case after running in it parallel.

The base case has the directories: 0.orig (with the respective fields), constant and system. Pretty standar. I run the simulation (PISO LES solver), decomposing the domain and reconstructing it, with the following Allrun script:

 #!/bin/sh
cd ${0%/*} || exit 1    # Run from this directory

# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions

# Make dummy 0 directory
mkdir 0

runApplication blockMesh
cp system/decomposeParDict.hierarchical system/decomposeParDict
runApplication decomposePar

cp system/decomposeParDict.ptscotch system/decomposeParDict
runParallel snappyHexMesh -overwrite -parallel

find . -type f -iname "*level*" -exec rm {} \;

ls -d processor* | xargs -I {} rm -rf ./{}/0
ls -d processor* | xargs -I {} cp -r 0.orig ./{}/0 $1

runApplication funkySetFields -case processor0 -time 0
mv log.funkySetFields log.funkySetFields1
runApplication funkySetFields -case processor1 -time 0
mv log.funkySetFields log.funkySetFields2
runApplication funkySetFields -case processor2 -time 0
mv log.funkySetFields log.funkySetFields3
runApplication funkySetFields -case processor3 -time 0
mv log.funkySetFields log.funkySetFields4

runParallel pisoFoam

runApplication reconstructParMesh -constant -mergeTol 1e-6

runApplication reconstructPar

It all goes splendid (with some issues, but the overall outcome is positive), but i have encounter a problem.

reconstructPar (or reconstructParMesh, not sure which of both) does not reconstruct time 0 (initial values), but it does it perfectly with the following time step.

I have already deleted the lines

# Make dummy 0 directory
mkdir 0

and re-run it all over again, but the problem persist; time 0 is not reconsrtucted.

Any thoughts? Thanks in advance for the time you might take reading my question.

[TobM]

reconstructPar has the option -withZero, use it like this:

reconstructPar -withZero

[Coke Rivas Ordenes]

Thank you very much! I tried -time 0 before asking, but i didn't encounter (or thought) the right way to use it so i received nothing but error messages . Thank you again, i will try it out very soon.

[cunconbkhp]

Hi tobi! I've run my case in oepnFoam 7 with the order below :

runApplication surfaceFeatures
runApplication blockMesh
runApplication decomposePar -copyZero
runParallel snappyHexMesh -overwrite -parallel
runParallel setRefinedFields

to recompose the case, i've trie reconstrucParmesh and reconstrucPar, event for each time step with a script :

#!/bin/bash
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;

but i have alaway the same error : processor0/0/polyMesh/pointProcAddressing not exist.

Could you please give me a hint!