CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Native Meshers: snappyHexMesh and Others (https://www.cfd-online.com/Forums/openfoam-meshing-snappyhexmesh/)
-   -   sHM - reconstructPar (https://www.cfd-online.com/Forums/openfoam-meshing-snappyhexmesh/81914-shm-reconstructpar.html)

Tobi November 10, 2010 09:12

sHM - reconstructPar
 
Hi all,

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

Code:

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 November 10, 2010 10:53

Why do you use -constant?

Have you tried it without?

Tobi November 10, 2010 11:16

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 January 28, 2011 09:43

Quote:

Originally Posted by Tobi (Post 282933)
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 January 28, 2011 10:04

Quote:

Originally Posted by sandy (Post 292657)
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

Code:

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

Code:

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 January 28, 2011 22:57

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 January 29, 2011 08:42

Quote:

Originally Posted by sandy (Post 292738)
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 January 19, 2016 09:05

Hello Tobi,

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

Quote:

--> 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 January 19, 2016 09:22

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 April 10, 2017 18:12

reconstructPar not reconstructing time 0
 
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:

Code:

#!/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

Code:

# 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 April 11, 2017 02:25

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

reconstructPar -withZero

Coke Rivas Ordenes April 11, 2017 07:45

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.


All times are GMT -4. The time now is 03:24.