ERROR: Invalid fix style (modify.cpp:677)

Submitted by vipul.shrivastava on Mon, 09/24/2012 - 14:39

Hello,
Can someone please help me out with this error mentioned below. I have run the same script on LIGGGHTS successfully but it's showing error in the CFDEM solverpiso. Even pair style i changed to gran/hooke/history, it wasn't working with gran/hertz/history. Thanks in advance.

// run_parallel_cfdemSolverPiso_backward_step_CFDDEM //

/home/su/CFDEM/su-PUBLIC-2.1.1/run/backward_step/CFD

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : cfdemSolverPiso -parallel
Date : Sep 24 2012
Time : 17:48:57
Host : "openfoam"
PID : 12891
Case : /home/su/CFDEM/su-PUBLIC-2.1.1/run/backward_step/CFD
nProcs : 2
Slaves :
1
(
"openfoam.12892"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
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 kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}

Reading g
Selecting locateModel standard
Selecting dataExchangeModel twoWayMPI
Starting up LIGGGHTS for first time execution
Executing input script '../DEM/in.tetmesh'
LIGGGHTS (Version LIGGGHTS-PUBLIC 2.1, compiled 2012-09-06-12:30:11 by su based on LAMMPS 20 Apr 2012)

communicate single vel yes

units si

region reg block -0.098 0.196 0 0.0588 0 0.01 units box

#region reg block 0 1 0 1 0 1 units box

#region reg block -8.99584406 9.0 -15.5 11.600 -62 367 units box

#region reg block -0.0095 0.0095 -0.18 0.14 -0.70 4.00 units box

#region reg block -1 1 -1 1 -0.05 2.5 units box

create_box 1 reg

Created orthogonal box = (-0.098 0 0) to (0.196 0.0588 0.01)
2 by 1 by 1 MPI processor grid

neighbor 1 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.9

fix m4 all property/global coefficientFriction peratomtypepair 1 0.05

fix m5 all property/global characteristicVelocity scalar 2.

fix m6 all property/global cohesionEnergyDensity peratomtypepair 1 500

#New pair style

pair_style gran/hooke/history #Hertzian without cohesion

pair_coeff * *

timestep 0.00001

fix gravi all gravity 9.81 vector 0.0 -1.0 0.0

#import triangular mesh from cad; the syntax is as follows:

#"fix customname all mesh/gran stl_filename wall_type scale_factor x_off y_off z_off"

fix cad all mesh/gran stepduct.stl 1 1 0. 0. 0. 0. 0. 0.

ERROR: Invalid fix style (modify.cpp:677)

alice's picture

alice | Mon, 09/24/2012 - 15:18

Hello Vipul,
the error occurs when an stl-file is incorporated. Try using the command mesh/surface instead of mesh/gran (see doc for details), some syntax changes have taken place recently...
Cheers,
Alice

vipul.shrivastava | Tue, 09/25/2012 - 13:04

Hello Alice,
I checked the documents and tried mesh/surface command but nw getting the error mentioned below. I compiled the CFDEM and LIGGGHTS again but still it's showing the same error. I don't think there is anything wrong with .stl files. Help would be highly appreciated.

..............................
..............................
fix cad all mesh/surface file stepduct.stl type1

ERROR on proc 0: Cannot open mesh file stepduct.stl (input_mesh_tri.cpp:74)
--------------------------------------------------------------------------
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 5948 on
node openfoam exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
------------------------------------------------------------------------

alice's picture

alice | Tue, 09/25/2012 - 14:40

Hello,
please be aware that the liggghts input script is run from the CFD-folder, so in case your stl-file is located anywhere else (e.g., in the DEM folder) you have to add the relative path (e.g., ../DEM/mystl.stl) to the file name so it can be found.
Cheers,
Alice