LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

Cluster install (Raijin)

Submitted by mardussa on Thu, 11/09/2017 - 00:49

Hi everyone,

I am trying to build LIGGGHTS 3.7 on the Raijin cluster at NCI in Aus. LIGGGHTS-3.6 is already installed on the cluster, but I am aiming to install and use CFDEM as I have on the Uni cluster (which is being decommissioned next month, hence moving to Raijin).

The cluster already has VTK 6.3 installed, however I have decided initially to compile without VTK by altering the Makefiles I use, removing -DLAMMPS_VTK and the VTK library and include locations.

Having problem during running a cohesive force model

Submitted by mmkamyabi on Tue, 11/07/2017 - 14:25

Hello everybody,
I'm trying to run an example of cohesive force model of washino/capillary/viscous. When I run this script (as below) the memory usage of the computer is increased illogically so that the run is dead. I decreased the amount of particles but it doesn't solve the problem. May you please help me to know what's the reason and what I can do to fix it?
Thank you

the input script is:

#Contact model example

atom_style granular
atom_modify map array

boundary m m m
newton off

communicate single vel yes

units si

Implementation of a time-variable value

Submitted by Kaikeyi on Mon, 11/06/2017 - 17:42

Hey,

I got a question about implementating a value in CFDEM Ligghts.

How do I implement a time-variable value, which shall be variable at every time step?

Is this a possible way to realize my question?
Csurface = fix_Csurface --> vector_atom
vector_atom = Csurface
fix_Csurface = NULL

Thanks for your helpful minds :)

rotate mesh with surface velocity

Submitted by Nils Borghoff on Mon, 11/06/2017 - 13:05

Hello everybody,

I want to simulate a belt conveyor which also rotates around a defined axis.
Normaly I used this command to give the surface of the belt velocity components

fix cv1 all mesh/surface/stress file belt_1.stl type 3 scale 0.001 surface_vel ${beltvX} ${beltvY} ${beltvZ}

Then I want ro rotate this mesh with following command

fix rotatecv1 all move/mesh mesh cv1 rotate origin 0.0 698.8 0.0 axis 0.0 0.0 1.0 period ${time1}

Error in running the cohesion force model

Submitted by mehdikamyabi on Sat, 11/04/2017 - 14:23

Dear all,
hello,
I'm new on LIGGGHTS so please forgive me if the question is stupid.
I was running the cohesive example of LIGGGHTS. When I changed the cohesive force as washino/capillary/viscous I received this error: "ERROR: Could not locate a fix/property storing value(s) for surfaceLiquidContentInitial as requested by cohesion_model washino/capillary/viscous. (../modify_liggghts.cpp:411)"

could you please help me to fix this?
Thank you in advance

the in.cohesion file is as:

#Contact model example

Fix wall for another pair_style is supported?

Submitted by fabregas04 on Thu, 11/02/2017 - 13:41

Hi,

I want to plant the mesh which have charge on it and interact with particles.
But by fix wall/gran, there is only pair_style for granular is supported.
How can I use another pair_style (e.g lj/coul from LAMMMPS) to fix the wall?

And another question is that when I put the wall by fix wall/gran, sometimes particles are stick to the mesh(wall).
Do you guys know why this happened?

IF LOOP IN LIGGGHTS

Submitted by AnjanaKittu on Thu, 11/02/2017 - 02:21

Hi,

I an new to LIGGGHTS and trying to figure out how to create an if loop. I looked at the documentation, but not sure how to proceed.
Here is my current code, not sure if this works.

label shearing1
variable xpos
variable strain_inc equal 0.0
variable max_height equal 0.0377

if ${xpos} < ${max_height}*0.01
#if ${xpos} < ${max_height}*.0025+${strain_inc} then

fix move all move/mesh cad1 linear 1e-04 0.0 0.0
fix move all move botwall linear 1e-04 0.0 0.0

Pages

Subscribe to RSS - LIGGGHTS® - User Forum