(CFDEM Tutorial) Only error occurs on 'cfdemPostproc/fillCylinder'

Thank you for your advice and I make 'post' folder in DEM part. After the calculation, the DEM section works well, but the part of CFD still doesn't work at all (the same error before).

Then, I saw the 'log_run_ligghts_fillCylinder_DEM' first (just normal I think),
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------
// run_liggghts_fillCylinder_DEM //

/home/bc/CFDEM/bc-PUBLIC-2.3.x/run/tutorials/cfdemPostproc/fillCylinder/DEM

LIGGGHTS (Version LIGGGHTS-PUBLIC 3.1.0, compiled 2015-03-24-10:14:51 by bc, git commit b3030a8fd1f722c77f69150c7b2d24b386eab848 based on LAMMPS 23 Nov 2013)
atom_style granular
atom_modify map array
communicate single vel yes

boundary m m m
newton off

units si

region reg block -0.015 0.015 -0.015 0.015 -0.001 0.0554 units box
create_box 1 reg
Created orthogonal box = (-0.015 -0.015 -0.001) to (0.015 0.015 0.0554)
1 by 1 by 1 MPI processor grid

neighbor 0.001 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5

#pair style
pair_style gran model hertz tangential history
pair_coeff * *

#timestep, gravity
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#walls
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.0553
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.0138 0. 0.

#cfd coupling
#fix cfd all couple/cfd

#particle insertion
region bc cylinder z 0.0 0.0 0.012 0.02 0.03 units box
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 200 radius constant 0.001
fix pdd1 all particledistribution/discrete 1. 1 pts1 1.0
Fix particledistribution/discrete (id pdd1): distribution based on mass%:
pts1: d=2.000000e-03 (max. bounding sphere) mass%=100.000000%
Fix particledistribution/discrete (id pdd1): distribution based on number%:
pts1: d=2.000000e-03 (max. bounding sphere) number%=100.000000%

fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. -1. insert_every once overlapcheck yes all_in yes particles_in_region 1000 region bc

#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere

#screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

#insert the first particles so that dump is not empty
run 1
Setting up run ...
INFO: Resetting random generator for region bc
INFO: Particle insertion ins: inserting every 0 steps
Memory usage per processor = 4.32245 Mbytes
Step Atoms KinEng 1 Volume
0 0 -0 0 5.076e-05
insertion: proc 0 at 0 %
insertion: proc 0 at 10 %
insertion: proc 0 at 20 %
insertion: proc 0 at 30 %
INFO: Particle insertion ins: inserted 305 particle templates (mass 0.000256) at step 1
- a total of 305 particle templates (mass 0.000256) inserted so far.
WARNING: Particle insertion: Less insertions than requested (../fix_insert.cpp:674)
1 305 0.00012777063 0 5.076e-05
Loop time of 0.0215509 on 1 procs for 1 steps with 305 atoms

Pair time (%) = 1.40667e-05 (0.065272)
Neigh time (%) = 0.000242949 (1.12732)
Comm time (%) = 2.86102e-06 (0.0132757)
Outpt time (%) = 1.09673e-05 (0.05089)
Other time (%) = 0.0212801 (98.7432)

Nlocal: 305 ave 305 max 305 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 1060 ave 1060 max 1060 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1060
Ave neighs/atom = 3.47541
Neighbor list builds = 1
Dangerous builds = 0
dump dmp all custom 1000 post/dump.liggghts_init id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

run 3000 upto
Setting up run ...
Memory usage per processor = 7.99072 Mbytes
Step Atoms KinEng 1 Volume
1 305 0.00012777063 0 5.076e-05
1000 305 0.00015403997 0 5.076e-05
2000 305 7.416768e-05 6.7679775e-06 5.076e-05
3000 305 2.0997362e-06 2.1338041e-06 5.076e-05
Loop time of 0.0945849 on 1 procs for 2999 steps with 305 atoms

Pair time (%) = 0.0474048 (50.1187)
Neigh time (%) = 0.0137155 (14.5007)
Comm time (%) = 0.000310898 (0.328697)
Outpt time (%) = 0.00357103 (3.77548)
Other time (%) = 0.0295827 (31.2764)

