CFD Online Discussion Forums - Allrun for SHM and LTSInterFoam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Allrun for SHM and LTSInterFoam (https://www.cfd-online.com/Forums/openfoam-solving/109369-allrun-shm-ltsinterfoam.html)

Allrun for SHM and LTSInterFoam

hi everybody
i try to run a calculation with this Allrun copy and modify from what i found on tutorial and on web
can some one tell me what doing this 2 lines?

find . -type f -iname "*level*" -exec rm {} \;
ls -d processor* | xargs -i cp -r 0.org/* ./{}/0/ $1

here the allrun script;

Code:

#!/bin/sh

cd ${0%/*} || exit 1    # run from this directory

# Source tutorial run functions

. $WM_PROJECT_DIR/bin/tools/RunFunctions

runApplication blockMesh

runApplication topoSet -dict cRefine.topoSetDict

runApplication refineMesh -overwrite -dict

rm log.refineMesh > /dev/null 2>&1

rm log.topoSet > /dev/null 2>&1

runApplication topoSet -dict cRerefine.topoSetDict

runApplication refineMesh -overwrite -dict



cp system/decomposeParDict-nonPar system/decomposeParDict

runApplication  decomposePar 

cp system/decomposeParDict-par system/decomposeParDict

runParallel snappyHexMesh 4  -overwrite



runApplication reconstructParMesh -constant

runApplication reconstructPar  # -constant -time 0

cp 0/alpha1.org 0/alpha1

runApplication setFields

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

ls -d processor* | xargs -i cp -r 0.org/* ./{}/0/ $1



runApplication setFields

rm log.decomposePar > /dev/null 2>&1

runParallel changeDictionary 4   -literalRE # Set boundary conditions

runApplication  decomposePar

runParallel renumberMesh 4 -overwrite #slight improvement in speed

runParallel  `getApplication` 4

runApplication reconstructParMesh -constant

runApplication reconstructPar

# ----------------------------------------------------------------- end-of-file# #

i fill to be nearly there, after losing some hairs, the error message in log.LTSInterfoam is;
[0] --> FOAM FATAL IO ERROR:
[0] size 679545 is not equal to the given value of 169885
[0]
[0] file: /mnt/disk2/laurent/boats/calculOF/laurent-2.1.x/hulls/a_04_refinemesh_parallel/processor0/0/alpha1 from line 18 to line 681403.
[0]

all suggestion or remarques will be very well come!
thanks you by advance for help!
LL

hello Foamers,
to make my question simple; if i want to run SHM and solver in parallel, do i need to reconstructPar, between SHM and the solver....
thanks you very much by advance

My guess is that you have still some processoN folder between after thw snappy reconstruct and the decomposition before the LTS and this might cause issues in cell number. I'd try to remove them in between.

Dear Laurent,
Regarding the two command lines:

Quote:

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

looks for all files with "level" in their name (regardless what comes before or after) and deletes these files

Quote:

ls -d processor* | xargs -i cp -r 0.org/* ./{}/0/ $1

copies all files and folders in 0.org/ to the 0 directory of each processor directory. Why the $1 is added at this end is not clear to me.

This second command seems pretty odd to me. DecomposePar will add the processor boundary conditions to each of the variables when decomposing but this command will overwrite the files again with the original 'undecomposed' ones. I expect this to be the cause of the error you get.

Besides this, the order of commands in your script also seems a bit odd to me. This is the order I would follow:

1. set up basic mesh (e.g. with blockMesh)
2. if necessary, refine/adjust mesh (e.g. refineMesh, snappyHexMesh)
3. set initial fields (e.g. copy 0.orig fields to 0, setFields, ...)
4. decomposePar
5. run the simulation (e.g. foamJob -p LTSinterFoam)
6. reconstructPar
7. Postprocessing ...

The message is, you should not make to many manipulations between decomposePar, running the simulation and reconstructPar (I would recommend not to make any). Try to modify your script according to this order and let me know if it solved the problem.

Regards,

L

hi;
first thanks you very much for answering, i just get the solver running!!! its very exciting!!!
in fact i saw that setFields is pretty smart: that set fields even when those are decomposed on multi processor!!!
and when i re-decompose, to run the solver, the repertories are re-write, so no need to delete them....
happy day to all LL