twoWayFiles-coupling

Dear all,
I'm a newbie, have just installed the CFDEM package (and was able to run correctly the first testcase settlingTestMPI).

Now I'm trying to run the twoSpheresGlowinski testcase, and I get exactly the same error described in the above posts: the calculation stops right in the beginning with the message:

wait for file ".../run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"

I looked at the scripts and it seems to me that the DEMrun() function calls the "liggghts" executable. I have no such file. Is it correct? How can I proceed?

Thanks,
Michele

UPDATE:
I created a symlink to the executable generated in the liggghts src (lmp_opensuse-mpi).
Starting the two scripts separately produce the following output.

DEM part:
##############################################################
michele@hod:~/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski> ./DEMrun.sh
rm: cannot remove `/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/log_run_cfdemSolverIB_shared-twoSpheresGlowinski-DEM': No such file or directory

// run_cfdemSolverIB_shared-twoSpheresGlowinski-DEM //

/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/DEM

##############################################################

and the CFD part:
##############################################################
michele@hod:~/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski> ./CFDrun.sh
rm: cannot remove `/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/log_run_cfdemSolverIB_shared-twoSpheresGlowinski-CFD': No such file or directory

// run_cfdemSolverIB_shared-twoSpheresGlowinski-CFD //

/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD

rm: cannot remove `couplingFiles/*': No such file or directory
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.x-e04b9fe59a09
Exec : cfdemSolverIB_shared
Date : Feb 08 2012
Time : 16:14:21
Host : "hod.lem3"
PID : 15283
Case : /home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/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

Create mesh for time = 0

Selecting dynamicFvMesh dynamicRefineFvMesh
Reading field p

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

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

Creating dummy density field rho = 1

Reading field phiIB

Reading field phiIB

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model laminar
Reading field interFace

Reading g
Selecting locateModel standard
Selecting dataExchangeModel twoWayFiles
Selecting voidFractionModel bigParticle

W A R N I N G - do not use in combination with differentialRegion model!

W A R N I N G - this model does not yet work properly!

Selecting averagingModel dilute
Selecting regionModel allRegion
Selecting meshMotionModel noMeshMotion

cfdem version: cfdem-2.3.0-beta-2011-12-05-13:51

If BC are important, please provide volScalarFields -imp/expParticleForces-
Selecting forceModel ShirgaonkarIB
Selecting forceModel Archimedes
Selecting liggghtsCommandModel execute ,provide dicts, numbered from 0 to n
liggghtsCommand run 10
firstCouplingStep = 1
lastCouplingStep = 1000000000
couplingStepInterval = 1

Starting time loop

Time = 0.0003

Selected 0 cells for refinement out of 4000.
Selected 0 split points out of a possible 0.
Courant Number mean: 0 max: 0
- evolve()

timeStepFraction() = 1
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
nr particles = 2
particle 0 is represented by 4cells
impForces = 0,0,0
impForces = 0,0,0
writing to file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/dragforce1"...
evolve done.
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.20265e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 3.20265e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 2.30776e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.611424, Final residual = 5.26185e-07, No Iterations 49
time step continuity errors : sum local = 1.09426e-15, global = 1.01763e-32, cumulative = 1.01763e-32
DICPCG: Solving for p, Initial residual = 0.733362, Final residual = 5.18069e-07, No Iterations 48
time step continuity errors : sum local = 1.07303e-15, global = -3.7673e-32, cumulative = -2.74967e-32
DICPCG: Solving for p, Initial residual = 0.643078, Final residual = 5.40593e-07, No Iterations 47
time step continuity errors : sum local = 1.10635e-15, global = 4.57046e-32, cumulative = 1.82079e-32
DICPCG: Solving for p, Initial residual = 0.719117, Final residual = 6.60655e-07, No Iterations 47
time step continuity errors : sum local = 1.04992e-15, global = -4.48172e-33, cumulative = 1.37262e-32
particleCloud.calcVelocityCorrection()
DICPCG: Solving for phiIB, Initial residual = 1, Final residual = 7.47446e-07, No Iterations 37
ExecutionTime = 7.77 s ClockTime = 9 s

Time = 0.0006

Selected 75 cells for refinement out of 4000.
Refined from 4000 to 4525 cells.
--> FOAM Warning :
From function dynamicRefineFvMesh::refine(const labelList&)
in file dynamicRefineFvMesh/dynamicRefineFvMesh.C at line 309
Cannot find surfaceScalarField phi_0 in user-provided flux mapping table
1
(
phi U
)
The flux mapping table is used to recreate the flux on newly created faces.
Either add the entry if it is a flux or use (phi_0 none) to suppress this warning.
Selected 0 split points out of a possible 75.
Courant Number mean: 3.86059e-06 max: 0.00125842
- evolve()

timeStepFraction() = 1
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
nr particles = 2
particle 0 is represented by 20cells
impForces = -0.00842491,-0.00842491,0.398186
impForces = 0.058032,0.058032,0.0797293
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/dragforce1"
writing to file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/dragforce1"...
evolve done.
DILUPBiCG: Solving for Ux, Initial residual = 0.324026, Final residual = 1.43092e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.324026, Final residual = 1.43092e-06, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.180016, Final residual = 1.03403e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.729219, Final residual = 8.26253e-07, No Iterations 51
time step continuity errors : sum local = 6.56144e-12, global = -2.4339e-23, cumulative = -2.4339e-23
DICPCG: Solving for p, Initial residual = 0.0505155, Final residual = 7.14931e-07, No Iterations 45
time step continuity errors : sum local = 4.17057e-12, global = 5.60924e-23, cumulative = 3.17534e-23
DICPCG: Solving for p, Initial residual = 0.00208332, Final residual = 7.94382e-07, No Iterations 27
time step continuity errors : sum local = 5.6962e-12, global = 6.57733e-25, cumulative = 3.24111e-23
DICPCG: Solving for p, Initial residual = 0.000115827, Final residual = 9.09118e-07, No Iterations 12
time step continuity errors : sum local = 6.33984e-12, global = 5.30479e-24, cumulative = 3.77159e-23
particleCloud.calcVelocityCorrection()
DICPCG: Solving for phiIB, Initial residual = 0.601124, Final residual = 9.28262e-07, No Iterations 37
ExecutionTime = 7.88 s ClockTime = 9 s

Time = 0.0009

Selected 264 cells for refinement out of 4525.
Refined from 4525 to 6373 cells.
--> FOAM Warning :
From function dynamicRefineFvMesh::refine(const labelList&)
in file dynamicRefineFvMesh/dynamicRefineFvMesh.C at line 309
Cannot find surfaceScalarField phi_0 in user-provided flux mapping table
1
(
phi U
)
The flux mapping table is used to recreate the flux on newly created faces.
Either add the entry if it is a flux or use (phi_0 none) to suppress this warning.
Selected 5 split points out of a possible 291.
Unrefined from 6373 to 6338 cells.
--> FOAM Warning :
From function dynamicRefineFvMesh::refine(const labelList&)
in file dynamicRefineFvMesh/dynamicRefineFvMesh.C at line 526
Cannot find surfaceScalarField phi_0 in user-provided flux mapping table
1
(
phi U
)
The flux mapping table is used to recreate the flux on newly created faces.
Either add the entry if it is a flux or use (phi_0 none) to suppress this warning.
Courant Number mean: 1.99677e-05 max: 0.00500506
- evolve()

timeStepFraction() = 1
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/x1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/v1"
wait for file "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
reading from file: "/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/couplingFiles/radius1"
nr particles = 2
#0 Foam::error::printStack(Foam::Ostream&) in "/home/michele/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/home/michele/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#2 in "/lib64/libc.so.6"
#3 Foam::polyMesh::findTetFacePt(int, Foam::Vector const&, int&, int&) const in "/home/michele/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#4 Foam::polyMesh::pointInCell(Foam::Vector const&, int, Foam::polyMesh::cellRepresentation) const in "/home/michele/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#5 Foam::polyMesh::findCell(Foam::Vector const&, Foam::polyMesh::cellRepresentation) const in "/home/michele/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#6 Foam::standardSearch::findCell(double** const&, double**&, double**&) const in "/home/michele/OpenFOAM/michele-2.1.x/platforms/linux64Gcc46DPOpt/lib/liblagrangianCFDEM_shared.so"
#7 Foam::cfdemCloudIB::evolve() in "/home/michele/OpenFOAM/michele-2.1.x/platforms/linux64Gcc46DPOpt/lib/liblagrangianCFDEM_shared.so"
#8
in "/home/michele/OpenFOAM/michele-2.1.x/platforms/linux64Gcc46DPOpt/bin/cfdemSolverIB_shared"
#9 __libc_start_main in "/lib64/libc.so.6"
#10
at /home/abuild/rpmbuild/BUILD/glibc-2.14.1/csu/../sysdeps/x86_64/elf/start.S:116
GNU Octave, version 3.4.3
Copyright (C) 2011 John W. Eaton and others.
This is free software; see the source code for copying conditions.
There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. For details, type `warranty'.

