Simulation with cfdemSolverPiso

Submitted by pipegaldames on Tue, 05/07/2013 - 19:01

Hi all, i try to use the solver cfdemSolverPiso for my own case, for this I modified the ErgunTestMPI example.
When i run my case output is:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.x-6ff339325ac7
Exec : decomposePar
Date : May 07 2013
Time : 12:32:47
Host : "pipegaldames"
PID : 20074
Case : /home/felipegaldames/CFDEM/felipegaldames-PUBLIC-2.2.x/run/hidrociclon/CFD
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Decomposing mesh region0

Create mesh

Calculating distribution of cells
Selecting decompositionMethod simple

Finished decomposition in 0.03 s

Calculating original mesh data

Distributing cells to processors

Distributing faces to processors

Distributing points to processors

Constructing processor meshes

Processor 0
Number of cells = 9610
Number of faces shared with processor 1 = 620
Number of faces shared with processor 2 = 215
Number of faces shared with processor 3 = 7
Number of processor patches = 3
Number of processor faces = 842
Number of boundary faces = 3158

Processor 1
Number of cells = 3922
Number of faces shared with processor 0 = 620
Number of faces shared with processor 3 = 115
Number of processor patches = 2
Number of processor faces = 735
Number of boundary faces = 1569

Processor 2
Number of cells = 3922
Number of faces shared with processor 0 = 215
Number of faces shared with processor 3 = 241
Number of processor patches = 2
Number of processor faces = 456
Number of boundary faces = 1460

Processor 3
Number of cells = 9610
Number of faces shared with processor 0 = 7
Number of faces shared with processor 1 = 115
Number of faces shared with processor 2 = 241
Number of processor patches = 3
Number of processor faces = 363
Number of boundary faces = 3323

Number of processor faces = 1198
Max number of cells = 9610 (42.0337% above average 6766)
Max number of processor patches = 3 (20% above average 2.5)
Max number of faces between processors = 842 (40.5676% above average 599)

Time = 0

Processor 0: field transfer
Processor 1: field transfer
Processor 2: field transfer
Processor 3: field transfer

End.

// run_parallel_cfdemSolverPiso_ErgunTestMPI_CFDDEM //

/home/felipegaldames/CFDEM/felipegaldames-PUBLIC-2.2.x/run/hidrociclon/CFD

rm: cannot remove `couplingFiles/*': No such file or directory
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.x-6ff339325ac7
Exec : cfdemSolverPiso -parallel
Date : May 07 2013
Time : 12:32:49
Host : "pipegaldames"
PID : 20089
Case : /home/felipegaldames/CFDEM/felipegaldames-PUBLIC-2.2.x/run/hidrociclon/CFD
nProcs : 4
Slaves :
3
(
"pipegaldames.20090"
"pipegaldames.20091"
"pipegaldames.20092"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Reading field p

Reading physical velocity field U
Note: only if voidfraction at boundary is 1, U is superficial velocity!!!

Reading momentum exchange field Ksl

Reading voidfraction field voidfraction = (Vgas/Vparticle)

Creating dummy density field rho

Reading particle velocity field Us

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model laminar

Reading g
Selecting locateModel engine
Selecting dataExchangeModel twoWayMPI
Starting up LIGGGHTS for first time execution
Executing input script '../DEM/in.hidro'
LIGGGHTS (Version LIGGGHTS-PUBLIC 2.3, compiled 2013-03-13-11:04:35 by felipegaldames based on LAMMPS 20 Apr 2012)
Created orthogonal box = (-0.5 -0.3 -0.2) to (0.5 0.4 0.2)
2 by 2 by 1 MPI processor grid
ERROR on proc 0: Cannot open mesh file hidrociclon.stl (input_mesh_tri.cpp:78)
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 20089 on
node pipegaldames exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.x-6ff339325ac7
Exec : reconstructPar
Date : May 07 2013
Time : 12:32:49
Host : "pipegaldames"
PID : 20094
Case : /home/felipegaldames/CFDEM/felipegaldames-PUBLIC-2.2.x/run/hidrociclon/CFD
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

--> FOAM FATAL ERROR:
No times selected

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

FOAM exiting

simulation finished? ...press enter to proceed

**************************************************************

Apparently the error relates with the Liggghts code, specifically when insert the stl file:

fix cad all mesh/surface/stress file hidrociclon.stl type 1 wear finnie
********************************************************

This fix is not well implemented? , how should be?

Thanks in advance.

Maryam | Thu, 05/09/2013 - 08:01

The fix must point to the address of the stl file:
fix cad all mesh/surface/stress file ../DEM/hidrociclon.stl type 1 wear finnie

Maryam

pipegaldames | Thu, 05/09/2013 - 20:59

Tranks Marayam, it works, but now other error appears in this line:

fix ins nve_group insert/stream seed 5330 distributiontemplate pdd1 &
nparticles 10000 massrate 1 insert_every 1000 overlapcheck yes all_in no vel constant -5.0 0.0 0.0 &
nsertion_face inface

ERROR: Fix insert/stream (id ins): unknown keyword or wrong keyword order (fix_insert_stream.cpp:107)

I try run this fix in a simulation only with LIGGGHTS and works well, but in CFDEM appear this error.
How should be this line in CFDEM?

Regards,
Felipe

Maryam | Thu, 05/16/2013 - 05:24

If this is exactly what you have in your input file, there is a typo in the fix: "insertion_face"

Maryam