Nlocal: 305 ave 305 max 305 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 1229 ave 1229 max 1229 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1229
Ave neighs/atom = 4.02951
Neighbor list builds = 57
Dangerous builds = 0
#write_restart liggghts.restart
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------

However, the CFD part still has some problems although I make the folder named 'post'.
When I look at the 'Allrun.sh' file, the folder '0' is related with the post-processing, but I'm not professional to OpenFOAM or LIGGGHTS, so I don't know what to do.

........(omit)....
#- clean up case
rm -r $casePath/CFD/0.*
........(omit).....

Also, Here is the 'run_cfdemPostproc_fillCylinder_CFD' files.
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------
// run_cfdemPostproc_fillCylinder_CFD //

/home/bc/CFDEM/bc-PUBLIC-2.3.x/run/tutorials/cfdemPostproc/fillCylinder/CFD

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.x-6f86ad2099b6
Exec : cfdemPostproc
Date : Apr 07 2015
Time : 16:53:45
Host : "ubuntu"
PID : 12111
Case : /home/bc/CFDEM/bc-PUBLIC-2.3.x/run/tutorials/cfdemPostproc/fillCylinder/CFD
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing 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 voidfraction field voidfraction = (Vgas/Vparticle)

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
Selecting locateModel standard
Selecting dataExchangeModel oneWayVTK
relativePath_"../DEM/post"
Selecting IOModel basicIO
--> FOAM Warning : IO model is in serial write mode, only data on proc 0 is written
Selecting probeModel off
Selecting voidFractionModel divided
Selecting averagingModel dense
Selecting clockModel off
start clock measurement at t >0.001
Selecting smoothingModel off
Selecting meshMotionModel noMeshMotion

CFDEMcoupling version: cfdem-2.9.0
, compatible to LIGGGHTS version: 3.1.0
, compatible to OF version and build: 2.3.x-commit-4d6f4a3115ff76ec4154c580eb041bc95ba4ec09

You are currently using:
OF version: 2.3.x
OF build: 2.3.x-6f86ad2099b6
CFDEM build: 64f5

If BC are important, please provide volScalarFields -imp/expParticleForces-
ignoring ddt(voidfraction)
Selecting forceModel noDrag
nrForceSubModels()=1
Selecting forceSubModel ImEx

--> FOAM FATAL ERROR:

request for volScalarField rho from objectRegistry region0 failed
available objects of type volScalarField are

9
(
voidfractionNext
UsWeightField_
nu
scalarViscosity
voidfractionPrev
ddtVoidfraction
p
voidfraction
sSmoothField
)

From function objectRegistry::lookupObject(const word&) const
in file /home/bc/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::GeometricField const& Foam::objectRegistry::lookupObject >(Foam::word const&) const at ??:?
#3 Foam::forceSubModel::forceSubModel(Foam::dictionary const&, Foam::cfdemCloud&, Foam::forceModel&) at ??:?
#4 Foam::ImEx::ImEx(Foam::dictionary const&, Foam::cfdemCloud&, Foam::forceModel&) at ??:?
#5 Foam::forceSubModel::adddictionaryConstructorToTable::New(Foam::dictionary const&, Foam::cfdemCloud&, Foam::forceModel&) at ??:?
#6 Foam::forceSubModel::New(Foam::dictionary const&, Foam::cfdemCloud&, Foam::forceModel&, Foam::word) at ??:?
#7 Foam::forceModel::setForceSubModels(Foam::dictionary&) at ??:?
#8 Foam::noDrag::noDrag(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
#9 Foam::forceModel::adddictionaryConstructorToTable::New(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
#10 Foam::forceModel::New(Foam::dictionary const&, Foam::cfdemCloud&, Foam::word) at ??:?
#11 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:?
#12
at ??:?
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14
at ??:?
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Thank you for reading my reply.

I think that we need to be a well-OpenFOAM user... The reason is the missing of the velocity(or pressure, Temperature,....) information in '0' folder.
Usually, the '0' folder has the information of initial condition (written in OpenFOAM manual).
Finally, the example case is not the perfect form (maybe the only 'dummy' file is located in '0' folder and you have to complete the information), so we have to start filling those things!.
Thank you for reading my reply.
-BC-