Octave was configured for "x86_64-suse-linux-gnu".

Additional information about Octave is available at http://www.octave.org.

Please contribute if you find this software useful.
For more information, visit http://www.octave.org/help-wanted.html

Read http://www.octave.org/bugs.html to learn how to submit bug reports.

For information about changes from previous versions, type `news'.

error: load: file `../../DEM/post/position_particle_1.txt' seems to be empty!
error: load: unable to extract matrix size from file `../../DEM/post/position_particle_1.txt'
error: called from:
error: /home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/octave/postproc.m at line 11, column 5

** (evince:16561): WARNING **: Error stating file '/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/octave/pos_y_two_part_rec_glow.eps': No such file or directory

** (evince:16561): WARNING **: Error setting file metadata: No such file or directory

** (evince:16571): WARNING **: Error stating file '/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/octave/vel_y_two_part_rec_glow.eps': No such file or directory

** (evince:16571): WARNING **: Error setting file metadata: No such file or directory
deleting data at: /home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski
rm: cannot remove `/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/0.*': No such file or directory
rm: cannot remove `/home/michele/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski/CFD/VTK': No such file or directory
done
michele@hod:~/OpenFOAM/michele-2.1.x/run/cfdemSolverIBCase_shared/twoSpheresGlowinski>
##############################################################

Now the CFD solver starts running but crashes at timestep 0.0009 s.
I'm working with Openfoam 2.1.x (but have already installed both 1.7.x and 2.0.x). Should I try a previous OpenFOAM release? Which one is preferrable?

Thank you
Michele.