I just started out with LIGGGHTS and after running some of the standard tutorials, I found this Bin Flow turorial, which should work but is not working for me at the moment
Included is my case folder and when running this I get the following error:
LIGGGHTS (Version LIGGGHTS-PUBLIC 3.3.0, compiled 2016-01-11-16:35:03 by bnaar, git commit 9b1fc57b77c483deaf2c0a6c69241c5d4f91e3c2 based on LAMMPS 23 Nov 2013)
Created orthogonal box = (-0.138 -0.138 -0.0045) to (0.138 0.138 0.43)
1 by 1 by 1 MPI processor grid
Fix particledistribution/discrete (id pdd): distribution based on mass%:
pts: d=3.000000e-03 (max. bounding sphere) mass%=100.000000%
Fix particledistribution/discrete (id pdd): distribution based on number%:
pts: d=3.000000e-03 (max. bounding sphere) number%=100.000000%
ERROR on proc 0: Cannot open mesh file hopper.stl (../input_mesh_tri.cpp:100)
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
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.
--------------------------------------------------------------------------
And as I am also not that skillful with Linux, I now do not know what the problem is as I believe that the tutorial is set up correctly.
All tips are welcome
| Attachment | Size |
|---|---|
 liggghts_tutorial.pdf | 1.02 MB |
 binflow.tar_.gz | 9.56 KB |
hunger | Thu, 01/14/2016 - 08:35
Did you consider the hint in
Did you consider the hint in the Download section of the tutorial?
"Important note: This tutorial does not reflect the syntax changes from LIGGGHTS 2.X to LIGGGHTS 3.X. Please refer to doc/liggghts_3.X_tutorial.html in your local installation directory for details on the syntax changes"
If you did not change the syntax, this could be the reason for your problem.
bnaar | Thu, 01/14/2016 - 18:14
Your hint indeed let me
Your hint indeed let me change the wall/gran and pair_style commands (so thank you!), but still the same error is appearing that Liggghts is not able to load in the mesh
I did change the mesh command now to:
fix cad1 all mesh/surface file mesh/hopper.stl type 2 scale 0.001
as the hopper.stl file was inside a subfolder mesh.
But the same error is appearing, hence, would you know what the problem could be?
Best regards,
Boris
bnaar | Thu, 01/14/2016 - 18:28
Okay, I already solved it,
Okay, I already solved it, there was a small typo somewhere, my bet!
Thank you.
Best,
Boris
VOny | Fri, 03/24/2017 - 23:45
Can you specify the typo you
Can you specify the typo you found? thanks
ThoriumDeLorean | Mon, 05/28/2018 - 21:00
Syntax changes and Typos
To run the Bin Flow tutorial with LIGGGHTS(R)-PUBLIC 3.X you must make the following syntax changes:
OLD 2.X Syntax:
# Use the imported mesh as granular wall
fix geometry all wall/gran/hertz/history mesh n_meshes 1 meshes cad1# Create stopper for funnel
fix stopper all wall/gran/hertz/history primitive type 2 zplane 0.0# Define the physics
pair_style gran/hertz/historyNEW 3.X Syntax:
# Use the imported mesh as granular wall
fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1# Create stopper for funnel
fix stopper all wall/gran model hertz tangential history primitive type 2 zplane 0.0# Define the physics
pair_style gran model hertz tangential historyThere are also two typos in the in.BinFlow script:
"
fox integrate all nve/sphere" --> "fix integrate all nve/sphere""
thermo_modigy norm no lost ignore" --> "thermo_modify norm no lost ignore"Hope this helps!
Dheeraj Minglani | Tue, 05/29/2018 - 08:26
Even after fixing getting error
fix pts all particletemplate/sphere 1 atom_type 1 density constant 1000 &
radius constant 0.0015
fix pdd all particledistribution/discrete 63243 1 pts 1.0
region factory cylinder z 0 0 0.10 0.4 0.41 units box
fix ins all insert/rate/region seed 123481 distributiontemplate pdd &
nparticles 150000 particlerate 150000 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -1.0 region factory ntry_mc 10000
# Import mesh from cad:
fix cad1 all mesh/surface file mesh/hopper.stl type 2 scale 0.001
# Use the imported mesh as granular wall
fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1
# Create stopper for funnel
fix stopper all wall/gran model hertz tangential primitive history type 2 zplane 0.0
# Define the physics
pair_style gran hertz tangential history
pair_coeff * *
### Detailed settings
# Integrator
fix integrate all nve/sphere
# Gravity
fix grav all gravity 9.81 vector 0.0 0.0 -1.0
# Time step
timestep 0.00001
# Thermodynamic output settings
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify norm no lost ignore
ERROR: Detected problem with model 'hertz' (and possibly subsequent arguments).
ERROR: The user must ensure that the contact model for this fix exists in the list of variants.
ERROR: For the list of variants see the 'style_tangential_model.h' and 'style_normal_model.h' file in the src directory of your LIGGGHTS installation.
ERROR: Fix wall/gran (id stopper): unknown contact model (/build/liggghts-QCz4in/liggghts-3.3.1+repack1/src/fix_wall_gran.cpp:168)
jsams | Wed, 05/30/2018 - 02:03
The primitive wallstyle
The primitive wallstyle command needs to be appended after the contact model definition. I think you might have accidentally mixed it up for the stopper.
Instead of the following:
fix stopper all wall/gran model hertz tangential primitive history type 2 zplane 0.0
use:
fix stopper all wall/gran model hertz tangential history primitive type 2 zplane 0.0