hello cfdem users,
I am trying to simulate gas-solid fluidized bed. Im learning cfd-dem coupling with the help of HopperEmpting testcase and changing the input scripts according to my case. I have packed the solid particles in the bed with the help of Liggghts and a liggghts.restart file was generated. Here's the coupling script that i used for DEM part of cfdem coupling.....
echo both
atom_style granular
atom_modify map array sort 0 0
boundary f f f
newton off
communicate single vel yes
units si
processors 1 1 2
read_restart ../DEM/liggghts.restart
neighbor 0.001 bin
neigh_modify delay 0
fix m1 all property/global youngsModulus peratomtype 5.e6 5.e8
fix m2 all property/global poissonsRatio peratomtype 0.45 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.9 0.9 0.9 0.9
fix m4 all property/global coefficientFriction peratomtypepair 2 0.5 0.5 0.5 0.5
pair_style gran/hertz/history 1 0 #Hertzian without cohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#frictional walls
fix xwalls all wall/gran/hertz/history 1 0 xplane 0 0.11 1
fix ywalls all wall/gran/hertz/history 1 0 yplane 0 0.11 1
fix zwalls all wall/gran/hertz/history 1 0 zplane 0 1 1
fix cfd all couple/cfd/force couple_every 100 mpi
fix integr nve_group nve/sphere
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
dump mydump all stl 1 post/dump.stl
dump dmp all custom 500 post/dump.pramod_restart id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
undump mydump
run 10000
but this is giving an error like this....
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 = 1
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 standard
Selecting dataExchangeModel twoWayMPI
Starting up LIGGGHTS for first time execution
Executing input script '../DEM/in.pramod_couple'
LIGGGHTS 1.4.5 based on lammps-10Mar10
Reading restart file ...
WARNING: Restart file used different # of processors
WARNING: Restart file used different 3d processor grid
orthogonal box = (0 0 0) to (0.11 0.11 0.3)
1 by 1 by 2 processor grid
20000 atoms
Resetting global state of Fix xwalls Style wall/gran/hertz/history from restart file info
Resetting per-atom state of Fix xwalls Style wall/gran/hertz/history from restart file info
Resetting global state of Fix ywalls Style wall/gran/hertz/history from restart file info
Resetting per-atom state of Fix ywalls Style wall/gran/hertz/history from restart file info
Resetting global state of Fix zwalls Style wall/gran/hertz/history from restart file info
Resetting per-atom state of Fix zwalls Style wall/gran/hertz/history from restart file info
INFO: nevery as specified in LIGGGHTS is overriden by calling external programResetting global state of Fix contacthistory Style contacthistory from restart file info
Resetting per-atom state of Fix contacthistory Style contacthistory from restart file info
--------------------------------------------------------------------------
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.
--------------------------------------------------------------------------
ERROR on proc 0: Cannot open dump file
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 3515 on
node ubuntu 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).
My question is that is there any problem with dump.stl file? As im new with with the cfdem, i dont know whether i have to create this dump.stl file or not. Is this the same file which was created when particle filling in the column was done??
Any reply regarding this would be appreciated.
thanks
cgoniva | Mon, 02/20/2012 - 09:03
"Cannot open dump file" most
"Cannot open dump file"
most probably the path you to your desired dump location is not correct, it has to defined relative to the location you start the coupling from.
please try: ../DEM/post/dump.stl
Cheers,
Chris
techblast333 | Mon, 02/20/2012 - 14:57
thanks for the reply sir but
thanks for the reply sir but i have already tried with ..DEM/post/dump.stl script also but the error is same. Is there any file missing? dump.stl file itself would be created or i have to create this file and which data will go to this file?
thanks.
techblast333 | Mon, 02/20/2012 - 20:17
sir, " cannot open dump
sir,
" cannot open dump file" error is not anymore and coupling between CFD and DEM is also taking place at the initial stages its giving errors:
Pair time (%) = 0.000333913 (30.2055)
Neigh time (%) = 0 (0)
Comm time (%) = 1.04801e-05 (0.948016)
Outpt time (%) = 4.84769e-05 (4.38518)
Other time (%) = 0.000712602 (64.4613)
Nlocal: 250 ave 500 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Nghost: 0 ave 0 max 0 min
Histogram: 2 0 0 0 0 0 0 0 0 0
Neighs: 813 ave 1626 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Total # of neighbors = 1626
Ave neighs/atom = 3.252
Neighbor list builds = 0
Dangerous builds = 0
Selecting voidFractionModel centre
centreVoidFraction constructor done
Selecting averagingModel dilute
Selecting regionModel allRegion
Selecting meshMotionModel noMeshMotion
cfdem version: cfdem-2.1.0-beta-2011-09-01-14:13
If BC are important, please provide volScalarFields -imp/expParticleForces-
Selecting forceModel DiFeliceDrag
Selecting momCoupleModel implicitCouple
Selecting liggghtsCommandModel execute ,provide dicts, numbered from 0 to n
liggghtsCommand run 100
firstCouplingStep = 1
lastCouplingStep = 1000000000
couplingStepInterval = 1
Starting time loop
Time = 0.001
Courant Number mean: 0 max: 0
- evolve()
Starting up LIGGGHTS
Executing command: 'run 100 '
run 100 Setting up run ...
Memory usage per processor = 10.9107 Mbytes
Step Atoms KinEng 1 Volume
100001 500 7.5576706e-08 8.5823724e-08 0.00363
100101 500 7.8238772e-08 8.8074865e-08 0.00363
Loop time of 0.13338 on 2 procs for 100 steps with 500 atoms
Pair time (%) = 0.0378126 (28.3495)
Neigh time (%) = 0 (0)
Comm time (%) = 0.000525416 (0.393924)
Outpt time (%) = 6.79287e-05 (0.0509287)
Other time (%) = 0.0949741 (71.2056)
Nlocal: 250 ave 500 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Nghost: 0 ave 0 max 0 min
Histogram: 2 0 0 0 0 0 0 0 0 0
Neighs: 813 ave 1626 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Total # of neighbors = 1626
Ave neighs/atom = 3.252
Neighbor list builds = 0
Dangerous builds = 0
LIGGGHTS finished
timeStepFraction() = 1
evolve done
update Ksl.internalField()
totaldragforceEuler calculus
totaldragforceEuler = 0
dv/dt =sum(ddt(voidfraction)) [0 0 -1 0 0 0 0] 0
[0] #0 Foam::error::printStack(Foam::Ostream&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #2 Uninterpreted:
[0] #3 void Foam::fvc::surfaceIntegrate >(Foam::Field >&, Foam::GeometricField, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[0] #4 Foam::tmp, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate >(Foam::GeometricField, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[0] #5 Foam::tmp, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate >(Foam::tmp, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[0] #6 Foam::fv::gaussDivScheme >::fvcDiv(Foam::GeometricField, Foam::fvPatchField, Foam::volMesh> const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[0] #7 Foam::tmp, Foam::Tensor >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div >(Foam::GeometricField, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
[0] #8 Foam::tmp, Foam::Tensor >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div >(Foam::GeometricField, Foam::fvPatchField, Foam::volMesh> const&) in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
[0] #9 Foam::incompressible::RASModels::kEpsilon::divDevReff(Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) const in "/home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
[0] #10
[0] in "/home/anoop/OpenFOAM/anoop-2.0.x/platforms/linuxGccDPOpt/bin/cfdemSolverPiso_shared"
[0] #11 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
[0] #12
[0] in "/home/anoop/OpenFOAM/anoop-2.0.x/platforms/linuxGccDPOpt/bin/cfdemSolverPiso_shared"
[ubuntu:04816] *** Process received signal ***
[ubuntu:04816] Signal: Floating point exception (8)
[ubuntu:04816] Signal code: (-6)
[ubuntu:04816] Failing at address: 0x12d0
[ubuntu:04816] [ 0] [0xc5a410]
[ubuntu:04816] [ 1] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x61) [0x8da7831]
[ubuntu:04816] [ 2] [0xc5a400]
[ubuntu:04816] [ 3] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam3fvc16surfaceIntegrateINS_6VectorIdEEEENS_3tmpINS_14GeometricFieldIT_NS_12fvPatchFieldENS_7volMeshEEEEERKNS5_IS6_NS_13fvsPatchFieldENS_11surfaceMeshEEE+0x898) [0x16e1b28]
[ubuntu:04816] [ 4] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam3fvc16surfaceIntegrateINS_6VectorIdEEEENS_3tmpINS_14GeometricFieldIT_NS_12fvPatchFieldENS_7volMeshEEEEERKNS4_INS5_IS6_NS_13fvsPatchFieldENS_11surfaceMeshEEEEE+0x38) [0x16e2018]
[ubuntu:04816] [ 5] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam2fv14gaussDivSchemeINS_6TensorIdEEE6fvcDivERKNS_14GeometricFieldIS3_NS_12fvPatchFieldENS_7volMeshEEE+0x75) [0x16e2a75]
[ubuntu:04816] [ 6] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam3fvc3divINS_6TensorIdEEEENS_3tmpINS_14GeometricFieldINS_12innerProductINS_6VectorIdEET_E4typeENS_12fvPatchFieldENS_7volMeshEEEEERKNS5_IS9_SC_SD_EERKNS_4wordE+0x6d) [0x324ecd]
[ubuntu:04816] [ 7] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam3fvc3divINS_6TensorIdEEEENS_3tmpINS_14GeometricFieldINS_12innerProductINS_6VectorIdEET_E4typeENS_12fvPatchFieldENS_7volMeshEEEEERKNS5_IS9_SC_SD_EE+0xc4) [0x325084]
[ubuntu:04816] [ 8] /home/anoop/OpenFOAM/OpenFOAM-2.0.x/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so(_ZNK4Foam14incompressible9RASModels8kEpsilon10divDevReffERNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEE+0xe4) [0x336984]
[ubuntu:04816] [ 9] cfdemSolverPiso_shared() [0x80605d3]
[ubuntu:04816] [10] /lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe6) [0x2042bd6]
[ubuntu:04816] [11] cfdemSolverPiso_shared() [0x805c761]
[ubuntu:04816] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 4816 on node ubuntu exited on signal 8 (Floating point exception)
I'm not able to find why this error is coming. Have i used any wrong model or solver or this error is because of wrong volume fraction or time step??
Anoop